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50TH  ANNIVERSARY  COMMEMORATED 


Appropriately,  the  closing  date  of  the  Symposium  marked  the 
50th  anniversary  of  a  significant  milestone  in  aviation.  On  8  May 
1919,  three  naval  aircraft  left  Long  Island  to  attempt  the  first 
crossing  of  the  Atlantic.  Of  the  three  Curtiss  flying  boats  that 
started  that  historic  3,925-nautical-mile  flight,  the  NC-4  (shown 
on  cover),  commanded  by  LCDR  A.  C.  Read,  was  successful, 
making  the  first  Atlantic  crossing  via  Newfoundland,  the  Azores, 
and  Portugal,  finally  arriving  at  Plymouth,  England. 

The  Secretary  of  the  Navy  designated  May  1969  a 
commemorative  period,  and  it  was  especially  appropriate  that  the 
Symposium  salute  the  aeronautical  pioneers  who  made  history  in 
May  1919. 
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FOREWORD 


These  Proceedings ,  published  in  five  volumes,  comprise  the  49  papers  presented 
at  the  Eighth  Navy  Symposium  on  Aeroballistics  held  at  the  Naval  Weapons  Center 
Corona  Laboratories,  Corona,  Calif.,  6,  7,  and  8  May  1969. 

This  symposium  was  the  eighth  in  a  series  begun  in  1950  under  the 
sponsorship  of  the  then  Bureau  of  Ordnance  Committee  on  Aeroballistics,  and 
currently  conducted  by  the  Naval  Aeroballistics  Advisory  Committee  as  sponsoring 
committee  for  the  Naval  Air  Systems  Command  and  the  Naval  Ordnance  Systems 
Command.  The  continuing  purpose  of  the  symposiums  has  been  to  disseminate  the 
results  of  aeroballistics  research  and  to  bring  the  research  findings  of  industry,  the 
universities,  and  government  laboratories  to  bear  upon  the  Navy’s  aeroballistics 
research  and  development  programs. 

Over  200  research  scientists  representing  more  than  72  organizations  attended 
this  eighth  symposium.  Sessions  1  and  2  covered  the  subjects  of  heat  transfer  and 
aerophysics,  nozzles  and  jet  effects;  Sessions  3  and  4  were  concerned  with 
aerodynamics  and  missile  stability;  and  Session  5  dealt  with  structures  and 
aeroelasticity,  and  external  carriage  and  store  separation. 

The  papers  in  these  Proceedings  have  been  reproduced  in  facsimile.  They 
appear  in  the  order  of  presentation  except  that  all  classified  papers  have  been  taken 
out  of  sequence  and  grouped  together  as  Volume  5,  a  confidential  volume.  Volumes 
1  through  4  are  unclassified.  This  is  Volume  1. 

Requests  for  or  comments  on  individual  papers  should  be  addressed  to  the 
respective  authors. 


RAY  W.  VAN  AKEN 
General  Chairman 
Symposium  Committee 


Published  by  the  Publishing  Division  of  the  Technical  Information  Department,  NWC;  first 
printing,  June  1969,  250  copies. 
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WELCOME 


To  conference  attendees: 

On  behalf  of  the  military  and  civilian  personnel  of  the  Naval  Weapons  Center  we 
extend  a  hearty  “welcome  aboard”.  It  is  an  honor  for  the  Naval  Weapons  Center  to 
be  the  host  organization  for  this  Symposium,  which  is  sponsored  by  the  Naval 
Aeroballistics  Advisory  Committee  on  behalf  of  the  Naval  Air  Systems  Command 
and  Naval  Ordnance  Systems  Command.  We  have  been  pleased  by  the  interest  shown 
in  the  meeting  and  hope  that  it  will  be  a  profitable  one. 


Dr.  Thomas  S.  Amlie 
Technical  Director,  NWC 


CAPT  M.  R.  Etheridge,  USN 
Commander,  NWC 
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The  charter  of  the  Naval  Aeroballistics  Advisory  Committee  (NAAC)  assigns  it  the 
responsibility  to  review,  at  least  annually,  the  current  and  projected  Navy  research 
and  development  effort  in  aeroballistics  and  to  make  recommendations  to  the 
joint-sponsoring  commands,  NAVAIR  and  NAVORD,  for  the  support  of  important 
and/or  critical  work  areas. 

It  is  encouraging  to  note  that  almost  half  of  the  papers  being  presented  at  this 
symposium  result  from  studies  prompted  by  your  recommendations  of  recent  years. 
The  wide  variety  of  the  subject  matter  attests  to  the  broad  interests  of  the  two 
command  sponsors  and  is  a  testimonial  to  the  broad  base  of  knowledge  which  our 
joint  aero  research  and  development  programs  are  generating.  The  wide 
representation  you  draw  from  Air  Force,  Army,  universities  and  industry  is  further 
indication  of  the  cross  fertilization  which  exists  and  in  large  part  must  be  credited 
to  the  aeroballistics  committee  activity. 

The  flexibility  of  the  NAAC  organization  in  addressing  Navy’s  needs  has  again 
been  demonstrated  by  the  organization  of  a  panel  on  separation  of  stores  from 
aircraft.  You  will  note  heavy  emphasis  on  store  carriage  and  separation  problems  in 
this  symposium  agenda.  With  the  assistance  of  this  panel  we  are  embarked  upon  a 
strong  program  seeking  solutions  to  these  problems,  striving  for  significant 
improvements  in  weapon  delivery,  safety,  and  accuracy.  Such  coordinated  effort 
brought  to  bear  on  this  complex  aircraft/missile  interface  problem  will  produce 
results  and  drive  home  recognition  that  aircraft  weapon  delivery  systems  are  systems, 
not  separate  development  programs  of  aircraft,  launchers,  and  weapons,  simply  glued 
together  in  production  as  a  package. 

While  it  is  necessary  to  find  solutions  to  current  problems,  it  is  equally  important 
that  we  maintain  long  range  research  and  exploratory  development  programs,  looking 
well  ahead  for  the  advanced  technology  which  will  be  needed  in  future  defense 
developments.  For  example,  the  series  of  papers  concerning  jet-interaction  control 
fall  in  this  category  and  point  to  attaining  technical  feasibility  of  a  finless  missile, 
with  its  promise  of  tactical  and  logistical  advantages. 

Much  of  the  remainder  of  the  program  illustrates  the  increasing  improvements  in 
experimental  techniques  and  the  ever  increasing  utilization  of  high  speed  computer 
technology.  These  combined  advances  should  go  a  long  way  in  improving  the 
analytical  and  theoretical  background  to  understand  the  behavior  of  fluid  systems, 
where  presently  we  depend  so  heavily  on  an  empirical  approach. 

It  may  be  desirable  to  point  out  that  in  the  May  1966  reorganization  of  the 
Navy  Department,  within  the  Navy  Material  Command  management  of  research  and 
the  early  stages  of  development  has  been  placed  on  an  organizational  par  with  the 
hardware  engineering  development  effort.  The  technology  funding  and  budget  items 
no  longer  fall  easy  prey  to  every  disaster  or  over-run  in  current  acquisition 
programs. 
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There  is  a  busy  and  attractive  future  ahead  in  aerodynamic  science  and 
technology.  As  scientific  fields  and  disciplines  go,  this  one,  despite  remarkable 
progress  of  the  past  50  years,  is  still  in  its  infancy.  There  is  room  for  youth,  new 
uninhibited  approaches,  innovations,  and  I’m  sure  there  are  great  surprises  of 
discovery  and  invention  ahead.  As  scientific  men  might  be  measured,  this  science  is 
really  just  leaving  the  first  generation  the  second  generation  is  now  well 
established— and  the  future  belongs  to  those  now  maturing  from  their  apprenticeship. 
I  suspect  we  are  just  beginning  the  period  where  some  of  the  science  is  old  enough 
to  be  forgotten  and  that  needs  careful  watch  in  this  arrogant  age. 

This  is  a  complex  physical  field— the  mathematics  is  tortuous— the  equations 
massive  and  complicated,  because  of  very  real  first,  second  and  third  order 
interactions.  In  the  past  this  has  forced  empiricism  to  a  high  degree— skillful  guess 
work  and  lucky  estimations  have  been  necessary— but  with  the  mathematical 
flexibility  and  near  omnipotence  of  the  massive  digital  computers  now  available,  the 
theory,  the  mathematics,  and  the  tedium  of  nonlinear  coefficients  can  be  properly 
assailed.  This  will  open  new  horizons  in  every  direction. 

Symposia  are  the  hallmark  of  our  scientific  and  engineering  community.  We  must 
love  them.  They  proliferate  at  an  astounding  rate.  Often  they  seem  to  be  drawn  to 
lovely  location  spots,  “The  high  rent  districts  attract  them.”  Now  managers,  as 
distinct  from  the  scientist/engineer  class,  tend  to  live  in  a  world  of 
cost-consciousness.  Time  and  again  I  note  instructions,  policies  and  the  like  setting 
restrictions,  limitations  on  attendance,  etc. 

Yet  we  need  this  communication.  The  body  of  knowledge  is  rapidly  enhanced  by 
communication  and  the  sharing  it  engenders.  As  we  get  better  educated  we  need 
more  and  more  of  this  communication— and  despite  the  miracles  of  printing,  data 
transmission,  telephone,  radio,  and  even  the  modern  copying  machines,  face  to  face 
human  communication  hasn’t  been  beat. 

So  I  exhort  each  of  you,  to  give  something,  to  take  something,  or  both,  from 
your  friends  at  this  meeting.  Make  it  profitable  for  yourself  and  for  someone  else. 
And  that  in  itself  advances  the  so-called  “state  of  the  art.” 

In  conclusion  let  me  suggest  the  strength  of  the  program  for  this  symposium  is 
an  excellent  indication  that  the  Naval  Aeroballistics  Advisory  Committee  remains  a 
strong  and  responsive  organization,  and  will  continue  to  exercise  a  great  influence  on 
the  aeroballistics  research  and  development  programs  of  the  Navy  and  the  nation. 
Both  NAVAIR  and  NAVORD,  integrated  organizationally  under  the  Chief  of  Naval 
Material,  intend  to  maintain  support  and  attain  progress  in  this  vital  scientific  area 
and  discipline. 
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INTRODUCTION 


The  base  and  lee  side  flows  of  slender  bodies  represent  funda¬ 
mental  fluid  mechanics  problems  which  derive  their  basic  features  from 
the  fact  that  flow  separation  exists.  The  necessity  of  studying  these 
regions  in  detail  has  been  dictated  by  the  effect  of  their  local  flow 
properties  on  guidance,  maneuverability,  communication  and  identifi¬ 
cation  of  high  speed  vehicles.  Design  of  vehicles  or  components  that 
will  function  in  the  separated  flow  regimes  is  only  possible  if  a 
reasonably  complete  picture  of  the  flow  field  is  available.  Due  to  the 
extreme  complexity  of  the  flow,  which  at  the  present  is  not  completely 
understood  physically,  there  are  formidable  difficulties  in  the 
analytical  solution  of  the  problem. 

A  configuration  that  demands  considerable  attention  is  the  slender 
body  at  angle  attack.  This  problem  is  of  extreme  practical  importance 
since  a  maneuverable  missile  will  not  be  at  zero  attitude  throughout 
its  entire  flight  trajectory.  In  fact,  in  such  vehicles,  the  angle  of 
attack  can  be  quite  large. 

Fig.  1  illustrates  a  tentative  picture  of  the  flow  field  around 
a  slender  body  at  angle  of  attack  in  super/hypersonic  speeds.  Experi¬ 
mental  measurements  in  Ref.  1  indicate  that  some  features  of  this 
model  might  be  correct,  others  (^uch  as  the  mixing  of  the  lee  side  flow 
with  the  base  flow  and  the  shape  and  characteristics  of  the  base  flow) , 
however,  are  based  on  intuition.  The  flow  field  shows  the  presence  of 
shock  waves  (besides  the  primary  shock  there  are  imbedded  secondary 
shocks  on  the  lee  side).  These  are  associated  with  strong  pressure 
gradients  which  give  rise  to  interactions  between  boundary  layers  and 
shock  waves.  These  interactions  lead  to  flow  separation.  The  incidence 
of  the  slender  body  to  the  primary  flow  causes  vortex  generation  which 
results  in  rolled  up  vortex  sheets  in  the  flow.  Recently  there  were 
indications  that  these  vortex  sheets  might  cause  the  reattachment  of 
the  flow  at  the  leeward  meridian. 

It  is  obvious  from  the  aforesaid  that  analysis  dealing  with  first 
principles  only  will  not  yield  results,  and  a  great  deal  of  empiricism 
must  be  employed  to  obtain  engineering  solutions  which  can  be  used  to 
predict  pressure  and  heat  transfer  in  the  leeward  and  base  regions  of 
a  general  body  of  revolution. 

A  comprehensive  review  of  the  available  base  flow  theories  in 
supersonic  flow  was  recently  presented  in  Ref.  2,  where  some  175  papers 
and  articles  were  reviewed.  Base  flow  theories  may  be  divided  into 
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four  main  groups,  namely  semiempirical  theories,  those  based  on  the 
Chapman-Korst  model,  integral  methods,  and  multimethod  base  flow 
theories  which  attempt  to  take  the  dynamics  of  the  recirculating  flow 
into  account. 

The  overwhelming  majority  of  the  experimental  information  reported 
in  the  literature  (in  the  form  of  pressure  and  heat  transfer  data  in 
the  base  region)  deals  with  either  two-dimensional  or  zero  angle  of 
attack  axisymmetric  geometries.  A  comprehensive  summary  of  the  experi¬ 
mental  reports  issued  over  the  past  several  years  on  this  subject  has 
been  presented  in  Ref.  3,  and  includes  near  wake  data  for  a  variety 
of  vehicle  configurations,  free  stream  Mach  numbers,  Reynolds  numbers, 
wall  temperature  ratios,  etc. 

The  the  past  few  years  a  few  investigators  obtained  some  data  in 
the  base  region  of  slender  cones  at  angles  of  attack  up  to  10°. 
Schlesinger  and  Martelluci  (Ref.  4)  tested  a  10°  half  angle  cone  at 
M  =6.0  and  at  a  10^  angle  of  attack.  The  free  stream  Reynolds  number 
is  sufficiently  large  so  that  fully  turbulent  flow  was  achieved  both  on 
the  cone  surface  and  in  the  near  wake.  Significant  changes  in  the  flow 
field  were  observed  due  to  the  angle  of  attack  of  the  body  compared  to 
the  zero  angle  of  attack  configuration. 

Schimdt  and  Cresci  (Ref.  5)  examined  the  flow  characteristics  in 
the  near  wake,  a  10°  half  angle  cone  at  M  =  7.7  and  at  10°  angle  of 
attack  for  laminar  flow.  They  obtained  radial  variations  of  pitot 
and  static  pressures  at  selected  axial  location  in  the  near  wake. 

Their  measurements  indicated  that  the  angle  of  attack  produced  a  region 
on  the  leeward  side  of  the  cone  surface  wherein  the  boundary  layer 
which  was  originally  laminar  on  the  windward  side  became  transitional 
or  turbulent  at  lower  Reynolds  numbers  than  expected  from  the  axisym¬ 
metric  flow  case.  This  effected  the  mixing  processes  and  therefore 
the  behavior  of  the  local  flow  conditions  in  the  wake.  In  addition, 
it  appeared  that  for  the  angle  of  attack  configuration  there  was  a 
large  scale  mixing  caused  by  the  vertical  inviscid  flow  above  the  cone 
surface.  Pitot  pressure  profiles  showed  that  the  symmetry  axis  in  the 
viscous  core  was  displaced  toward  the  leeward  side  of  the  cone  by  about 
one  tenth  of  the  base  diameter.  The  stagnation  pressure  ratio  increased 
by  a  factor  of  five  above  the  laminar  axisymmetric  flow  conditions. 

Much  of  the  published  experimental  information  concerning  base 
flow  has  been  clouded  by  the  uncertainty  introduced  by  model  support 
interference.  The  sting  model  support  system,  in  common  usage  in  wind 
tunnel  measurements,  is  bound  to  distort  the  flow  field  to  some  degree; 
consequently,  the  reliability  of  the  resulting  data  might  be  question¬ 
able. 


It  is  the  objective  of  the  experimental  program  currently  in 
progress  in  the  Aerodynamics  Laboratory  of  NSRDC  to  obtain  reliable, 
interference-free  data  of  separated  flow  regimes  surrounding  slender 
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bodies  at  angles  of  attack  for  a  wide  range  of  speeds  (subsonic  to 
hypersonic)  and  angles  of  attack  (up  to  60°) .  These  data  would 
facilitate  engineering  solutions  which  can  be  used  to  predict  fluid 
flow  properties  in  the  leeward  and  near  wake  regions.  Furthermore, 
it  is  hoped  that  the  final  outcome  of  the  program  would  contribute  to 
the  basic  understanding  of  the  physical  mechanisms  which  govern  the 
fluid  dynamic  behavior  in  the  separated  flow  regimes. 

To  date,  the  major  emphasis  has  been  placed  on  establishing  the 
basic  parameters  of  the  problem,  developing  and  testing  methods  of 
model  support  in  the  wind  tunnel  which  would  have  only  a  small  effect 
on  the  separated  flow,  and  evaluating  suitable  instrumentation  for  the 
final  experiments.  To  this  end  a  series  of  preliminary  experiments 
were  carried  out  at  M  =  6.22  and  9.89,  and  unit  Reynolds  number  of 
1  x  10^/ft.,  for  angles  of  attack  up  to  60°  using  sting  mounted,  band 
supported,  and  free  drop  models.  In  these  experiments  surface  and  base 
heat  transfer,  base  pressure  and  wake  survey  tests  were  conducted 
together  with  oil  flow,  schlieren,  and  high-speed  motion  picture 
visualization  techniques.  The  results  of  these  tests  are  outlined  in 
the  present  report. 

A  considerable  amount  of  work  was  done  on  the  development  of  a 
direction  insensitive  miniature  total  pressure  probe,  to  be  used  in 
subsequent  surveys  of  separated  flow  regions,  and  on  the  evaluation  of 
a  multi-channel  pressure  telemetry  system  installed  in  the  base  of  a 
small  model.  The  results  of  this  work  are  reported  in  Ref.  20. 


TEST  APPARATUS  AND  EXPERIMENTAL  TECHNIQUES 

WIND  TUNNEL  FACILITY  AND  DATA  ACQUISITION  SYSTEM 

All  of  the  experiments  described  herein  were  conducted  in  the 
NSRDC  hypersonic  facility.  High  speed  acquisition  of  all  test  data  is 
provided  by  a  high-speed  fifty-channel  Beckman  Model  210  data  acquisi¬ 
tion  system.  Details  of  the  facility  can  be  found  in  Ref.  6. 

MODELS  AND  MODEL  MOUNTING  SYSTEMS 

Several  models  with  identical  outside  dimensions  were  used  in  the 
test  program.  The  basic  model  configuration  consisted  of  a  9°  half 
angle,  six  inches  long,  sharp  cone.  All  the  models  were  fabricated 
from  type  416  stainless  steel  with  mirror  surface  finishes  and  geometric 
tolerances  of  ±0.001  inch. 

Sting  Mounted  Model 

This  configuration  was  used  to  obtain  heat  transfer  and  base 
pressure  data.  Some  wake  and  lee  side  surveys  and  oil  flow  information 
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relating  to  separation  lines  were  also  performed  with  this  model. 

Four  Chromel-Alumel  thermocouples  were  installed  on  the  surface  of 
the  9°  cone  three  inches  from  the  tip  and  90°  apart  to  measure  surface 
heat  transfer.  In  addition,  one  chromel-alumel  thermocouple  was  mounted 
in  the  base  of  the  model. 

Base  pressures  were  measured  at  two  locations,  180°  apart  at  0.612 
inches  from  the  model  axis.  The  model  support  assembly  was  attached  to 
the  sector  at  the  base  of  the  sting.  Fig.  2  shows  the  sting  mounted 
cone  assembly  in  the  test  section. 


Band  Supported  Models 

In  order  to  accommodate  a  model  support  with  minimal  effect  on  the 
separated  flow  regimes,  a  horizontal  carriage  system  was  designed  and 
installed  in  the  hypersonic  tunnel. 

The  model  support  itself  consisted  of  two  7-inches  long,  0.029 
inch  thick,  and  0.50  inches  wide  bands  constructed  of  a  special  steel 
alloy  (uniloy  A-286)  which  did  not  lose  its  tension  strength  at  high 
temperatures  (higher  than  90,000  psi  at  1500°  F) .  Fig.  3  shows  a 
photograph  of  the  band  support  system. 


PRESSURE  INSTRUMENTATION 


Total  Pressure  Transducers 

Two  calibrated  Statham  pressure  transducers  0-2  and  0-5  psia  were 
used  for  the  Mach  and  various  wake  surveys.  These  were  accurate  and 
repeatable  to  within  ±0.57o  full  scale. 


Base  Pressure  Transducers 

Because  base  pressure  values  at  hypersonic  speeds  are  quite  low, 
conventional  pressure  transducers  would  be  inadequate.  Therefore  two 
low  pressure,  Datametrics  type  1014  Electronic  Manometer  and  type  511-3 
Barocel  pressure  sensor  systems  were  used  for  these  measurements. 

These  systems  were  capable  of  measuring  pressures  between  0  and  1  psia 
on  seven  consecutive  scales  from  0.001  to  1.0  psi  full  scale  with 
accuracy  and  linearity  of  ±0.1%  of  full  scale  in  each  range.  The 
National  Bureau  of  Standards  calibrated  both  systems  to  ±0.05%  full 
scale  accuracy. 
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interference  in  supersonic  speeds.  It  has  been  confirmed  that  the 
base  pressure  is  strongly  influenced  by  sting  interference.  Based  on 
Refs.  8  to  12  the  measured  critical  sting- to-base  diameter  and  sting 
length-to-base  diameter  ratios  for  laminar  flow  were  plotted  and  faired 
in  Figs.  5  and  6  as  functions  of  the  Mach  number.  In  the  present 
experiments  structural  considerations  limited  the  maximum  sting  length- 
to-base  diameter  ratio  to  9.3,  and  the  minimum  sting  diameter-to-base 
diameter  ratio  to  0.32.  These  were  used  in  both  tested  Mach  numbers. 

At  M  =  6.22  the  sting  geometry  is  not  too  far  from  the  extrapolated 
values  of  about  9.7  and  0.20.  At  M  =  9.89  the  data  of  Figs.  5  and  6 
will  not  yield  meaningful  results,  and  therefore  the  zero  angle  of 
attack  base  pressure  measurement  is  compared  with  free  flight  model 
tests  under  similar  conditions.  This  comparison  with  references  shows 
that  the  measured  data  in  the  present  experiments  are  in  fair  agreement 
with  the  reported  data  values. 

There  is  no  information  in  the  literature  about  sting  effects  at 
angle  of  attack  and  therefore  the  data  presented  here  should  only  be 
regarded  as  preliminary. 


Pressure-Time  Response 

Since  the  measured  base  pressure  values  were  quite  low,  considerable 
time  was  required  to  reach  steady  state.  Computations,  based  on  Ref.  13, 
showed  that  in  the  present  apparatus,  depending  on  the  measured 
pressure,  the  time  response  ranges  from  20  to  45  seconds  for  1  percent 
accuracy.  Subsequent  measurements  showed  these  time  estimates  to  be 
within  ±87,  of  experimental  value. 


Pressure-Error  Due  to  High  Temperature 

Model  surface  pressure  measurements  leeward  side,  and  base 
pressure  readings  are  subject  to  error  due  to  a  gaseous  flow  phenomenon 
which  occurs  at  the  boundary  of  a  surface  and  a  gas  at  low  pressure 
and  high  temperature.  Ref.  14  presents  a  computational  procedure  for 
the  error  estimation.  According  to  this  calculation,  the  maximum 
error  at  M  =  9.89  is  around  2.57,,  the  average  error  however  is  less 
than  17>.  The  data  are  not  corrected  for  this  error. 

ACCURACY  ESTIMATES 

Due  to  slight  instrumentation  inaccuracy  certain  errors  exist. 

It  is  estimated  that  quantities  which  are  computed  on  the  basis  of 
pitot  pressure  measurements,  such  as  Mach  number  and  wake  surveys, 
have  a  probable  error  of  ±1%.  The  individual  temperature  measurements 
have  a  probable  error  of  ±27,.  The  error  of  base  pressure  measurements 
below  <y  =  20°  is  tentatively  estimated  at  ±257,,  this  increases  to  ±507, 
at  a  =  40°.  Beyond  this  pitch  angle  the  sting  effects  completely 
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RESULTS  AND  DISCUSSION 

SURFACE  TEMPERATURE  DISTRIBUTION 

A  sample  of  the  measured  temperature-time  data  is  presented  in 
Fig.  7.  The  time  history  shows  that  the  initial  transient  after 
model  injection  die  out  in  about  4  seconds  and  thereafter  the  slopes 
are  nearly  constant.  The  highest  heating  rate  occurs  at  the  most 
windward  side  which  is  signified  by  the  highest  temperature-time  slope. 
The  temperature- time  rate  90°  from  the  windward  meridian  is  equal  on 
each  side  throughout  the  angle  of  attack  range  showing  that  the  flow  is 
symmetrical  about  the  meridian  plane.  The  heating  rate  on  the  lee  side 
apparently  is  higher  than  in  the  base  region. 

Fig.  8  shows  the  temperature- time  slopes  on  the  windward,  leeward 
meridians  and  at  two  locations  on  the  base  region  for  M  =  9.89.  The 
slopes  in  the  meridian  planes  stay  relatively  constant  below  an  angle 
of  attack  of  10°;  however,  beyond  that  angle  they  increase  rapidly. 

It  is  interesting  to  note  the  nonuniformity  of  the  base  temperature 
distribution  beyond  10°  angle  of  attack.  As  the  angle  of  attack 
increases  the  wrapped  shock  heats  up  the  near  wake  and  base  region  in 
the  windward  half  of  the  base  while  the  temperature- time  rate  is 
relatively  low  in  the  base  area  far  from  the  shock.  In  this  area  the 
base  temperature  slope  is  nearly  constant  up  to  27°  of  angle  of  attack 
and  beyond  this  the  increase  is  still  moderate. 

Another  significant  point  is  the  large  difference  in  the  tempera¬ 
ture  slope  between  the  base  region  and  the  leeward  meridian.  Recent 
tests  at  NOL  indicated  that  the  rolled  up  vortex  sheets  generated  at 
angles  of  attack  might  be  reattaching  close  to  the  leeward  meridian 
thereby  causing  a  local  peak  in  heat  transfer.  This  could  then  be  the 
reason  for  the  large  temperature  slope  difference.  The  results  were 
similar  in  the  M  =  6.22  speed  range. 


BASE  PRESSURE  MEASUREMENTS 

Base  pressure  measurements  were  made  with  the  sting  mounted  model 
at  several  angles  of  attack  between  0  and  60°  for  both  M  =  6.22  and 
9.89.  A  sample  of  the  base  pressure- time  history  data  is  presented  in 
Fig.  9.  The  settle-out  time  in  that  particular  run  was  about  20 
seconds.  The  base  pressures  at  both  locations  at  <y  =  0  were  nearly  the 
same.  Furthermore,  this  is  true  even  at  angles  of  attack  as  shown  in 
Figs.  10  and  11,  where  the  base  pressure  ratios  (nondimens ionalized  by 
the  free  stream  static  pressure)  are  presented  as  functions  of  the 
angle  of  attack  for  M=  6.22  and  9.89.  This  observation,  at  least  for 
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Oi  =  0,  was  confirmed  by  other  investigators  who  demonstrated  that  the 
base  pressure  distribution  on  blunt  based  axisymmetric  bodies  exhibited 
a  slight  maximum  in  the  geometric  center  with  an  axially  symmetric 
decrease  toward  the  edge  of  the  base. 

In  view  of  the  axial  symmetry  of  the  base  pressure  distribution, 
it  is  not  unreasonable  that  the  obtained  data  at  zero  angle  of  attack 
show  similar  values,  particularly  since  the  sensing  orifices  were 
located  at  180°  apart  and  equi-distant  from  the  center  axis. 

As  shown  in  Figs.  10  and  11,  the  base  pressure  ratio  is  nearly 
constant  below  about  15°  angle  of  attack  and  then  increases.  This  was 
also  confirmed.  Ref.  15,  for  a  5°  cone  angle.  Beyond  about  40° ,  because 
of  the  sting  effects,  no  steady  base  pressure  value  could  be  reached. 

A  fair  amount  of  confidence  in  the  data  obtained  here  may  be  gained 
by  comparing  it  to  the  work  of  other  investigators.  Based  on  a  large 
number  of  experimental  measurements  conducted  in  the  M  =  7.7  to  19.0 
speed  range,  both  for  laminar  and  turbulent  flows,  an  empirical  base 
pressure  correlation  was  developed  recently  for  a  10°  half  angle  sharp 
cone  at  zero  angle  of  attack.  Fig.  12  shows  the  base- to-f rees tream 
static  pressure  ratio  as  a  function  of  both  the  freestream  Mach  and 
Reynolds  numbers.  The  correlation  equation  may  be  expressed  as: 


0.75 


(5) 


where  the  base  diameter  is  the  characteristic  length  in  the  Reynolds 
number  computation.  The  measured  base  pressure  ratios  at  M  =  6.22  and 
zero  angle  of  attack  are  apparently  close  to  the  predicted  value  of 
Equation  (5).  In  our  case  M  A/Re  =  0.53  and  consequently  from  Fig.  L 
P^/P  «-0.21.  This  is  shown  in  Fig.  10. 

B  oo 

Theoretical  calculation  in  Ref.  16  (showing  the  effect  of  cone 
angle  bluntness  ratio  and  Mach  number  on  the  base  pressure  ratio) 
predicts  base  pressure  ratio  values  of  0.17  for  M  =  6.22  and  0.19  for 
M  =  9.89.  Both  are  in  good  agreement  with  the  measured  values  at 
a  =  0  considering  all  the  uncertainties.  These  values  are  shown  in 
Figs.  10  and  11. 

CONE  SHOULDER  AND  NEAR  WAKE  SURVEY 

The  data  from  this  survey  has  not  been  completely  analyzed  yet 
but  preliminary  indications  are  that  the  band  support  does  not  disturb 
the  near  wake  significantly.  However,  a  more  thorough  examination  of 
the  data  and  some  additional  measurements  are  forthcoming. 
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TWO -DIMENSIONAL  WAKE  SURVEY 

The  objectives  of  this  part  of  the  investigation  were  to  determine 
the  effect  of  the  band  support  on  the  flow  field  and  demonstrate  that 
it  is  smaller  than  the  effect  of  the  wire  support  which  was  claimed  to 
have  no  interference  with  the  near  wake  flow.  According  to  Ref.  17,  a 
critical  wire  diameter- to-base  diameter  ratio  range  exists  below  which 
there  is  no  interference  due  to  the  wire  support.  They  found  that  this 
critical  diameter  ratio  is  in  between  0.0031  and  0.0062.  In  the  present 
case  the  critical  wire  size  would  be  0.012  inches  in  diameter.  To 
begin  with,  this  wire  would  support  only  about  0.5  lb.  and  in  the 
present  configuration  the  maximum  load  is  on  the  order  of  15  lbs.  This 
would  necessitate  the  use  of  30  wires  which  is  unfeasible.  A  more 
realistic  wire  size  would  be  the  0.02  inch  diameter  which  would  exceed 
the  prescribed  size  but  which  would  have  the  load  bearing  capacity  of 
about  5  lbs.  Both  of  these  sizes  were  tested  together  with  the  band 
support  and  indeed  it  turned  out  that  the  band  support  has  the  smallest 
disturbing  effect  on  the  flow;  its  shock  strength  is  lower  than  either  of 
the  wires.  The  schlieren  picture  of  the  wake  of  the  band  support  is 
shown  in  Fig.  13  at  M  =  6.22,  while  for  the  same  Mach  number  the  wake 
shock  for  the  0.02  inch  diameter  wire  is  shown  in  Fig.  14.  The 
quantitative  results  of  the  wake  survey  at  M  =  6.22  are  shown  in 
Figs.  15  to  17.  The  data  are  presented  in  terms  of  nondimensionalized 
quantities.  The  vertical  distance  is  nondimensionalized  by  the  thick¬ 
ness  of  the  band  or  the  wire  diameters  respectively,  and  the  pitot 
pressure  is  divided  by  the  freestream  pitot  pressure.  At  freestream 
conditions,  P  /P  is  unity.  When  the  value  is  larger  the  probe 
encounters  a  2  2co  shock  region,  conversely  where  it  is  smaller  the 
wake  region  is  traversed.  Fig.  15  shows  the  vertical  wake  traverses  of 
the  band  support  at  six  horizontal  locations  behind  the  trailing  edge 
of  the  band.  Fig.  16  shows  six  traverses  for  the  0.02  wire  and  Fig.  17 
gives  three  wake  traverses  for  the  0.01  wire.  In  all  three  cases,  the 
wake  widens  downstream.  The  shock  locations  closely  correspond  to  the 
locations  in  the  schlieren  photographs. 

FLOW  VISUALIZATION 

Oil  Flow 

The  surface  streamlines  appear  to  indicate  a  strong  cross  flow  at 
angle  of  attack  which  was  also  noted  by  several  authors  (see  Refs.  1 
and  15  for  example).  They  seem  to  form  an  inflection  line  at  the 
azimuthal  angle  corresponding  to  the  minimum  surface  pressure  and  then 
became  tangent  to  the  separation  line  a  few  degrees  behind  the  minimum 
pressure  line.  Fig.  18  illustrates  the  oil  flow  patterns  in  this 
experiment.  Protractor  measurements  indicated  that  beyond  20  angle  of 
attack  at  both  tested  Mach  numbers,  flow  separation  occurs  at 
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approximately  142°  as  shown  in  Fig.  19.  To  within  the  accuracy  of  the 
present  data,  flow  separation  appears  to  occur  along  a  conical 
generator.  The  appearance  of  a  separation  line  very  close  to  a  local 
minimum  in  pressure  is  consistent  with  the  data  presented  in  Ref.  18 
and  Ref.  1.  The  oil  flow  pattern  exhibited  the  same  behavior  for  both 
the  sting  and  band  supported  models. 


Schlieren  Photos 

Schlieren  photos  indicated  that  the  incipient  shock  on  the  wind¬ 
ward  side  of  the  body  becomes  highly  curved  downstream  of  the  model 
at  high  angles  of  attack.  The  shock  strength  is  very  small  on  the 
leeward  side  shock  and  approximates  a  shock  wave.  A  shock  is  visible 
on  the  underside  of  the  sting  at  high  angles  of  attack  as  shown  in 
Fig.  20.  This  shock  interacts  with  the  incipient  shock  on  the  windward 
side.  The  appearance  of  the  shock  wave  on  the  sting  indicates  that  the 
near  wake  flow  is  supersonic. 

High  Speed  Motion  Pictures 

High  speed  movies  were  taken  of  the  motion  of  the  band  supported 
model.  Analysis  of  these  motion  pictures  indicated  that  the  natural 
frequency  of  yawing  oscillation  of  the  band  supported  model  was  about 
5  cps,  and  no  asymmetric  forces  acted  on  the  band  supported  model  at 
angle  of  attack  and  consequently  it  was  stable  and  at  zero  degree  yaw 
angle  in  the  yaw  plane. 


CONCLUSIONS 

A  preliminary  experimental  investigation,  concerning  the  aerothermo- 
dynamic  characteristics  of  hypersonic  flows  around  a  highly  inclined 
9°  half  angle  sharp  cone  at  Mach  numbers  of  6.22  and  9.89,  has  been 
conducted.  In  addition,  some  work  was  performed  in  developing  suitable 
instrumentation  for  a  detailed  investigation  of  separated  flow  regimes 
in  all  speed  ranges  around  highly  inclined  axisymmetric  bodies.  From 
this  investigation  the  following  conclusions  are  made: 

1.  The  base  pressure  measurements  showed  an  increase  with 
angle  of  attack.  The  zero  angle  of  attack  base  pressure  values  agree 
well  with  those  reported  in  the  literature. 

2.  The  base  temperature  distribution  and  slopes  are  effected 
by  the  angle  of  attack,  but  they  are  lower  than  the  surface  temperatures 
around  the  model. 

3.  The  two-dimensional  wake  survey  showed  that  the  band 
support  has  a  smaller  effect  on  the  flow  field  than  even  thin  wires. 
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4.  The  location  of  the  measured  separation  line  is  in  good 
agreement  with  the  reported  values  of  other  investigators. 

5.  No  significant  forces  act  in  the  yaw  plane  of  angle  of 

attack. 


FUTURE  PLANS 


A  detailed  investigation  of  the  leeward  side  and  near  wake  regions 
of  slender  bodies  at  high  angles  of  attack  in  the  super/hypersonic 
speed  range  is  planned  in  the  near  future.  The  support  system  (discuss¬ 
ed  in  the  present  paper)  and  the  advanced  instrumentation  developed 
(see  Ref.  20)  will  be  utilized  in  the  program. 
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APPENDIX 


DERIVATION  OF  ANGLE  OF  ATTACK  CORRECTION 

Using  the  equation  of  a  simply  supported  cantilevered  beam  with 
a  concentrated  load  on  its  free  end,  the  maximum  deflection  angle 
can  be  expressed  in  the  present  conditions  as: 

to  =  5.68  x  10"4P  (1) 

where  P  =  concentrated  load  (lbs). 

Since  the  model  is  a  9°  half  angle  cone,  the  concentrated  load 
acting  at  the  end  of  the  sting  (which  coincides  with  the  center  of 
the  model)  may  be  expressed  as  the  normal  component  of  the  aerodynamic 
force  on  the  cone  or 

P  =  CN  qS  (2) 

where 

q  =  dynamic  pressure  (psi) 

S  =2.84  in^ 

=  the  normal  force  coefficient  and  according  to 
Ref.  7  it  is  a  linear  function  of  the  angle  of 
attack  (q) o  For  a  sharp  9°  cone  this  function 
may  be  expressed  up  to  a  =  60°  as: 

CN  =  0.04a  (3) 

where  a  is  in  degrees. 

Combining  Equations  (5)  and  (7)  and  rewriting  it  in  a  more 
convenient  form  gives : 

Aa  =  0.00372 

This  equation  was  used  in 
for  sting  deflection. 


P  Pto 
to 


(4) 


the  data  reduction  program  to  correct 


18 


8th  Navy  Symposium  on  Aeroballistics 


SEPARATED  REGION: 


Vol.  1 


LOW  PRESSURE  AND  HEAT  TRANSFER 


BASE  REGION: 


LOW  PRESSURE  LOW  HEAT  TRANSFER  LOW  VELOCITY 

(U)  Fig.  1.  Flow  Field  of  a  Slender  Body  at  Angle  of  Attack  in 

Super/Hypersonic  Flow 


(U)  Fig.  2.  Sting  Mounted  Model  in  the  Test  Section 
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(U)  Fig.  3.  Band  Supported  Model  in  the  Test  Section 


(U)  Fig.  4.  Band  Support  Installed  at  Nozzle  Exit 
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(U)  Fig.  5.  Critical  Sting-to-Base  Diameter  Ratio  as 
Function  of  Mach  Number  in  Laminar  Flow 


(U)  Fig.  6.  Critical  Sting  Length- to-Base  Diameter  Ratio  as 
Function  of  Mach  Number  in  Laminar  Flow 
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(U)  Fig.  7.  Temperature-Time  History  at  ca  =  35°  and  M  =  9.89 


WINDWARD  MERIDIAN 

leeward  meridian 


ANGLE  OF  ATTACK  (DEGREES) 


(U)  Fig.  8. 


Temperature-Time  Rate  as  Function  of  Angle  of  Attack 
at  M  =  9.89 
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(U)  Fig.  9.  Base  Pressure-Time  History  at  <y  =  28°  and  M  =  9.89 


(U)  Fig.  10. 


Base  Pressure  Ratio  as  Function  of  Angle  of 
Attack  at  M  =  6.22 


23 


8th  Navy  Symposium  on  Aeroballistics 


Vol.  1 


(U)  Fig.  11.  Base  Pressure  Ratio  as  Function  of  Angle  of 

Attack  at  M  =  9.89 


m 

(U)  Fig.  12.  Base  Pressure  Correlation  for  10°  Half  Angle 

Sharp  Cone 
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(U)  Fig.  13.  Wake  of  Band  Support  at  M  =  6.22 


(U)  Fig.  14.  Wake  of  a  0.02"  Dia.  Wire  at  M  =  6.22 
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(U)  Fig.  15.  Wake  Survey  of  the  Band  Support  at  M  =  6.22 


(U)  Fig.  16.  Wake  Survey  Behind  a  0.02"  Dia.  Wire  at  M  =  6.22 
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(U)  Fig.  17.  Wake  Survey  Behind  a  0.01"  Dia.  Wire  at  M  -  6.22 
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(U)  Fig.  19.  Flow  Separation  Angle  as  Function  of  Angle  of  Attack 


(U)  Fig.  20.  Schlieren  Photograph  of  Cone  at  <y  =  49°  and  M  =  6.22 
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Paper  No.  2 

AN  EXPERIMENTAL  INVESTIGATION  OF  THE  HYPERSONIC 
AERODYNAMIC  CHARACTERISTICS  OF  SLENDER  BODIES 
OF  REVOLUTION  AT  HIGH  ANGLES  OF  ATTACK* 

(U) 

(Paper  UNCLASSIFIED) 


by 

Robert  Feldhuhn,  Allen  Winkelmann 
and  Lionel  Pasiuk** 

U.S.  Naval  Ordnance  Laboratory 
White  Oak,  Silver  Spring,  Md.  20910 


ABSTRACT.  (U)  An  experimental  investigation  of  the 
aerothermodynamic  properties  associated  with  axisymmetric 
bodies  at  large  angles  of  incidence  has  been  conducted  in 
two  of  the  wind  tunnels  at  the  Naval  Ordnance  Laboratory 
(NOL) .  During  this  study,  surface  pressure,  heat  transfer 
and  static  force  measurements  were  obtained  with  a  slender 
5°  half-angle  cone.  Static  force  measurements  were  also 
obtained  with  a  2/3  power-law  body  and  two  ducted  cone 
configurations.  In  addition,  flow  visualization  experi¬ 
ments  and  a  limited  series  of  Pitot-tube  surveys  provided 
some  information  concerning  the  separated  flow  field  on 
the  leeward  side  of  a  yawed  cone. 


*  The  work  described  herein  has  been  sponsored  by  the 
Advanced  Research  Projects  Agency  under  ARPA  order  number 
905,  Program  Code  No.  6E30,  as  part  of  Project  Defender. 
**Aerospace  Engineer,  Weapon  Dynamics  Division,  Naval 
Ordnance  Systems  Command 
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R, 


R 


n 


=  base  radius 
=  tip  radius  of  curvature 

=  distance  from  the  tip  measured  along  the  surface 
of  the  cone 


=  Stanton  number 


pUC 


w 


T 

\ 

U 

X 

X 


cp 


=  calculated  stagnation-streamline  Stanton  number 
=0  for  a  yawed  cylinder 

=  free-stream  stagnation  temperature 

=  wall  temperature 

=  velocity 

=  distance  from  tip  of  model  measured  along  the 
axis  of  symmetry 

=  distance  of  the  center  of  pressure  from  tip  of 
model  measured  along  axis  of  symmetry 

a  =  angle  of  attack 

y  =  ratio  of  specific  heats 

0  =  boundary -layer  momentum  thickness 

9c  =  cone  half-angle 

p  =  density  of  gas 

T  =  time 

(f>  =  azimuthal  angle  measured  from  most  windward 

meridian  plane 

Subscripts 

e  =  local  conditions  at  the  outer  edge  of  the 

boundary  layer 

°°  =  free-stream  conditions 
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INTRODUCTION 

(U)  Increased  interest  in  the  development  of  maneu¬ 
verable  hypersonic  vehicles  has  provided  the  impetus  for  the 
present  investigation.  Slender  axisymmetric  conical  shapes 
and  slender  power -law  shapes  are  of  interest  from  this 
standpoint  as  center  of  pressure  variation  and  aerodynamic 
drag,  respectively,  can  be  minimized  with  these  configura¬ 
tions  (refs.  1-4).  However,  if  a  body  is  required  to 
achieve  large  lateral  accelerations,  axisymmetric  configura¬ 
tions  must  sustain  a  large  angle  of  attack  and  maintain  a 
very  high  velocity  in  order  to  compensate  for  their  low 
lift  effectiveness  characteristics.  Both  of  these  alterna¬ 
tives  lead  to  rather  substantial  increases  in  convective 
heat  transfer.  Furthermore  the  possibility  of  sustaining 
large  angles  of  attack  immediately  confronts  the  aero- 
dynamicist  with  the  fact  that  there  will  be  large  regions 
of  separated  flow  on  the  leeward  side  of  the  vehicle.  A 
description  of  the  detailed  characteristics  of  the  leeward 
flow  field  is  necessary  in  order  to  evaluate  the  control 
effectiveness  and  wake  flow-field  characteristics  associated 
with  inclined  bodies  of  revolution. 

(U)  It  is  the  purpose  of  this  investigation  to  provide 
some  fundamental  experimental  measurements  concerning 
aerothermodynamic  properties  of  flow  fields  around  highly 
inclined  bodies  of  revolution  at  hypersonic  speeds.  The 
study  is  intended  to  provide  data  in  the  following  areas 
pertinent  to  the  design  of  an  axisymmetric  maneuverable 
hypersonic  vehicle. 

(1)  heat  transfer 

(2)  performance  and  stability 

(3)  separated  flow-field  phenomena 

(U)  The  present  paper  summarizes  the  results  that  we 
have  obtained  and  described  in  detail  in  references  5  and 
6.  Our  investigation  of  separated  flow  phenomena  is 
continuing  and  we  are  presently  planning  to  conduct  a  series 
of  detailed  flow-field  surveys  in  the  flow  field  on  the 
leeward  side  of  an  inclined  cone. 
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Test  Facilities 

(U)  Two  facilities  at  NOL  were  utilized  during  the 
course  of  the  experimental  investigation.  Force,  pressure, 
heat-transf er  and  flow -visualization  experiments  were 
performed  in  the  Hypersonic  Tunnel  at  a  Mach  number  of  6. 
These  experiments  were  conducted  over  a  range  of  Reynolds 
numbers  that  was  large  enough  to  ensure  that  data  was 
obtained  for  laminar  and  turbulent  flow  conditions, 
respectively.  An  investigation  of  phenomena  associated 
with  the  separated  flow  field  on  the  leeward  side  of  an 
inclined  cone  was  conducted  in  the  Supersonic  Tunnel  No.  2 
at  a  free-stream  Mach  number  of  5  and  a  nominal  free-stream 
Reynolds  number  per  foot  of  5  x  10^. 

(U)  An  outline  of  the  test  program  is  given  in 
Table  I. 

Model  Descriptions  and  Instrumentation 

(1)  Conical  Configurations 

(U)  Individual  force,  pressure  and  heat- transfer 
models  whose  exterior  geometry  was  a  5°  half-angle  cone 
with  a  2.5-inch  base  diameter  were  fabricated  from  stain¬ 
less  steel.  Interchangeable  sharp  and  slightly  blunted 
nose  tips  (Rn  =  0.031  inch)  were  constructed  for  these 
models.  A  sketch  giving  the  dimensions  of  this  configuration 
is  shown  in  Figure  1.  The  distance  of  the  pressure 
orifices  and  thermocouple  junctions  from  the  tip  is  also 
tabulated  in  Figure  1.  The  static  force-measurements  were 
carried  out  with  a  slightly  blunted  cone  with  the  same 
exterior  dimensions  as  the  model  shown  in  Figure  1.  The 
force  model,  designated  henceforth  as  configuration  6FN, 
was  mounted  on  an  internal  water-cooled  multi-component 
balance . 

(2)  2/3  Power-Law  Configuration  (6FX-1) 

(U)  A  model  of  a  slightly  blunted  axisymmetric 
body  whose  outside  radius  varies  essentially  with  the  axial 
coordinate  raised  to  the  2/3  power  was  constructed  and 
mounted  on  an  internal  balance.  A  sketch  of  this  configura¬ 
tion  is  shown  in  Figure  2. 
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(3)  Ducted-Cone  Configurations  (6FX-2-1  and  6FX-2-2) 

(U)  Two  ducted  cone  force  models  were  also  con¬ 
structed  to  investigate  the  possibility  of  substantially 
increasing  the  lateral  maneuverability  of  bodies  of  revolu¬ 
tion.  Sketches  of  these  configurations  are  shown  in 
Figures  3a  and  3b. 

(4)  Data  Recording 

(U)  The  outputs  of  the  transducers  (i.e., 
pressure  transducers,  thermocouples,  and  force  balances) 
were  recorded  on  multi-channel  magnetic  tape  data  recording 
systems.  Individual  channels  were  monitored  on  X-Y  plotters 
during  the  experiments . 

Testing  Techniques  and  Data  Reduction 

(U)  Each  model  employed  in  these  experiments  was 
supported  on  a  sting  whose  axis  was  coincident  with  the 
axis  of  the  model.  During  most  of  the  experiments,  the 
model  sting  assembly  was  attached  to  a  pitch  sector  that 
oriented  the  model  at  the  desired  angle  of  attack.  Pressure 
and  heat  transfer  distributions  in  the  circumferential 
direction  were  obtained  by  employing  a  remote  control 
roll-indexing  device  in  conjunction  with  the  pitch  sector. 
Models  mounted  on  a  dog-leg  sting  roll-indexing  device 
assembly  (Fig.  4)  were  used  for  some  flow  visualization 
experiments.  This  mounting  arrangement  facilitated  photo¬ 
graphic  recording  of  the  phenomena  from  different  orienta¬ 
tions  with  respect  to  the  model  while  the  model  maintained 
a  constant  angle  of  attack. 

(U)  Stagnation  pressure  and  temperature  were  measured 
in  the  plenum  chamber  upstream  of  the  nozzle  throat.  Free- 
stream  conditions  were  calculated  from  the  measured  stagna¬ 
tion  conditions  by  assuming  that  the  flow  within  the  nozzle 
behaves  as  a  calorically  perfect  gas  experiencing  an 
isentropic  expansion  .  The  Mach  number  in  the  test  section 
was  determined  from  the  measured  ratio  of  the  Pitot  pressure 
to  the  stagnation  pressure. 

(1)  Surface  Pressure  Measurements 


(U)  Surface  pressure  distributions  were  obtained 
in  the  Hypersonic  Tunnel  with  the  slightly  blunted  5°  half¬ 
angle  cone  pressure  model  mounted  in  a  roll-indexing  device. 
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The  cone  was  pitched  to  the  desired  angle  of  attack  after 
flow  was  established  in  the  tunnel  and  then  rolled  about 
its  axis  in  order  to  obtain  data  in  the  circumferential 
direction.  The  outputs  from  pressure  transducers  that  were 
connected  with  one  windward  and  one  leeward  pressure  tap, 
respectively,  were  monitored  on  an  X-Y  plotter  in  order  to 
minimize  measurement  errors  associated  with  the  time 
response  of  the  pressure  measurement  system. 

(U)  Pressure  distributions  along  the  most  wind¬ 
ward  and  most  leeward  meridian  generators  were  obtained 
with  a  sharp  5°  half-angle  cone  in  the  Supersonic  Tunnel 
No.  2.  Data  was  obtained  at  different  angles  of  attack  by 
pitching  the  model  and  waiting  until  the  measured  pressure 
reached  its  steady  value. 


(U)  Testing  times  for  these  experiments  were  of 
the  order  of  minutes.  The  pressure  transducers  were 
calibrated  just  before  or  just  after  each  test. 

(2)  Heat  Transfer  Experiments 

(U)  Heat-transfer  measurements  were  made  by  the 
transient  calorimeter  technique.  In  these  experiments  the 
model  was  positioned  in  the  test  section  of  the  Hypersonic 
Tunnel  at  the  desired  orientation  and  then  covered  with  a 
retractable  cooling  shield.  The  model  was  cooled  by 
spraying  liquid  nitrogen  over  its  surface  until  the  wall 
temperature  reached  some  predetermined  level.  Data  at  two 
initial  average  wall  to  total  temperature  ratios  (Tw/Tc  ~ 
0.35  and  0.55)  were  obtained  in  this  manner  during  the  high 
Reynolds  number  experiments  (M  =  6 ,  Re^/FT  =  21  x  10^) , 
while  data  at  one  average  wall  temperature  ratio  (Tw/Tc  ~ 
0.55)  were  obtained  during  the  low^Reynolds  number  experi¬ 
ments  (M  =  5.93  Re  /FT  =  2.4  x  10°) .  Flow  was  established 

oo  oo  , 

in  the  wind  tunnel  with  the  model  within  the  cooling  shield. 
The  shield  was  retracted  after  steady  flow  was  established, 
at  the  desired  supply  pressure,  and  the  temperature 
variation  with  time  was  recorded. _  The  heat-transfer  coef¬ 
ficient  was  calculated  at  a  time  t  =  Ty+T2/2  from  the 
following  relationship: 
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where:  Tq  =  stagnation  temperature 

Cw  =  CWo(l+eT)  =  specific  heat  of  17-4PH 
stainless  steel 

Values  of  Cwo  and  e  for  each  interval  of  time  were  deter¬ 
mined  from  curve  fits  to  available  measurements  of  the 
specific  heat  of  17-4PH  stainless  steel. 

( 3 )  Force  and  Moment  Measurement 

(U)  Static  force  and  moment  measurements  were 
obtained  in  the  Hypersonic  Tunnel  with  models  mounted  on 
water-cooled,  multi-component,  internal  strain-gage 
balances.  The  balances  were  calibrated  before  the  experi¬ 
ments  by  placing  weights  at  different  positions  along  the 
balance.  A  calibration  was  also  made  to  account  for  the 
angular  deflection  of  the  balance  and  the  apparent  weight 
of  the  model. 

(U)  The  deduced  lift  and  drag  data  presented 
herein  include  the  drag  associated  with  the  base  pressure. 

(U)  The  volume,  reference  areas,  and  reference 
lengths  of  the  configurations  used  for  the  data  reduction 
are  tabulated  in  Table  II. 

( 4 )  Pitot-Tube  Surveys 

(U)  Pitot-tube  surveys  were  obtained  on  the 
leeward  side  of  a  yawed  cone  by  traversing  a  Pitot  tube 
perpendicular  to  the  free-stream  direction.  The  axis  of 
the  probe  was  aligned  with  the  flow  direction  in  the 
undisturbed  free  stream. 

( 5 )  Flow  Visualization  Techniques 

(U)  During  the  course  of  these  experiments, 
several  standard  wind-tunnel  flow-visualization  techniques 
were  employed  to  illustrate  phenomena  associated  with 
inclined  bodies  of  revolution. 

(U)  (a)  Schlieren  Photography  -  Schlieren  photo¬ 

graphs  were  taken  throughout  the  course  of  the  experiments 
(knife  edge  horizontal) .  A  series  of  photographs  were 
obtained  with  the  sharp  conical  model  mounted  on  a  dog-leg 
sting.  The  dog-leg  sting  was  affixed  to  a  roll-indexing 
device.  Once  the  axis  of  the  roll-indexing  device  was 
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aligned  with  the  flow  direction  in  the  undisturbed  free 
stream,  the  model  could  be  rolled  about  this  axis  while 
maintaining  a  constant  angle  of  attack  (Fig.  4). 

(U)  (b)  Surface  Oil  Flow  Pattern  -  A  qualitative 

measure  of  the  direction  of  the  shear  stress  at  the  surface 
of  the  cone  was  determined  in  the  Hypersonic  Tunnel  by 
observing  the  oil  pattern  that  was  formed  after  the  surface 
was  exposed  to  the  incident  stream.  A  mixture  of  lampblack 
and  silicone  oil  was  applied  to  the  surface  at  random 
points.  The  shield  used  in  the  heat- transfer  experiments 
was  placed  over  the  model  so  that  the  starting  procedure  in 
the  wind  tunnel  did  not  disturb  the  flow  pattern.  Once 
flow  was  established,  the  shield  was  withdrawn  for  approxi¬ 
mately  25  seconds.  Prior  to  the  shut  down  of  the  tunnel, 
the  shield  was  replaced  in  order  to  preserve  the  pattern. 

(U)  Measurements  of  separation  points  were  made 
at  three  stations  along  the  model  with  protractors. 

(U)  (c)  Surface  Sublimation  Experiments  -  A  qualita¬ 

tive  measurement  of  the  surface  shear  stress  and  heat 
transfer  on  the  leeward  side  of  a  sharp  cone  was  provided  by 
the  surface  sublimation  experiments  conducted  in  the  Supersonic 
Tunnel  No.  2.  A  sharp  teflon  cone  (0C  =  5°,  base  diameter 
=  2.5  inches)  was  painted  black  and  then  coated  with 
azobenzene.  Azobenzene  is  an  orange  colored,  low  melting 
point  (<  200°F)  organic  compound  which  is  suitable  for 
illustrating  regions  of  high  shear  and  heat  transfer  in 
supersonic  tunnels  that  are  operated  at  stagnation  tempera¬ 
tures  above  150°F. 

(U)  The  coated  model  was  mounted  in  the  wind 
tunnel  on  the  offset-sting  assembly  (Fig.  4)  and  exposed  to 
the  incident  stream  for  approximately  five  minutes.  The 
transient  start-up  and  shut-down  periods  did  not  alter  the 
pattern . 

(U)  (d)  Vapor-Screen  Photographs  -  A  qualitative 

indication  of  the  aerothermodynamic  properties  within  the 
flow  field  on  the  leeward  side  was  provided  by  the  vapor- 
screen  photographs.  A  sharp  5°  half-angle  cone  was 
mounted  on  the  dog-leg  sting  roll  device  assembly  in  the 
Supersonic  Tunnel  No.  2  (Fig.  4).  Flow  was  established  at 
the  desired  stagnation  pressure  without  heating  the  air  to 
prevent  air  liquefaction.  A  thin  slit  of  light  was  passed 
across  the  test  section  and  the  resultant  pattern  was 
photographed . 
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RESULTS  AND  DISCUSSION 


Pressure  Experiments 

(U)  Some  typical  surface  pressure  distributions 
recorded  on  the  slender  cones  during  these  experiments  are 
shown  in  Figures  5,  6  and  7.  The  circumferential  pressure 
distribution  shown  in  Figure  5  was  obtained  in  the  Hypersonic 
Tunnel  at  angles  of  attack  up  to  approximately  55°.  The 
data  were  fairly  well  predicted  by  the  Newtonian  relationship, 
Cp  =  2cos26  (where  6  is  the  angle  between  the  local  unit 
normal  and  the  free-stream  velocity  vector) ,  on  the  wind¬ 
ward  side  of  the  cone.  The  data  illustrate  the  fact  that 
the  pressure  distribution  on  the  leeward  side  of  a  5° 
half-angle  cone  is  rather  insensitive  to  angle  of  attack 
once  the  angle  of  attack  exceeds  20°.  The  data  also  illus¬ 
trate  the  existence  of  a  relative  minimum  in  the  pressure 
distribution  at  an  azimuthal  angle  of  approximately  125°. 

(U)  An  examination  of  this  pressure  data  along  with 
the  results  of  flow- visualization  experiments  indicates 
that  at  large  angles  of  attack,  flow  separation  occurs  in 
a  region  of  adverse  pressure  gradient  at  an  azimuthal  angle 
of  approximately  132°. 

(U)  The  data  measured  in  the  Hypersonic  Tunnel  were 
recorded  at  nominal  free-stream  Reynolds  numbers  per  foot 
of  2  x  106  and  20  x  106.  While  data  on  the  windward  side 
showed  very  little  dependence  upon  Reynolds  number,  the 
data  obtained  on  the  leeward  side  of  a  slender  slightly 
blunted  cone  did  indicate  a  consistent  dependence  upon  the 
Reynolds  number.  This  dependence  is  illustrated  in 
Figure  6.  In  order  to  demonstrate  conclusively  the  conical 
nature  of  the  flow  field  on  the  windward  side  and  the 
existence  of  a  pressure  gradient  on  the  leeward  side,  the 
pressure  distribution  along  the  most  windward  and  the  most 
leeward  generators  of  an  inclined  sharp  cone  (0C  =  5°) 
were  measured  in  the  Supersonic  Tunnel  No.  2.  The  results 
of  these  experiments,  presented  in  Figure  7,  illustrate 
these  facts. 


(U)  The  observation  of  a  constant  pressure  along  the 
most  windward  meridian  generator  is  a  manifestation  of  the 
conical  nature  of  the  inviscid  flow  field  on  the  windward 
side.  The  flow  field  on  the  leeward  side  shows  a  dependence 
upon  the  free-stream  Reynolds  number  based  upon  distance 
from  tip.  The  results  of  our  measurements  were  fairly  well 
correlated  with  those  of  other  experiments  by  plotting  the 
static-pressure  ratio  (P/P^)  as  a  function  of  the  viscous 
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interaction  parameter  (M^  /  /Re^  )  (Fig.  8) .  At  the  present 
time  there  are  no  analytical  tredlments  which  adequately 
predict  the  leeward  flow  field  of  a  conical  body  at  large 
angles  of  attack. 

Heat  Transfer  Experiments 

(U)  The  heat-transfer  measurements  were  obtained  on 
a  slightly  blunted  (Rn/Rb  =  0.025)  slender  cone  (6C  =  5°) 
in  the  Hypersonic  Tunnel  at  a  nominal  number  of  6.0.  Some 
typical  data  obtained  in  these  experiments  are  presented 
in  Figures  9  thru  13. 

(1)  Zero  Angle  of  Attack 

(U)  The  data  obtained  at  zero  angle  of  attack 
with  the  slightly  blunted  cone  are  presented  in  Figure  9. 

All  of  the  data  obtained  at  the  lower  unit  Reynolds  number 
indicated  that  the  boundary  layer  was  laminar.  The  data 
obtained  at  the  higher  unit  Reynolds  number  indicate  that 
boundary -layer  transition  occurs  for  these  conditions  at 
free-stream  Reynolds  numbers  between  10  x  10^  and  14  x  10^. 

(U)  Since  the  tip  of  the  cone  was  slightly 
blunted  in  these  experiments,  the  local  flow  properties 
can  no  longer  be  determined  from  solutions  to  the  inviscid 
flow  over  a  sharp  circular  cone.  A  calculation  of  the 
laminar  boundary -layer  growth  on  a  slightly  blunted  cone 
was  accomplished  using  the  integral  method  proposed  by 
Wilson  (ref.  8)  to  account  for  variable  Mach  number  condi¬ 
tions  external  to  the  boundary  layer.  The  variations  of 
local  momentum-thickness  Reynolds  number,  local  Reynolds 
number  based  upon  distance  from  the  tip  and  local  Mach 
number  along  the  cone  as  calculated  by  Wilson's  method 
are  presented  in  Figures  14,  15  and  16.  At  zero  angle  of 
attack,  the  conditions  over  which  transition  occured  in 
our  experiments  are  tabulated  below: 

10  x  106  <  Re  <  14  x  106 
00 ,  s 

4.45  <  Me  <4.7 

5.4  x  106  <  Re  <  9.5  x  106 
e ,  s 

950  <  Re  Q  <  1250 

e ,  0 
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(U)  A  calculation  of  the  turbulent  heat  transfer 
to  a  sharp  cone  as  calculated  by  the  Spalding-Chi  method 
(ref.  9)  is  shown  for  the  purposes  of  comparison.  Our 
measurements  of  turbulent  heating  rates  are  approximately 
10-20  percent  higher  than  those  predicted  by  this  method. 

(2)  Cone  at  Angle  of  Attack 

(U)  Typical  heat-transfer  distributions  along 
the  most  windward  streamline  at  angles  of  attack  of  15°, 

20°,  30°,  45°  and  54°  are  illustrated  in  Figures  10a  and 
10b.  A  comparison  of  these  measurements  with  those  obtained 
at  zero  angle  of  attack  (Fig.  9)  indicates  that  boundary 
layer  transition  occurs  on  the  stagnation  streamline  closer 
to  the  tip  of  a  highly  yawed  cone  than  it  does  on  an  unyawed 
cone.  The  variation  of  transition  Reynolds  number  (i.e., 
based  upon  free-stream  properties)  with  angle  of  attack 
is  shown  in  Figure  17  for  one  wall  temperature  ratio  and 
one  unit  Reynolds  number. 

(U)  Most  available  data  concerning  the  variation 
of  boundary-layer  transition  along  the  windward  streamline 
with  angle  of  attack  have  been  limited  to  small  angles  of 
attack  (references  10,  11  and  12).  The  data  presented  in 
references  10  and  11  indicate  that  transition  moves  rearward 
on  the  windward  meridian  of  a  cone  for  small  angles  of 
attack  while  the  data  presented  in  reference  12  indicate 
that  transition  moves  forward  along  the  windward  meridian  of 
a  sharp  cone  for  small  angles  of  attack.  A  method  which 
predicts  the  location  of  transition  on  simple  shapes  is  not 
available  at  the  present  time.  For  a  yawed  cone,  the 
problem  is  complicated  by  the  presence  of  cross-flow  and 
varying  Mach  number  and  unit  Reynolds  number  conditions  at 
the  edge  of  the  boundary  layer . 

(U)  The  measured  laminar  heat-transf er  distri¬ 
butions  along  the  stagnation  streamline  compare  favorably 
with  those  predicted  from  calculations  of  the  laminar  heat- 
transfer  distribution  on  the  stagnation  streamline  of  a 
yawed  cylinder  whose  axis  is  parallel  to  the  most  windward 
generator  of  the  cone  (ref.  13) 

(U)  Some  representative  heat- transfer  distri¬ 
butions  in  the  circumferential  direction  are  shown  in 
Figures  11,  12  and  13.  For  the  purposes  of  presentation, 
the  data  are  normalized  by  the  calculated  laminar  heat 
transfer  coefficient  at  the  stagnation  streamline  of  a 
yawed  cylinder  (ref.  13) .  The  data  shown  in  Figure  11  are 
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representative  of  the  laminar  heat-transfer  distribution  on 
the  windward  side  of  a  yawed  circular  cone.  The  peak  heat 
transfer  rates  occur  along  the  stagnation  streamline.  A 
second  local  maximum  in  the  heat  transfer  rate  occurs  along 
the  most  leeward  meridian  generator  of  the  cone.  This 
observation  is  probably  related  to  the  fact  that  the  most 
leeward  generator  acts  like  a  rear  stagnation  point  to  the 
separated  flow.  At  the  present  time  there  is  no  adequate 
method  of  predicting  this  increase  of  heat  transfer  on  the 
leeward  generator.  Pitot  surveys,  surface  oil-flow 
patterns,  surface-sublimation  experiments  and  schlieren 
photographs  all  provide  information  which  indicate  that 
there  are  regions  of  attached  flow  near  the  most  leeward 
meridian  generator  of  the  cone.  These  results  will  be 
presented  in  a  later  section  where  separated  flow  phenomena 
will  be  discussed. 

(U)  The  heat-transf er  data  shown  in  Figure  12, 
beyond  x  =  1  inch,  are  representative  of  data  that  were 
obtained  where  the  boundary  layer  on  the  windward  side  was 
turbulent.  Heat-transfer  coefficients  are  maximum  along 
the  stagnation  streamline  and  decrease  in  magnitude  around 
the  periphery  of  the  cone. 

(U)  Data  shown  in  Figure  13  were  obtained  in  a 
region  on  the  cone  where  boundary-layer  transition  occurred. 
These  results  indicate  maxima  in  the  heat  transfer  distri¬ 
bution  approximately  15°  from  the  most  windward  generator. 
One  explanation  for  this  observation  is  that  transition 
can  occur  closer  to  the  tip  on  generators  that  are  displaced 
from  the  most  windward  generator. 

Force  and  Moment  Measurements 

(U)  In  order  to  provide  an  indication  of  the 
performance  and  stability  characteristics  of  axisymmetric 
maneuverable  hypersonic  vehicles,  a  series  of  static -force 
tests  were  conducted  with  a  slightly  blunted  cone 
(configuration  6FN ) ,  a  slightly  blunted  2/3  power-law  body 
(configuration  6FX-1) ,  and  two  ducted  cone  configurations 
(6FX-2-1  and  6FX-2-2)  in  the  Hypersonic  Tunnel.  The 
measured  lift  and  drag  data  along  with  the  lift  to  drag 
ratios  for  these  configurations  are  shown  in  Figures  18 
thru  21.  These  results  indicate  that  the  maximum  lift  to 
drag  ratio  of  the  ducted  conical  configurations  are  as 
much  as  20  to  50  percent  larger  than  the  maximum  lift  to 
drag  ratios  of  the  simple  bodies  of  revolution. 
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(U)  The  measured  distance  between  the  center  of  pressure 
and  the  tip  for  each  of  these  configurations,  normalized  by 
the  length  of  the  body,  is  shown  in  Figure  22.  The  slightly 
blunted  cone  (configuration  6FN)  exhibited  the  smallest 
variation  in  the  location  of  the  center  of  pressure.  Our 
measurements  indicated  that  the  center  of  pressure  of  a 
slightly  blunted  cone  varied  by  less  than  0.3  percent  of  the 
body  length  over  an  angle  of  attack  variation  from  2  degrees 
to  54  degrees.  Furthermore,  measurements  which  we  have 
made  indicate  that  the  center  of  pressure  location  of  a 
slightly  blunted  cone  can  be  predicted  by  the  following 
relationship : 


R  3 


[1  "  (rT)  cos  9c] 


X 

_££ 

L 


3cos20„  Rn  ^  2 

[1  ~  cos  6c] 


Rn  (l-sin0c) 

K 


COS0 


1  - 


R 

n 

R. 


(l-sin0c) 
COS  0  o 


(2) 


The  value  computed  from  this  formula  is  illustrated  in 
Figure  22  for  the  purposes  of  comparison. 


(U)  When  assessing  the  presented  lift  and  drag 
data,  it  is  important  to  note  that  lift  to  drag  ratio  is 
not  the  only  figure  of  merit  for  a  maneuverable  hypersonic 
vehicle.  An  examination  of  the  equations  of  motion  for  an 
unpowered  vehicle  indicate  that  two  additional  dimensional 
figures  of  merit  are  the  ballistic  coefficient,  w/CDA,  and 


the  wing -loading  parameter  CLA/w.  For  a  given  altitude  and 
velocity,  the  longitudinal  and  lateral  accelerations  are 
proportional  to  the  ballistic  coefficient  and  the  wing¬ 
loading  parameter,  respectively.  In  order  to  calculate 
these  parameters,  a  preliminary  design  analysis  must  be 
conducted  to  obtain  some  estimate  of  the  weight  of  the 
vehicle.  Since  this  is  beyond  the  scope  of  the  present 
investigation,  the  only  attempt  to  compare  the  maneuver¬ 
ability  characteristics  will  be  made  by  showing  the  drag 
factor  at  zero  lift  (Volume/CDA) ,  the  lift  factor  (CLA/ 
Volume)  and  the  lift-to-drag  ratios.  The  ballistic 
coefficient  and  the  wing  loading  parameter  can  be  found 
from  these  factors  once  the  packaging  density  (fT*  = 

we is  determined. 

TTA-N  I  UTMA  9 
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(U)  The  drag  factors  at  zero  lift  for  these  configura¬ 
tions  are  shown  in  Figure  23.  Configurations  6FN  and  6FX-1 
have  essentially  the  same  length,  base  radius  and  drag 
coefficient.  The  larger  volume  afforded  by  the  power-law 
body  accounts  for  the  30  percent  increase  of  the  drag  factor 
of  this  body  beyond  that  of  a  cone  (see  Table  II) .  The  drag 
coefficients  at  zero  lift  of  configurations  6FX-2-1  and 
6FX-2-2  are  2  to  2-1/2  times  that  of  the  simple  conical 
body  (Figs.  18,  20  and  21).  These  increases  are  not 
proportionately  reflected  in  the  drag  factor  at  zero  lift 
because  the  volumes  of  configurations  6FX-2-1  and  6FX-2-2 
are  between  1-1/2  to  2  times  that  of  the  simple  conical 
body  (see  Table  II) .  In  calculating  the  drag  factor  of 
configuration  6FX-2-1  one  notes  that  the  increase  in  drag 
coefficient  was  more  than  offset  by  the  increase  in  volume. 

(U)  A  comparison  of  the  lateral  maneuverability 
characteristics  of  these  configurations  is  shown  in  Figure  24. 
The  results  indicate  that  the  ducted  conical  configurations 
can  sustain  at  least  twice  the  lift  factor  and  a  considerably 
larger  lift-to-drag  ratio  than  the  simple  bodies  of  revolution 
can  maintain  at  a  given  angle  of  attack.  The  ducted  conical 
configurations  achieve  a  given  lift  factor  at  a  smaller 
angle  of  attack  than  the  simple  bodies  of  revolution. 
Furthermore,  for  a  given  lift  factor,  the  lift-to-drag  ratio 
of  the  ducted  configurations  is  generally  as  large  or  larger 
than  the  lift-to-drag  ratio  of  the  simple  bodies  of  revolution. 
For  the  conditions  of  our  experiments,  these  results  indicate 
that  the  maneuverability  characteristics  of  bodies  of 
revolution  can  be  augmented  by  the  addition  of  a  duct.  A 
decision,  concerning  whether  or  not  such  configurations  are 
practical  from  the  standpoint  of  a  maneuverable  hypersonic 
vehicle,  requires  an  understanding  of  design  trade-offs 
which  was  beyond  the  scope  of  our  investigation. 

Separated  Flow  Phenomana  on  a  Yawed  Cone 

(U)  A  complete  understanding  of  the  flow  field  on  the 
leeward  side  of  a  yawed  cone  is  not  presently  available. 
Uncertainties  concerning  the  boundaries  between  viscous  and 
inviscid  regions  of  the  flow  and  the  entropy  distribution  on 
the  leeward  side  prevent  one  from  constructing  a  valid  flow- 
field  model.  It  is  our  intention  to  perform  a  detailed  series 
of  flow-field  surveys  which  will  provide  some  useful  infor¬ 
mation  concerning  these  problems.  However,  some  of  the  data 
which  we  have  already  obtained  do  provide  information  rele¬ 
vant  to  the  construction  of  a  valid  flow-field  model.  This 
information  was  obtained  from  surface-pressure  distributions, 
flow- visualization  experiments,  and  Pitot-tube  surveys. 


43 


8th  Navy  Symposium  on  Aeroballistics 


Vol.  1 

(1)  Surface  Pressure 

(U)  It  was  previously  noted  that  the  pressure 
distribution  along  the  most  leeward  generator  of  a  highly 
yawed  cone  is  dependent  upon  the  viscous  interaction 
parameter  (M^/y/Re^  x)  (Fig.  8).  A  satisfactory  explanation 
of  this  observation ' is  not  presently  available.  However, 
this  result  does  cause  one  to  question  previous  observations 
which  indicate  that  the  flow  field  on  a  highly  yawed  conical 
configuration  is  "essentially  conical"  (References  14,  15 
and  16) .  If  the  flow  on  the  leeward  side  of  a  yawed  cone 
is  determined  by  the  interaction  between  a  viscous  flow  and 
an  inviscid  flow,  as  this  correlation  seems  to  suggest,  the 
assumption  that  the  flow  field  is  conical  appears  to  be 
unjustifiable . 

(U)  A  second  indication  of  non-conical  flow 
characteristics  is  illustrated  by  a  schlieren  photograph 
that  was  obtained  with  a  sharp  5°  half-angle  cone  at  angle 
of  attack  of  20°  (Fig.  25) .  The  photograph  indicates 
non-conical  flow  phenomena  as  the  density  gradients  do  not 
exist  along  straight  lines  emanating  from  the  vertex. 

(2)  Flow-Visualization  Experiments 

(U)  In  order  to  demonstrate  some  separated  flow 
phenomena  on  a  highly  yawed  cone  some  standard  flow-visuali¬ 
zation  experiments  were  conducted. 

(U)  Oil-flow  and  surface-sublimation  experiments 
enabled  us  to  locate  flow  separation  lines  on  the  surface 
of  the  cone.  To  within  the  accuracy  of  the  experiments 
(A<J>  =  ±5°)  ,  the  measured  separation  lines  were  found  to  be 
along  conical  generators.  The  primary  separation  line  is 
shown  in  the  photograph  obtained  during  the  surface  subli¬ 
mation  experiment  (Fig.  26).  The  appearance  of  regions  of 
large  shear  and  heat  transfer  on  the  leeward  side,  of  a 
yawed  cone  was  also  illustrated  in  a  surface-sublimation 
experiment  (Fig.  27).  This  observation  is  consistent  with 
the  heat-transfer  measurements  presented  previously  and 
the  experiments  of  Tracy  (ref.  14)  and  Rainbird  (refs.  15 
and  17) .  The  direction  of  the  shear  near  the  most  leeward 
meridian  was  illustrated  in  a  surface  oil  flow  experiment 
(Fig.  28).  Measured  azimuthal  location  of  the  separation 
lines  are  shown  in  Figure  29. 

(U)  An  attempt  was  made  to  illustrate  some 
properties  of  the  separated  flow  field  between  the  surface 
and  the  shock  wave  on  the-  leeward  side  of  the  yawed  cone. 
Schlieren  photographs  and  vapor  screen  photographs  were 
utilized  in  this  regard.  Conventional  schlieren  photo¬ 
graphs  did  illustrate  the  attached  flow  on  the  leeward  side 
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of  the  cone  at  an  angle  of  attack  of  15°  (Fig.  30).  This 
is  indicated  by  the  appearance  of  a  density  gradient  near 
the  base  which  appears  like  a  Prandtl-Meyer  expansion  fan 
that  one  observes  on  flat-based  bodies  in  a  supersonic 
flow . 


(U)  At  larger  angles  of  attack  (Fig.  31) , 
schlieren  photographs  indicate  the  existence  of  two  density 
gradients  on  the  leeward  side  of  a  cone.  Pitot-tube  surveys 
(Fig.  39)  indicate  that  gradients  do  exist  on  the  leeward 
side  of  a  yawed  cone  in  the  plane  of  symmetry.  From  these 
measurements  we  conclude  that  there  are  apparently  three 
distinct  regions  within  the  flow  field  on  the  leeward  side 
of  a  highly  yawed  cone. 

(U)  Schlieren  photographs  obtained  with  the 
cone  mounted  on  a  dog-sting  assembly  offered  views  from 
different  orientations  with  the  model  mounted  at  a  fixed 
angle  of  attack.  A  typical  photograph  obtained  with  this 
arrangement  is  shown  in  Figure  32.  By  analyzing  a  series  of 
photographs  obtained  from  different  orientations,  we  were 
able  to  define  the  shockwave  in  the  cross-flow  plane 
(Fig.  33  and  34).  The  analysis  employed  to  define  a 
3-dimensional  disturbance  from  a  conventional  schlieren 
photograph  is  described  in  reference  18.  While  our  attempts 
to  define  the  leeward  flow  field  solely  from  such  schlieren 
photos  were  not  entirely  successful  ,  large  density  gradients 
were  found  to  exist  in  the  plane  of  symmetry  and  off  the 
plane  of  symmetry.  We  believe  that  the  gradients  that  exist 
off  the  plane  of  symmetry  are  associated  with  the  existence 
of  imbedded  shock  waves  in  the  flow  field  which  we  observed 
with  vapor  screen  techniques. 

(U)  Photographs  of  the  vapor  screen  pattern 
obtained  using  the  set-up  illustrated  in  Figure  4  are 
shown  in  Figures  35  thru  37.  The  enveloping  bow  shock  wave 
is  illustrated  in  all  of  the  photographs.  Imbedded  shock 
waves  are  discernible  in  the  photograph  shown  in  Figure  37. 
Furthermore,  the  absence  of  crystallized  particles  near 
the  body  in  Figures  36  and  37  suggest  the  existence  of  a 
closed  region  on  the  leeward  side  of  a  yawed  cone  which  is 
characterized  by  a  high  temperature. 

(3)  Pitot-Tube  Surveys 

(U)  Additional  information  concerning  the  leeward 
flow  was  obtained  from  Pitot-tube  surveys.  Two  Pitot  surveys 
that  were  obtained  in  the  plane  of  symmetry  on  the  leeward 
side  of  a  yawed  cone  are  shown  in  Figures  38  and  39. 
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The  pressure  measured  at  the  surface  is  shown  for  the 
purposes  of  comparison.  In  the  survey  recorded  with  the 
model  at  an  angle  of  attack  of  15°,  one  notes  that  the 
Pitot  pressure  at  a  point  0.134  inch  above  the  surface 
is  approximately  10  times  the  pressure  measured  by  the 
static  tap  at  the  surface.  This  result  should  be  compared 
with  the  survey  taken  at  30°  angle  of  attack.  For  the 
latter  condition  the  surface  pressure  is  in  close  agreement 
with  the  Pitot  pressure.  The  apparent  difference  between 
the  Pitot  pressure  and  the  surface  pressure  at  15°  angle 
of  attack  is  associated  with  the  presence  of  an  energetic 
attached  flow  near  the  leeward  meridian  generator.  One 
might  speculate  that  the  attached  flow  is  actually  a 
reattached  flow.  The  reattachment  might  be  the  result  of 
a  pair  of  vortices  formed  from  the  flow  which  separates 
along  the  primary  separation  line  (ref.  17). 

(U)  At  the  present  time,  no  attempt  has  been 
made  to  distinguish  between  viscous  and  inviscid  regions 
of  the  flow.  Our  experiments  have  neither  confirmed  nor 
negated  the  existence  of  a  vortical  singularity  within  the 
flow  field.  A  determination  of  the  existence  or  non-existence 
of  the  vortical  singularity  as  well  as  a  mapping  of  the 
velocity  and  thermodynamic  properties  in  the  flow  field  on 
the  leeward  side  of  a  yawed  cone  are  two  specific  objectives 
of  our  forthcoming  flow- field  experiments. 


CONCLUSION 

(U)  An  experimental  investigation  concerning  the 
aerothermodynamic  characteristics  of  flows  around  highly 
inclined  bodies  of  revolution  has  been  conducted  in  the 
wind  tunnels  at  NOL.  From  this  investigation,  the  following 
conclusions  concerning  hypersonic  flow  around  inclined 
axisymmetric  bodies  are -made: 

(1)  (U)  At  large  angles  of  attack,  boundary-layer 
transition  on  the  windward  streamline  occurs  closer  to  the 
tip  than  at  zero  angle  of  attack. 

(2)  (U)  The  measured  static  longitudinal  aerodynamic 
characteristics  indicate  that  significant  increases  in 
maneuverability  can  be  achieved  by  the  addition  of  a  duct 
around  a  slender  cone.  The  measured  lift-to-drag  ratio  of 
the  ducted  conical  bodies  was  as  much  as  50  percent  greater 
than  the  maximum  lift-to-drag  ratio  of  the  simple  bodies  of 
revolution  for  the  specific  Mach  number  and  Reynolds  number 
conditions  of  the  present  tests.  The  increase  in  the  lift- 


46 


8th  Navy  Symposium  on  Aeroballistics 


Vol.  1 

to-drag  ratio  is  realized  solely  because  of  the  increase  in 
lift  effectiveness  of  these  configurations. 

(3)  (U)  The  location  of  the  center  of  pressure  of  a 
slightly  blunted  cone  exhibited  a  variation  of  less  than 

0 . 3  percent  of  the  body  length  over  an  angle  of  attack 
variation  from  2  degrees  to  54  degrees. 

(4)  (U)  Schlieren  photographs  and  pressure  distributions 
along  the  most  leeward  meridian  generator  demonstrated  the 
existence  of  non-conical  regions  on  the  leeward  side  of  a 
yawed  cone.  The  measured  pressure  distribution  has  been 
correlated  with  other  data  by  plotting  the  data  as  a  function 
of  the  viscous  interaction  parameter  (Mot  /  /Re^  ~) 

(5)  (U)  Flow- visualization  techniques,  Pitot-tube 
surveys  and  heat- transfer  measurements  indicate  that  there 
are  regions  of  attached  flow  on  the  leeward  side  of  a 
yawed  cone  (0C  =  5°  and  a  <  20°).  These  regions  are 
characterized  by  local  maximum  values  of  heat  transfer  and 
surface  shear  stress. 

(6)  (U)  Vapor-screen  photographs  indicate  the  existence 
of  imbedded  shock  waves  on  the  leeward  side  of  the  yawed 

cone . 
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(U)  TABLE.  II 


Configuration 

Reference  Length 

Reference  Area 

Volume 

6FN 

13.967  in. 

4.91  sq. in . 

23.38  cu.in. 

6FX-1 

14.027  in. 

4.91  sq . in . 

29.57  cu.in. 

6FX-2-1 

16.800  in. 

4.26  sq . in . 

47.81  cu.in. 

6FX-2-2 

16.800  in. 

4.26  sq . in. 

36.36  cu.in. 
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(U)  Fig.  3a  Ducted  Cone  -  Configuration  6FX-2-1 
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(U)  Fig.  3b  Ducted  Cone  -  Configuration  6FX-2-2 
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(U)  Fig.  4  Vapor  Screen  Set-up  (Top  View) 
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q>  -  AZMUTHAL  ANGLE 

(U)Fig.  5b  Circumferential  Pressure  Distribution 
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(U)  Fig.  6  Static  Pressure  Variation  Along  the  Most  Leeward  Meridian  Plane-  a  =  20° 
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(U)  Fig.  7  Surface  Pressure  Distribution  Along  5-Degree  Cone 
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(U)  Fig.  8  Correlation  of  Pressure  Distribution  Along  Leeward 
Meridian  Generator 
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(U)  Fig.  10a  Heat  Transfer  Distribution  Along  Stagnation  Streamline 
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(U)Fig.  11  Circumferential  Heat  Transfer  Distributions 
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(U)  Fig.  12  Circumferential  Heat  Transfer  Distributions 
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(U )  Fig.  13  Circumferential  Heat  Transfer  Distributions 
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(U)  Fig.  14  Effect  of  Bluntness  on  Momentum  Thickness  Reynolds  Number 
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(U)  Fig.  15  Effect  of  Bluntness  on  Local  Reynolds  Number 
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(U)  Fig.  16  Effect  of  Bluntness  on  Local  Mach  Number 
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(U)  Fig.  18  Lift  and  Drag  Coefficients  and  Lift  to  Drag  Ratio-Configuration  6  FN 
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CENTER  OF  PRESSURE  -  Xc.p. / L 
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LIFT-TO-DRAG  RATIO  -  L/D  LIFT  FACTOR  "  CLA/ VOLUME 


8th  Navy  Symposium  on  Aeroballistics 


Vol.  1 


79 


8th  Navy  Symposium  on  Aeroballistics 


Vol.  1 


80 


(U)  Fig.  25.  Schlieren  Photograph  of  an  Inclined  Slender  Cone  (a  =  20°) 
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(U)  Fig.  26  Azobenzene  Coated  Model  -  Side  View  -  («=  40°) 
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(U)  Fig.  27  Azobenzene  Coated  Model  -  Leeward  Side  ( <*=  15°) 
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(U)  Fig.  29  Separation  Measurements 
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(U)  Fig.  30  Sch Keren  Photograph  of  an  Inclined  Slender  Cone  (<*=  15°) 
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(U)  Fig.  31  Schlieren  Photograph  of  an  Inclined  Slender  Cone  (  « =  30°) 
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(U)  Fig.  32  Color  Schlieren  Photograph  (black  and  white  print) 
of  an  Inclined  Slender  Cone  («=  40°,  co  =  7 0° ) 
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(U)  Fig.  33  Shock  Wave  Shape  in  a  Cross-Flow  Plane  -  a  =  15° 


88 


8th  Navy  Symposium  on  Aeroballistics 


Vol.  1 

MACH  WAVE 
*  PROJECTION 


t 


(U)  Fig.  34  Shock  Wave  Shape  in  a  Cross-Flow  Plane  -  a  =  40° 
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(U)  Fig.  35  Vapor  Screen  Photograph  -«■  0°  (Camera  ED 
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(U ) Fig.  36  Vapor  Screen  Photograph  -  a*  15°  (Camera IT) 
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(U)  Fig.  37  Vapor  Screen  Photograph  4(f  (Camera  I) 


HEIGHT  (INCHES) 


8th  Navy  Symposium  on  Aeroballistics 


(U)  Fig.  38  Pitot  Tube  Survey 
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(U)  Fig.  39  Pitot  Tube  Survey 
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Paper  No.  3 

ANALYTICAL  INVESTIGATION  ON  LAMINAR  FLOW  FIELD 
AND  HEAT  TRANSFER  ON  LEEWARD  SIDE  OF  A 
SHARP  NOSED  HYPERSONIC  CONE  AT 
LARGE  ANGLE  OF  ATTACK 
(U) 

(Paper  UNCLASSIFIED) 
by 

Paul  K.  Chang,  Mario  J.  Casarella, 
and  Russell  A.  Smith 
The  Catholic  University  of  America 
Washington,  D.C.  20017 


SUMMARY 


(U)  At  a  large  angle  of  attack,  the  flow  over  a  sharp  nosed  hyper¬ 
sonic  cone  separates  on  the  leeward  side,  and  forms  downstream  a  large 
viscous  separated  region  between  the  wrapped  shock  and  cone  surface. 

The  complex  physical  aspects  of  the  attached  boundary  layer  region  and 
the  downstream  separated  flow  have  been  reviewed  and  consequently  the 
problems  to  be  solved  have  been  identified.  An  insight  and  under¬ 
standing  of  the  separated  flow  phenomena  have  been  gained  by  the  per¬ 
formed  analysis,  which  is  based  upon  a  proposed  analytical  model. 

(U)  For  convenience,  the  viscous  flow  field  on  the  leeward  side  is 
divided  into  the  following  regions;  separation  point,  constant  pressure 
mixing,  imbedded  shock  reattachment,  far  downstream  flow,  and  re¬ 
circulating  flow,  where  spiral  vortices  may  be  formed.  By  investigat¬ 
ing  each  of  these  regions  of  the  laminar  separation,  a  simplified 
analysis  of  the  overall  method  has  been  developed  to  determine  a  dis¬ 
torted  body  shape.  This  computational  procedure  is  based  upon  the 
matching  of  the  pressure  distribution  in  the  leeward  meridianal  plane. 

(U)  For  this  analysis,  the  configuration  of  the  distorted  body  shape 
is  considered  as  the  boundary  of  inviscid  and  viscous  flow  field.  The 
boundary  conditions  of  viscous  and  inviscid  flow  are  then  to  be  matched 
along  the  distorted  body  shape.  A  generalized  design  procedure  of  the 
distorted  body  shape  involving  the  laminar  flow  separation  for  given 
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flow  conditions,  such  as  free  stream  Mach  number,  angle  of  attack,  cone 
half  angle  and  wall  temperature  has  been  formulated  by  using  the  limit¬ 
ed  available  experimental  data.  Two  sample  numerical  calculations  have 
been  carried  out  to  determine  the  distorted  body  shape  and  the  values 
of  the  properties  of  flow  and  temperature  along  the  outer  boundaries 
of  the  distorted  body  shape. 

(U)  In  addition,  an  attempt  has  been  made  to  evaluate  the  laminar 
convective  heat  transfer  in  the  neighborhood  of  rear  stagnation  based 
upon  the  proposed  laminar  flow  model.  It  appears  that  the  coefficient 
of  laminar  convective  heat  transfer  may  be  computed  by  the  available 
analysis  for  laminar  flow  over  a  slender  ellipse  rather  than  a  circular 
cylinder . 


(U)  a 
b 
h 
k 
M 

P 

u 


X-x 


Y-y 

Y 

x 


y 

T 

a 


SYMBOLS 

Exponential  parameter  for  velocity  profile 

1 1  I!  II  If  It 

Convective  heat  transfer  coefficient 
Conductivity  of  heat 
Mach  number 
Pressure 

Streamwise  velocity  component 
Reverse  velocity  component 
u*  =  ud/u3 

Transformed  coordinate 
Transformed  coordinate 
Y  =  Y  -Yd 

Cartesian  physical  coordinate,  also  distance  along  the 
surface  measured  from  rear  stagnation  point 

Cartesian  physical  coordinate 

Temperature 

Angle  of  attack 
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SUBSCRIPTS: 

(U)  d 

dividing  stream  line 
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outer  edge  of  shear  layer 

o 

initial  condition  at  separation 

r 

reverse  flow 

s 

rear  stagnation 

stag 

stagnation 

w 

wall 

z 

zero  velocity  line 

3 

flow  property  at  the  edge  of  shear  layer 
immediately  downstream  of  imbedded  shock 
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INTRODUCTION 

(U)  The  considered  range  of  investigation  is: 

Free  stream  Mach  number:  5-15 
Altitude:  10,000  -  100,000  ft. 

Angle  of  Attack:  0°  -  60°  and  possibly  higher 

(U)  The  flow  field  around  a  slender  body  of  revolution  at  large 
angle  of  attack  and  at  hypersonic  free  stream  velocity  is  schematically 
represented  in  Figure  1.  The  viscous  boundary  layer  on  the  body  sur¬ 
face  separates  due  to  the  streamwise  adverse  pressure  gradient  and  the 
separated  flow  region  is  extended  into  the  shock  layer  within  the 
outer  wrapped  shock.  Since  extension  of  this  viscous  separated  region 
becomes  larger  as  the  angle  of  attack  increases,  the  aerodynamic  per¬ 
formance  of  slender  hypersonic  vehicle  as  affected  by  the  character¬ 
istics  of  the  flow  and  heat  transfer  at  high  angles  of  attack  are 
strongly  dependent  on  the  features  of  viscous  separation. 

(U)  A  limited  number  of  experimental  investigations  have  been  carried 
out  for  a  cone  at  relatively  large  angles  of  attack,  in  which  detailed 
measurements  were  attempted  in  the  separated  flow  region  on  the  leeward 
meridian.  Tracy  [1]  performed  the  extensive  experiments  on  the  cone 
separated  region  in  laminar  flow  at  free  stream  Mach  Number  7.95,  and 
Rainbird  [2]  has  investigated  the  turbulent  separated  flow  of  a  cone 
up  to  numerical  Mach  Number  of  4.25. 

(U)  The  phenomena  of  separated  flow  downstream  of  a  circular  cone  at 
angle  of  attack  exhibits  special  unique  features  compared  to  other  sep¬ 
arated  flow  phenomena  which  have  been  investigated  extensively  in  the 
past.  As  shown  in  Figure  2,  the  separated  flows  downstream  of  a 
circular  cylinder  and  a  wedge  of  infinite  span  have  much  in  common  with 
that  behind  a  circular  cone,  but  the  essential  nature  of  the  flow 
structures  are  distinctly  different,  because  the  separation  downstream 
of  a  cone  is  a  closed  one  in  the  cross  sectional  plane  perpendicular  to 
the  central  axis,  although  the  separated  flow  downstream  of  a  circular 
cylinder  is  open,  similar  to  wedge  as  well  as  cylinder,  having  a  long 
trail  of  wake. 
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(U)  An  important  special  feature  of  the  leeward  side  of  cone  flow  is 
the  so-called  "vortical  singularity"  introduced  by  Ferri  [3] .  This 
point  may  be  thought  of  as  a  point  of  convergence  for  the  inviscid 
streamline  within  the  wrapped  shock  that  one  considers  solely  in  the 
cross-sectional  plane.  Although  this  topic  has  been  discussed  and 
analyses  performed  proving  its  existence,  it  appears  that  further 
investigations  involving  the  matching  of  inviscid  and  viscous  flow 
regions  are  needed  for  further  understanding. 

(U)  Since  the  flow  structure  of  leeward  separated  flow  downstream  of 
a  circular  cone  is  different  from  that  of  circular  cylinder  and  wedge 
as  discussed  previously,  the  analytical  solution  of  Crocco  and  Lees 
[4]  for  the  latter  cases  matching  the  near  and  far  wake  flow  at  the 
critical  point,  is  not  applicable  for  the  cone  separated  flow.  Thus, 
a  new  approach  is  needed  for  an  analytical  solution  to  the  cone  sep¬ 
arated  flow.  Hence,  a  first  attempt  is  made  to  formulate  a  laminar 
separated  flow  model  which  represents  the  most  important  aspects  of 
the  complex  flow  structure  on  the  leeward  side  of  the  cone.  Then,  a 
simplified  analysis  is  pursued  which  would  both  evaluate  the  flow 
model,  involving  heat  transfer,  by  confirming  existing  experimental 
results,  and,  in  addition,  give  insight  as  to  what  improvements  can 
be  made  in  the  postulated  model. 

(U)  Based  upon  this  model,  characteristics  of  the  attached  and  sep¬ 
arated  shear  layers  of  the  cone  are  computed,  and,  as  shown  in  Figure 
3,  the  distorted  body  shape  is  evaluated  based  upon  the  displacement 
thickness  of  the  shear  layer.  In  order  to  confirm  the  distorted  body 
shape,  a  heat  transfer  analysis  is  made  and  compared  with  the  experi¬ 
mental  data.  If  the  values  of  properties  flow  as  well  as  heat  transfer 
of  inviscid  and  viscous  flow  are  matched  along  the  outer  edge  of  the 
distorted  body,  then  the  analytical  solution  of  the  complete  shock 
layers  are  obtainable. 


the  separated  flow  field  is  divided  into  six 


(U)  For  convenience, 
regions,  which  are: 

Region  1. 
Region  2. 
Region  3. 
Region  4. 
Region  5. 
Region  6. 


Separation  point  region. 

Constant  pressure  free  shear  mixing. 
Imbedded  shock  region. 

Reattachment  point  region. 

Swirling  vortex  region. 

Recirculating  flow  region. 


(U)  An  analysis  of  each  region  is  made  from  which  a  design  procedure 
is  to  be  proposed  for  the  distorted  body  shape.  Some  of  the  analyses 
for  the  respective  regions  made  certain  assumptions  for  which  improve¬ 
ments  may  be  incorporated  whenever  new  experimental  data  are  available. 
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(U)  The  details  for  the  flow  model,  procedure  of  calculation,  results 
of  computation,  are  referred  to  in  references  [5]  and  [6]. 

2.  SZMELXEXEH  ANALYST  £  OF  LAMINAR  SEPARATED  -FL.01L  TN  THE  LEFMAEH 

OF  A  HYPERSONIC  GONE.. 

(U)The  simplified  analysis  of  laminar  separated  flow  in  the  leeward 
side  of  a  hypersonic  cone  is  presented  in  this  section.  Then  in 
Section  3,  simplified  analysis  of  convective  heat  transfer  is  pre¬ 
sented  . 

2a.  SIMPLIFIED  ANALYSIS  OF  VISCOUS  FLOW  BEHAVIOR  IN  THE  LEEWARD 

SIDE  OF  HYPERSONIC  CONE. 

(U)  The  analysis  is  based  on  the  proposed  model  constructed  on 

the  static  pressure  distribution  along  the  leeward  meridian. 

The  assumptions  made  in  the  proposed  computational  scheme  are: 

(1)  Leeward  separated  region  as  well  as  outer  wrapped  shock. 

(2)  The  effect  of  axial  flow  field  to  the  flow  phenomena 
of  cross-flow  components  of  the  separated  shear  layer 
is  negligible. 

(3)  The  qualitative  features  of  the  flow  behavior  are  not 
significantly  affected  by  the  heat  transfer. 

(4)  The  edge  of  the  shear  layer  and  positive  outer  wrapped 
shock  on  the  leeward  meridian  are  known. 

(U)  Assumption  (2)  implies  that  the  proposed  model  is  essentially 
two-dimensional  with  cross-flow  simplifying  the  calculation.  The 
applicability  of  such  computations  are  to  be  judged  finally  in  com¬ 
parison  to  the  pertinent  results. 

(U)  The  analyses  of  six  different  regions,  separation  point,  constant 
pressure  free  shear  mixing,  imbedded  shock,  reattachment,  swirling 
vortex  and  recirculating  flow  are  presented  separately,  using  the 
sketch  of  Figure  4. 

Region  1.  Separation  point  region: 

(U)  The  properties  of  flow  at  the  separation  point  are  referred  to  as 
the  initial  values  for  the  leeward  separated  region.  They  are  cal¬ 
culated  by  the  three-dimensional  attached  flow  analysis,  up  to  sep¬ 
aration,  as  presented  in  Ref.  5  and  6.  A  simplified  two-dimensional 
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analysis  such  as  Cohen  and  Reshotko  [7]  can  also  be  applied;  however, 
a  proper  correction  will  be  necessary  to  compensate  for  the  discrep¬ 
ancies  between  the  two  medthods. 

Region  2.  Constant  pressure-free  mixing  region. 

(U)  The  analysis  of  constant  pressure  free  mixing  determines  the 
streamwise  velocity  distribution,  location  of  dividing  streamline  as 
well  as  zero  velocity  lines. 

(U)  The  boundary  layer  equations  are  reduced  to  the  incompressible 
form  by  using  Howarth  transformation 

Y  ,  P 

X  =  x,  Y  =  /  (  -  )  dy 

o  o  J 

Ke 

where  (x,  y)  coordinates  originate  at  the  separation  point  and  x  is  the 
distance  measured  in  the  downstream  direction  and  y  is  the  normal  to  x. 

(U)  The  sketch  of  velocity  profile  for  constant  pressure  mixing  region 
is  shown  in  Figure  5.  Consider  a  velocity  profile  of  the  form 

u 

-  =  f(n)  =  (1  -  na)b  a>  l  (i) 


where 


n  = 


Y  +  52 


Y  =  Y  -  Yh  (5)  ,  C  = 


vex 


A  =  6  t  +  6  r 


The  boundary  conditions  are 


u  =  u 

e 

at 

Y  = 

6i 

or 

n  =  l 

u  =  o 

at 

Y  = 

6 

z 

or 

Z3 

II 

-Z5 

N 

M 

ii 

3 

at 

Y  = 

62 

or 

n  =  0 

The  constants  a  and  b  are  selected  in  such  a  way  that  the  velocity 
profile  has  an  initial  displacement  thickness  6  and  momentum  thickness 

0  at  the  separation  point,  in  agreement  with  tfie  initial  values. 

o 


102 


8th  Navy  Symposium  on  Aeroballistics 


by  introducing 


Vol.  1 


=  1  -  A  •  f(n) 


where 


The  velocity  along  the  dividing  streamline  is  then  given  by 


u 

e 


where 


=  l  -  A  f(nd) 


nd  =  62  /  a 

The  location  of  the  zero  velocity  line  is 

1  1 

Yz  +  5  2  1  Z  T 

X  f(n  )  =1  or  nz  =  -  +  [  1  -  ( - )b  ] 

A  * 

In  the  physical  plane 

Y  =  (n  -  n  ,)  A 
z  z  d 

The  location  of  dividing  stream  line  can  be  computed  by  the  expression 


Yd  (5) 

6*~~ 


=  -  0.3746 


■)  yir 


(U)  A  unique  solution  is  obtained  by  matching  the  shear  stress  at  the 
dividing  streamline  and  maintaining  a  mass  balance  of  the  recirculating 
flow  with  the  flow  between  the  dividing  stream  line  and  zero  velocity 
line . 


(U)  For  the  simplest  form  of  the  approximation,  py  is  considered  con¬ 
stant  and  Prandtl  number  as  unity.  Furthermore,  it  is  assumed  that 
boundary  layer  approximations  are  valid  in  the  mixing  layer,  provided 
that  the  dividing  stream  line  is  prescribed  since  it  depends  on  the 
interaction  with  the  outer  flow  and  is,  therefore,  related  to  the  higher 
order  terms. 
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(U)  By  substituting  the  velocity  profile  given  in  equation  (1) ,  the 
integrated  form  of  the  transformed  boundary  layer  equations,  is  ob¬ 
tained  for  the  shear  flow  above  the  dividing  stream  line  as 

d  1  6^  X 

-  (ax/  f(n)  [  l  -  X  f(n)  ]  dy  } - 2 - f  (n  )  (2) 

d£  nd  A  a 

and  for  the  shear  flow  below  the  dividing  stream  line  as 


-77  {  a  x  /nd  [  1  -  f (n)  i  [  l  -  Xf (n)  ]  dy  }  =  — - £1  (n.) 

Q-s  o  A  ^ 

(U)  By  matching  the  shear  stress  along  the  dividing  streamline  at 
each  positive  x,  and  expression  for  the  thickness  of  the  mixing  region 
related  to  the  other  parameter  is  obtained  as: 


(3) 


where 

=  /  f(n)  [  1  -  Af(n)  ]  dn 

^d 

nd 

^2  =  f  [  1  -  f(n)  ]  [  1  -  Af(n)  ]  dn 

0 

In  addition,  from  the  required  condition  of  mass  flow  between  the 
dividing  stream  line  and  zero  velocity  line  to  be  equal  to  the  reverse 
mass  flow  of  the  lower  stream,  one  obtains 


A 


f(n)  dn 


(4) 


(U)  A  unique  numerical  solution  for  the  shear  layer  flow  characteristics 
is  obtained  by  numerically  integrating  equation  (2)  along  with  equations 
(3)  and  (4)  for  the  unknown  A,  A  and  ri^.  These  results  can  be  rep¬ 
resented  in  dimensionless  form  by  introducing  variable 

1  =  x^e 

Rex  6 2  u 

o  e 
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(U)  Figure  6  shows  the  velocity  profile  and  the  value  of  the  dividing 
streamline  velocity  of  ud/ue  vs.  £  obtained  by  numerical  integra¬ 
tion.  The  streamwise  u^  velocity  along  the  dividing  streamline 

u^  is  in  good  qualitative  agreement  with  that  obtained  by  Denison. 

and  Baum  [8]. 

(U)  The  displacement  thickness  above  the  zero  velocity  line  is  ob¬ 
tained  by 

*  1 

6  (O  =  f  [  1  -  Xf (n)  ]  dq 


(U)  A  calculation  shows  that  the  displacement  thickness  in  the  con¬ 
stant  pressure  region  appears  to  remain  approximately  constant,  however, 
the  location  of  the  zero-velocity  line  is  significantly  affected  by 
the  temperature  field. 

Region  3.  Imbedded  shock  region. 

(U)  The  pressure  rise  due  to  imbedded  shock  contributes  to  the  re¬ 
attachment  of  flow  along  the  dividing  streamline.  For  simplicity, 
the  imbedded  shock  is  regarded  as  a  straight  oblique  shock.  Thus, 
the  flow  properties  downstream  of  the  imbedded  shock  are  calculated 
by  using  the  well-known  oblique  shock  compression  based  upon  the 
properties  upstream  of  the  imbedded  shock,  assuming  the  pressure  rise 
is  equal  to  the  pressure  difference  between  the  leeward  surface  pressure 
and  minimum  surface  pressure.  For  sufficiently  large  angles  of  attack, 
the  imbedded  shock  is  scattered  near  the  shear  layer.  It  can  be  ex¬ 
pected  that  as  the  angle  of  attack  becomes  quite  large  the  wrapped 
shock  will  start  to  open  along  the  leeward  plane  and  a  necking  region 
develops .  The  imbedded  shock  may  approach  the  so-called  ’’wake  shock” 
near  the  necking  region  similar  to  that  observed  in  the  hypersonic 
wake  for  a  cylinder.  The  value  of  either  a  or  a/02  and  other  para¬ 
meters  at  which  the  wrapped  shock  starts  to  open  remains  an  important 
unknown  for  which,  to  the  best  of  our  knowledge,  no  data  are  available. 

Region  4.  Reattachment  point  region. 

(U)  The  reattachment  process  is  assumed  to  be  isentropic  compression 
of  streamwise  flow  along  the  dividing  stream  line. 

(U)  At  the  reattachment  point,  which  may  be  considered  the  rear  stag¬ 
nation  point,  the  velocity  is  zero  where  the  dividing  stream  line  and 
the  zero  velocity  line  meet.  The  streamwise  velocity  profile  may  be 
considered  similar  to  that  of  separation  involving  9u/3y  -  0  at  y  =  0 . 
Harper  [9]  shows  that  for  the  incompressible  two-dimensional  flow,  with 
its  direction  perpendicular  to  a  flat  plate,  the  flow  phenomena  is 
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inviscid,  in  the  neighborhood  of  the  flat  plate,  although  upstream  and 
downstream  on  the  wall  from  this  inviscid  zone,  the  viscous  flow  pre¬ 
vails.  This  indicates  that  for  incompressible  flow  near  the  reattach¬ 
ment,  the  compression  process  is  isentropic.  Thus,  for  hypersonic  cone, 
the  reattachment  process  is  assumed  to  be  isentropic  compression  of 
streamwise  flow  along  the  dividing  streamline  and  the  stagnation 
pressure  at  the  reattachment  is  computed  by  the  following  equation 
given  by  Weiss  [10] 

Y-l  9  Ud  2  To  y/Cy-1) 

p  stag  =  p  [  1  +  — - M  (  -  )  (  -  )  ]  (5) 

3  Z  3  u,  T  , 


where  subscript  3  refers  to  flow  properties  at  the  edge  of  shear 
layer  immediately  downstream  of  the  imbedded  shock  and  d  refers  to 
the  properties  along  the  dividing  stream  line.  The  temperature  ratio 
in  this  equation  is  computed  by 


w 


+  (  1  - 


Lw 


)  + 


(  Y-l  ) 


~  M3  u3  (  l-u3)  (6) 


where  u. 


.  ud/u3 


For  the  successful  evaluation  of 


at  the 


reattachment,  necessary  to  determine  the  reliable  value  of  u 

which  depends  upon  the  value  of  5  at  the  location  of  the  imbedded 
shock  as  shown  in  Figure  6. 


Region  5.  Swirling  vortex  region  requiring  experimental 
investigation. 

(U)  The  flow  phenomena  in  this  region  is  quite  complicated  and  least 
understood.  From  Figure  4,  it  is  seen  that  due  to  the  closed  boundary 
of  the  separated  flow  region,  the  mass  rate  of  flow  above  the  dividing 
stream  line  must  flow  away  perpendicular  to  the  meridianal  plane.  The 
velocity  profile  at  reattachment  which  is  similar  to  that  of  separation 
has  an  inflection  point  indicating  the  instability  of  the  flow,  tending 
to  form  a  vortex.  Thus,  it  appears  that  in  the  region  of  downstream 
reattachment,  a  swirling  vortex  is  formed.  The  size  of  the  swirling 
vortex  flow  region  is  approximately  determined  by  the  conservation 
of  mass  flow.  A  critical  question  is  the  location  of  the  vortical 
singularity  or  edge  of  shear  layer  on  the  leeward  meridian.  At  this 
stage  of  the  investigation,  semi-empirical  data  are  used  to  locate  the 
farthest  downstream  edge  of  separated  region  in  the  leeward  meridian. 
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Region  6.  Recirculating  flow  region. 

(U)  In  this  region,  the  magnitude  of  recirculating  or  reverse  flow 
is  small.  Nevertheless,  this  reverse  flow  plays  a  role  as  will  be 
demonstrated  in  Section  3,  by  computing  the  convective  heat  transfer 
in  the  rear  stagnation  region  on  the  body  surface.  The  recirculating 
flow  region  is  considered  a  closed  region  bounded  by  the  dividing 
stream  line,  leeward  meridian,  and  body  surface. 

2b.  DESIGN  PROCEDURE  FOR  THE  DISTORTED  BODY  SHAPE. 


(U)  The  solution  of  the  separated  flow  field  in  the  leeward  plane 
is  obtained  by  matching  the  flow  field  upstream  and  downstream  at  the 
reattachment.  For  the  numerical  computation,  a  pressure  distribution 
along  the  leeward  meridian  is  given  as  sketched  in  Figure  7. 

(U)  This  diagram  serves  for  two  purposes  simultaneously.  One  purpose 
is  for  the  constructing  of  the  distorted  body  shape,  and  the  other  is 
the  establishment  of  boundary  conditions  for  the  outer  boundary  of  the 
distorted  body.  The  latter  is  used  to  match  the  inviscid  and  viscous 
flow  solutions.  In  order  to  obtain  the  final  solution,  a  part  of  the 
procedure  of  calculating  the  pressure  distribution  on  this  diagram 
requires  an  iteration  process.  This  is  because  the  positioning  of  the 
imbedded  shock  is  not  known  apriori.  In  Figure  7,  the  positions  4  and 
5  on  the  leeward  meridianal  plane  are  fixed  using  semi-empirical  data. 
The  point  of  reattachment  is  determined  as  follows: 

(U)  The  slope  of  the  dividing  streamline  at  separation  is  assumed 
equal  to  that  of  the  cross-flow  component  of  the  external  inviscid 
velocity.  This  slope  is  extended  from  the  separation  point  by  a 
straight  line  until  it  intercepts  the  meridianal  plane,  thus  locating 
the  reattachment  point,  point  3.  The  magnitude  of  the  static  pressure 
at  5  is  computed  from  shock  wave  analysis.  The  value  of  static  press¬ 
ure  at  4  may  be  determined  because  the  static  pressure  gradient  be¬ 
tween  4  and  5  is  approximately  constant  for  layer  angles  of  attack 
as  measured  by  Tracy  [1].  Hence,  by  drawing  a  straight  line  between 
5  and  4,  representing  the  constant  static  pressure  gradient,  the 
pressure  at  point  4  is  determined  in  the  pressure  diagram. 

(U)  Now,  consider  the  left  hand  side  of  the  diagram.  Point  1  is 
already  fixed  because  the  static  pressure  at  the  rear  stagnation  point 
is  determined  using  the  surface  pressure  distribution.  Next,  assume 
the  location  of  the  imbedded  shock,  point  2,  which  lies  between  the 
rear  stagnation  point  and  its  reattachment  point.  The  pressure  at 
point  2  on  the  pressure  diagram  is  fixed  by  assuming  the  pressure  is 
constant  between  point  1  and  2. 
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(U)  Then,  draw  a  straight  line  between  point  2  and  4  passing  through 
point  3,  again  assuming  a  constant  static  pressure  gradient.  With  this 
procedure,  the  diagram  is  complete.  However,  since  the  position  of  the 
imbedded  shock  has  been  assumed,  the  diagram  is  thus  far  arbitrary. 

The  correct  diagram  is  obtained  by  a  trial  and  error  method,  in  which 
the  computed  stagnation  pressure  at  point  3  must  match  with  that  ob¬ 
tained  on  the  pressure  diagram. 

2c.  NUMERICAL  EXAMPLES 


(U)  Two  examples  of  the  distorted  body  shape  computation  using  Tracy ?s 
measurement  of  a  =  20°  and  24°  at  M^  =  7.95  for  cone  half  angle  10° 
are  presented  in  Figure  8  and  9. 

(U)  The  results  of  these  two  cases  show  that  agreeable  matching  of 
the  upstream  and  downstream  flows  may  be  made  by  the  proposed  iterative 
computational  scheme.  The  displacement  thickness  is  defined  by 
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)  dy 
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by  taking  the  streamwise  coordinate  along  the  zero  velocity  line 
based  upon  the  computed  values  of  flow  properties  up  to  the  reattach¬ 
ment  point.  By  adding  this  displacement  thickness  above  the  zero 
velocity  line,  the  distorted  body  shape  is  constructed. 

3 .  Analysis  of  laminar  convective  heat  transfer  at  the  rear  stag¬ 

nation  region. 

(U)  An  attempt  is  made  to  analyze  convective  heat  transfer  at  rear 
stagnation  region  of  a  cone,  based  upon  the  proposed  model  presented. 
This  analysis  may  serve  for  two  purposes:  one  purpose  is  to  determine 
whether  the  proposed  model  which  is  based  on  pressure  is  applicable 
for  heat  transfer  computation  and  another  purpose  is  to  establish  a 
computational  procedure  for  heat  transfer. 

(U)  At  first  a  comparative  study  of  heat  transfer  at  the  rear  stag¬ 
nation  region  at  subsonic  and  hypersonic  speeds  is  made. 

(U)  Figures  10  and  11  show  the  flow  field  and  measured  convective 
heat  heat  transfer  around  the  circular  cylinder  at  subsonic  air  speed. 
Figure  11  shows  the  proposed  flow  model  and  measured  convective  heat 
transfer  along  the  circumferential  surface  of  a  hypersonic  cone  at 
angles  of  attack. 
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(U)  Roshko  [11]  ,  after  reviewing  concepts  in  the  separated  flow  at 
finite  Reynolds  number,  emphasized  the  existence  of  the  closure  of  the 
separated  flow  region  which  contains  the  vortices  as  sketched  in  Fig¬ 
ure  10.  The  reverse  flow  near  the  center  line  combines  with  the  for¬ 
ward  flow  in  the  free  shear  layer  to  form  the  vortex  pattern.  Al¬ 
though  for  subsonic  separated  flow,  unsteadiness  involving  vortex 
shedding  and  turbulence  are  readily  observable,  the  closure  is,  in 
general,  not  obvious.  While  for  supersonic  flow,  closure  of  separated 
region  due  to  the  reattachment  is  unmistakable. 

(U)  Roshko  [11]  used  Hama’s  photographs  of  subsonic  separated  flow 
behind  a  circular  cylinder  to  observe  the  decrease  of  downstream  width 
of  separated  region  indicating  the  existence  of  the  closure  of  sep¬ 
arated  flow  further  downstream.  Hama’s  picture,  taken  in  a  shorter 
time  compared  to  the  vortex  shedding  period  shows  the  vortex  pattern, 
while  another  picture  taken  with  a  longer  time  span  compared  to  vortex- 
shedding  period,  indicates  apparently  the  closure  point  in  the  separated 
region  far  downstream.  In  1958,  Russell  [12]  found  that  if  the  stream¬ 
lines  representing  the  mean  values  of  velocities  measured  in  the  sub¬ 
sonic  separated  flow  are  constructed,  then  the  surface,  by  discontinuity, 
does  not  extend  to  the  infinity,  but  closes. 

(U)  From  Figures  10  and  11,  it  is  seen  that  although  similarities  of 
reverse  flow  structure  in  closed  recirculating  flow  region  and  pressure 
distribution  exist,  between  subsonic  and  hypersonic  separated  flows  if 
the  swirling  vortex  in  the  hypersonic  case  is  not  considered,  there 
exist  also  differences  in  size  and  geometrical  configuration  of  closed 
recirculating  flow  region.  Next  consider  the  convective  heat  transfer. 
Figure  10  shows  the  Giedt’s  [13]  measurement  along  the  cylinder  surface 
at  subsonic  speed.  Since  the  flow  conditions  at  Reynolds  number  larger 
than  140,000  are  turbulent,  one  lowest  curve  representing  Re  =  70800 
is  laminar  flow  case,  to  be  compared  with  the  hypersonic  laminar  flow 
heat  transfer.  After  studying  the  heat  transfer  behavior,  causing  high 
rate  of  heat  transfer  at  rear  stagnation  region,  Giedt  [13]  conceived 
that  the  reverse  flow  boundary  layer  is  formed  on  the  leeward  surface 
of  the  cylinder  and  its  build  up  continues  toward  the  point  of  separa¬ 
tion.  However,  it  was  not  possible  to  prove  this  concept  by  his  ex¬ 
periment,  because  the  attempted  measurement  of  the  build  up  of  the 
boundary  layer  was  not  conclusive  due  to  the  small  response  of  the 
probe  in  the  separated  flow  region. 

(U)  The  convective  heat  transfer  measured  by  Tracy  [1]  on  the  cone  of 
10°  half  angle  exposed  to  =  7.95  are  shown  in  Figure  11.  At  the 
rear  stagnation  point,  a  pronounced  increase  of  heat  transfer  is 
noticed  similar  to  the  case  of  subsonic  speed. 
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(U)  From  this  brief  comparative  study,  the  following  hypothesis  is 
made : 

The  phenomena  of  the  flow  and  heat  transfer  in  the  vicinity 
of  rear  stagnation  point  are  essentially  those  of  the  boundary  layer 
formed  by  reverse  flow. 

3a .  Proposed  analysis  for  the  rear  stagnation  region. 

(U)  For  the  analysis,  it  is  proposed  to  apply  the  Falkner-Skan  sol¬ 
ution  of  two  dimensional  incompressible  laminar  boundary  layer  and 
energy  equations  given  by 

f,"+ff"+i(i-f,2)=o 
e"  +  f  e'  =o 
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6  is  a  constant  and  is  given  by 
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where  the  superscript  (T)  refers  to  differentiation  with  respect  to  ri . 
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(U)  The  subscripts  e  and  w  are  referred  to  outer  edge  of  shear  layer 
and  wall  respectively.  The  coordinate  x  is  measured  along  the  body 
surface  from  the  stagnation  point  and  y  is  perpendicular  to  x  where  C 
is  constant. 


(U)  The  convective  heat  transfer  coefficient  h  around  the  leeward 
surface  of  a  cone  may  be  computed  then  by 
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(U)  The  applicability  of  the  Falkner-Skan  solution  to  the  flow 
downstream  of  separation  is  shown  by  Stewartson  [14].  The  numerical 
solution  is  obtained  rapidly  by  these  equations,  in  particular  using 
the  already  tabulated  results  such  as  published  in  the  monograph  of 
Evans  [15]. 

(U)  The  Falkner-Skan  solutions  are  applicable  to  similarity  con¬ 
ditions  involving  constant  values  of  g  9  but,  in  the  vicinity  of 
stagnation  point  on  the  blunt  body  surface,  g  is  not  constant,  and 
its  value  decreases  from  unity  at  the  stagnation  point,  requiring  the 
correction.  However,  in  the  vicinity  of  the  stagnation  point,  higher 
order  correction  terms  formulated  by  the  product  of  derivative  of  ^ 
with  respect  to  integration  of  (ue/u00)x  with  respect  to  x,  where 
x  is  distance  measured  from  stagnation  point  and  integral  for  (ue/u0o)x 
with  respect  to  x  are  small,  thus,  solutions  involving  only  the  first 
order  term  are  sufficiently  accurate  [16].  The  presented  Falkner-Skan 
solutions,  equations  (8),  (9),  (10),  are  the  first  order  solutions. 

3b ,  Numerical  examples  of  convective  heat  transfer  in  the  rear 
stagnation  region. 

(U)  For  the  numerical  evaluation  of  heat  transfer  in  the  vicinity 
of  rear  stagnation  point  on  a  hypersonic  cone,  the  property  values 
computed  from  the  proposed  flow  model  and  the  Falkner-Skan  solutions 
are  used.  This  calculation  serves  also  to  confirm  the  applicability  of 
proposed  flow  model  to  heat  transfer,  if  the  predictions  are  in  agree¬ 
ment  with  experiment  data.  Although  limited,  the  following  numerical 
examples  show  that  the  predictions  of  heat  transfer  based  upon  the 


111 


8th  Navy  Symposium  on  Aeroballistics 


Vol.  1 

proposed  hypothesis  of  flow  phenomena  in  the  vicinity  of  rear  stag¬ 
nation  are  consistent  with  the  experimental  data. 

(U)  Tracy  [1]  measured  laminar  convective  heat  transfer  around  the 
circumferential  surface  of  a  cone  of  half  angle  10°  at  =  7.95  at 
two  different  angles  of  attack  a  =  20°  and  24°.  For  convenience,  these 
two  test  conditions  are  designated 

Case  I  for  a  =  20° 


Case  II  for  a  =  24° 

(U)  The  results  of  numerical  computations  of  convective  heat  transfer 
coefficients  h  at  the  rear  stagnation  point  are  listed  in  the  fol¬ 
lowing  table  with  the  measured  data  of  Tracy  [1], 

Table  I  Convective  heat  transfer  coefficient  at  rear  stagnation 
point . 
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Measured 
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2.08  2.32 


(U)  The  discrepancy  among  the  computed  values  and  experimental  data 
may  be  attributed  to  effects  of  three  dimensionality  of  flow,  approx¬ 
imate  values  obtained  from  the  flow  model,  viscosity  and  heat  transfer 
effects . 
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(U)  The  predicted  behavior  of  hg  in  the  vicinity  of  the  rear  stag-* 
nation  point  on  the  hypersonic  cone  is  shown  in  Figure  12.  The 

variation  of  g  in  the  vicinity  of  rear  stagnation  point  on  cone  which 
corresponds  to  the  measured  data  of  pressure  and  convective  heat  trans¬ 
fer  indicates  that  the  phenomena  of  flow  as  well  as  heat  transfer  on 
hypersonic  cone  are  smaller  to  those  on  slender  ellipse  stagnation 
region. 

(U)  Apparently  the  geometry  of  the  recirculatory  flow  region  near  the 
body  surface  is  an  important  factor  for  the  stagnation  phenomena. 

Since  the  zero  velocity  line  is  affected  by  heat  transfer,  the  varia¬ 
tion  of  geometry  of  recirculatory  flow  region  due  to  heat  transfer  is 
to  be  investigated. 

(U)  As  concluding  remarks  on  the  heat  transfer  effects  on  the  pro¬ 
posed  flow  it  may  be  said  that  the  geometrical  boundary  and  boundary 
conditions  of  the  distorted  body  may  be  used  for  the  approximate  com¬ 
putation  of  heat  transfer. 


CONCLUSIONS 


(U)  After  studying  the  phenomena  of  the  separated  laminar  flow  region, 
a  simplified  two-dimensional  laminar  analysis  and  design  procedure  has 
been  presented  to  compute  a  distorted  body  shape.  This  analysis  accounts 
for  the  effect  of  the  separated  viscous  flow  at  large  angles  of  attack. 
This  distorted  body  shape  is  based  upon  the  static  pressure  distribu¬ 
tion  along  the  leeward  meridional  plane.  The  limitation  of  the  ob¬ 
tained  results  is  due  to  a  lack  of  experimental  data  on  properties  of 
the  leeward  side  flow  phenomena.  A  simplified  analysis  of  convective 
heat  transfer  in  the  rear  stagnation  region  based  upon  the  proposed 
flow  model  yields  a  reasonable  agreement  with  experimental  data.  This 
has  shed  more  light  toward  the  understanding  of  the  complex  separated 
flow. 
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FIGURE  1  SCHEMATIC  DIAGRAM  OF  FLOW  FIELD  ABOUT  SLENDER  BODY  OF  REVOLUTION 
AT  ANGLE  OF  ATTACK 
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FIGURE  2  COMPARISON  OF  THE  SEPARATED  FLOW  FIELD  FOR  VARIOUS  GEOMETRIC 
CONFIGURATIONS 
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FIGURE  3  SCHEMATIC  DIAGRAM  OF  THE  DISTORTED  BODY  SHAPE 
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REGION  1  yV  REGION  2 
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DIVIDING  STREAM  LINE 
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FIGURE  5  VELOCITY  PROFILE  FOR  CONSTANT  PRESSURE  MIXING  REGION 
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FIGURE  6  NUMERICAL  RESULTS  FOR  CONSTANT  PRESSURE  SHEAR  MIXING 
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FIGURE  7  SKETCH  OF  THE  CROSS-SECTIONAL  PLANE  ALONG  THE  LEEWARD  SIDE 
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FIGURE  8  SAMPLE  COMPUTATION  OF  DISTORTED  BODY  SHAPE (BASED  ON  DATA  OF  TRACY  a=20°) 


SCALE  3  . 

WRAPPED  5C  =10 

OUTER  SHOCK  a  =24 


8th  Navy  Symposium  on  Aeroballistics 


Vol.  1 


rvi 


O 

<r 

CM 


<3 


124 


FIGURE  9  SAMPLE  COMPUTATION  OF  DISTORTED  BODY  SHAPE  (BASED  ON  DATA  OF  TRACY 
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EXPERIMENTAL  VALUES  (TRACY) 
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FIGURE  13  CONVECTIVE  HEAT  TRANSFER  ALONG  THE  CONE  SURFACE 
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Paper  No.  4 

PRESSURE  DISTRIBUTION  ON  BODIES  AT  LARGE  ANGLE  OF  ATTACK 

(U) 

(Paper  UNCLASSIFIED) 
by 

Howard  R.  Kelly 
Naval  Weapons  Center 
China  Lake,  Calif.  93555 


ABSTRACT.  A  simple  engineering  method  is  proposed  for  finding 
pressure  distributions  about  certain  bodies  of  revolution  at  moderately 
large  angles  of  incidence.  It  improves  upon  the  accuracy  of  previous 
methods,  and  is  more  amenable  to  machine  calculation.  Its  simplicity 
makes  it  useful  for  subroutines  to  find  approximate  pressure  distributions 
for  heat  transfer  calculations. 
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INTRODUCTION 


The  prediction  of  pressure  distribution  has  become  increasingly 
important  in  modern  missile  technology.  One  of  the  more  urgent  needs 
for  pressure  data  is  to  aid  in  the  calculation  of  heat  transfer. 

Increasing  speed  of  aircraft  and  missiles  magnifies  the  heating  problem, 
and  increasing  complexity  of  missile  components  tends  to  make  them 
more  vulnerable  to  heating  damage.  In  addition,  the  requirement  of  high 
maneuverability,  requiring  high  angles  of  incidence,  adds  to  the 
complexity  of  the  problem. 

The  first  step  toward  solving  the  problem  of  pressure  prediction  is 
to  improve  our  ability  to  predict  the  pressure  on  bodies  of  revolution. 
Then,  hopefully,  the  effects  of  fins  and  other  appendages  can  be  added 
to  those  of  the  basic  body.  The  flow  field  about  bodies  of  revolution 
in  axial  flow  is  already  fairly  well  known.  An  attempt  is  made  here  to 
improve  the  prediction  of  pressure  distribution  for  bodies  at  moderate 
to  large  angles  of  incidence,  with  flow  separation  on  the  leeward  side 
of  the  body. 

EXISTING  THEORETICAL  METHODS 

A  number  of  theoretical  methods  of  computing  pressure  exist,  each 
with  its  own  limitations,  and  some  of  these  will  be  discussed  separately 
in  succeeding  paragraphs.  The  earliest  attempts  at  theoretical  methods 
included  first-order  approximations,  such  as  that  of  Tsien^,  and  linear¬ 
ized  theory,  first  considered  by  Lighthill.^  The  applicability  of  these 
is  limited  to  very  slender  bodies  at  moderate  Mach  numbers,  and  very 
small  angles  of  incidence.  These  predictions  do  not  compare  well  with 
exact  theory  for  practical  cases,  so  they  are  of  little  interest  here. 

A  much  better  result  was  obtained  by  Van  Dyke^  in  his  study  of  second 
order  theory,  and  an  application  of  his  methods  is  discussed  below, 
since  it  is  useful  at  moderate  supersonic  Mach  numbers.  Flow  at  high 
supersonic  Mach  numbers,  usually  called  hypersonic  flow,  requires  a 
different  approach,  such  as  Newtonian  flow  theory.^  The  only  really 
general  theory,  that  applies  to  all  supersonic  Mach  numbers,  is  the 
method  of  characteristics,  as  described  by  Ferri.5  This  has  been  developed 
for  the  case  of  small  angles  of  incidence,  and  is  considered  to  be  the 
most  accurate  method  available.  This  accuracy  requires  very  long  and 
tedious  calculations,  so  the  method  is  primarily  used  as  a  standard  of 
comparison  for  simpler  methods,  and  is  little  used  for  practical  calcula¬ 
tions  . 
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Many  other  methods  have  been  developed  for  bodies  of  special  shapes. 

A  very  accurate  method  for  axial  flow  over  cones  is  due  to  Taylor  and 
Maccoll^,  and  an  extension  to  cones  at  small  angles  of  incidence  was  made 
by  Stone^.  The  results  of  these  theories  were  tabulated  by  Kopal^’^’W. 
These  results  are  quite  useful  for  cases  where  tabulated  values  can  be 
used.  The  search  for  a  practical  method  for  non-conical  bodies  led  to  the 
shock-expansion  method  of  Eggers  and  Savin-^5  later  modified  by 

Syvertson  and  Dennis^  and  by  Fenter^. 

It  is  evident  that  an  exact  theoretical  description  of  the  flow 
about  a  body  at  incidence  is  often  very  difficult  if  not  impossible. 

Such  an  exact  description  is  not  needed,  however,  for  missile  design 
problems.  Therefore,  an  engineering  approximation  is  sought,  using 
available  theory  and  extending  its  domain  with  semi-empirical  methods. 

The  primary  aim  is  simplicity,  so  the  results  will  be  useful  as  an  input 
to  computer  programs  for  prediction  of  heat  transfer  effects.  The  present 
method  will  be  restricted  to  the  simplest  case,  that  of  a  sharp-nosed 
body . 


PREVIOUS  ENGINEERING  APPROXIMATION 


One  of  the  first  good  engineering  methods  for  supersonic  and 
hypersonic  pressure  distributions  is  due  to  Vendemia1-5 .  He  discusses 
some  of  the  above-mentioned  theories  in  detail  and  chooses  the  most 
useful  for  several  classes  of  bodies.  Since  the  tangent  ogive  is  of 
primary  interest  here,  we  shall  discuss  only  his  method  for  such  a  body 
shape . 


Vendemia  calculates  the  pressure  distribution  along  the  meridian 
lines  (or  generating  lines)  of  the  body  of  revolution.  A  starting 
value  at  the  windward  side  of  the  nose  tip  is  found  from  cone  tables, 
and  pressures  calculated  on  various  meridian  lines  at  the  tip  from 
the  Generalized  Newtonian  Theory,  such  as  used  by  Lees-^: 


C 

P 


C 

p 

max 


.  2 

sm 


6 


(1) 


where  6  is  the  angle  between  the  free  stream  and  the  tangent  to  the 
body  surface.  Then,  for  a  given  meridian,  the  Generalized  Newtonian 
method  is  used  to  compute  pressure  coefficients  along  the  meridian  until 
an  appropriate  matching  point  is  reached,  where  it  can  be  matched  to  the 
Shock  Expansion  Method11.  In  this  way,  the  pressure  distribution  is 
computed  along  each  meridian  to  the  nose-body  junction. 

Vendemia  then  assumes  that  the  pressure  coefficient  along  the 
cylindrical  body  tends  to  approach  zero  asymptotically.  He  computes 
values  of  this  decaying  pressure  coefficient  from  the  method  of  Fenter1^, 
where 
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Ax/ i 

C  =  C  K  (2) 

P  P(x/£  =  l)6 

In  this  equation  A  is  the  nose  length,  x  is  axial  distance  along  the 
body,  measured  from  the  nose  tip,  Ax  the  axial  distance  measured  from 
the  nose-body  junction,  and  K  the  hypersonic  similarity  parameter,  which 
is  the  ratio  of  the  Mach  number  to  the  nose  fineness  ratio: 


M 

o 

SL/d 


(3) 


The  engineering  approximation  of  Vendemia  seems  to  work  very  well 
when  compared  with  experimental  data  from  Perkins  and  Jorgensen-^,  for 
an  angle  of  incidence  of  five  degrees  (a  =  5°) .  Some  discrepancies 
appear,  however,  for  the  case  of  Ot  =  10°.  He  attributes  this  to  the 
presence  of  separated  flow.  This  is  quite  true,  but  his  assumption 
that  the  pressure  coefficient  will  tend  to  zero  along  the  body  is  also 
unrealistic.  Another  disadvantage  of  his  method  is  that  extensive 
interpolation  in  tables  is  required,  especially  at  the  high  angle  of 
incidence.  In  addition,  the  Kopal  tables  used  do  not  always  cover  the 
range  of  values  of  Mach  number  and  cone  angle  needed  for  practical  cases. 
A  new  improved  engineering  method  is  needed  to  treat  the  pressure 
distribution  along  the  body  more  realistically  and,  at  the  same  time,  be 
more  amenable  to  simple  machine  calculation. 


A  NEW  ENGINEERING  METHOD 


CROSSFLOW  ON  A  CYLINDER 


Prediction  of  Lift 

In  1954,  the  present  author  developed  a  method  of  predicting  forces 
on  a  body  of  revolution  at  large  angles  of  incidence-^ .  This  was  based 
on  the  assumption  that  the  axial  flow  and  cross  flow  are  essentially 
independent,  and  that  the  cross  flow  develops  along  the  body  in  much  the 
same  way  that  the  flow  develops  on  a  cylinder  moving  side-ways  through  a 
fluid,  when  started  impulsively  from  rest.  It  is  well  known  that  a 
cylinder  moving  sideways  through  a  perfect  (non-viscous)  fluid  will 
experience  no  drag.  The  actual  drag  on  a  cylinder  in  a  viscous  fluid  is 
due  to  flow  separation  and  shedding  of  vortices.  The  various  theories 
of  lift  on  bodies  of  revolution  at  incidence  neglect  the  effect  of 
viscosity  and  usually  predict  a  linear  dependence  of  lift  on  angle  of 
incidence.  The  actual  lift  is  non-linear  with  angle  at  large  angles, 
due  to  flow  separation  in  the  actual  viscous  flow.  The  viscous  contribution 
to  the  normal  force  on  a  body  of  revolution  was  correlated  with  the 
empirical  drag  data  for  a  cylinder  as  measured  by  Schwabe-1-” . 
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The  success  of  the  prediction  of  lift  by  semi-empirical  methods 
suggests  that  a  similar  procedure  may  be  used  to  predict  pressure 
distributions.  While  the  measurements  of  Schwabe  were  made  in  terms  of 
flow  speeds,  they  were  also  plotted  in  terms  of  pressure  distribution, 
and  these  pressures  were  used  in  an  attempt  to  improve  the  prediction  of 
pressure  on  the  lee  side  of  bodies  at  angle  of  incidence.  The  results 
were  not  satisfactory  for  several  reasons. 

First,  the  Schwabe  measurements  were  made  for  incompressible  flow 
at  low  Reynolds  number.  Other  data  for  fully  developed  cross  flow 
indicate  that  the  pressure  varies  with  both  Mach  number  and  Reynolds 
number . 

Second,  the  necessary  computations  are  quite  involved.  Interpolation 
must  be  made  in  functions  of  both  time  and  azimuthal  angle.  If  one 
also  were  to  include  Mach  number  and  Reynolds  number,  the  corrections 
will  not  be  simple  enough  to  be  practical. 

Last,  the  pressure  correction  cannot  be  expected  to  be  as  successful 
as  the  correction  applied  to  the  lift  force.  This  is  primarily  because 
the  lift  is  the  integral  of  the  pressure,  so  that  errors  in  pressure 
distribution  will  sometimes  be  cancelled  out  in  the  integration  process. 

In  addition,  the  pressure  near  the  ninety-degree  position  in  azimuth,  as 
measured  from  the  windward  side,  will  contribute  little  or  nothing  to 
the  lift.  For  these  reasons,  it  is  much  more  difficult  to  obtain  a 
completely  accurate  pressure  correction  than  to  obtain  a  correction  for 
the  lift  force.  This,  along  with  lack  of  simplicity,  ruled  out  this 
method  of  correction. 

The  Asymptotic  Pressure 

A  simpler  approach  to  the  problem  is  suggested  by  Vendemia^,  as 
described  above.  We  can  modify  his  method  of  exponential  decay  so  that 
the  equilibrium  pressure  is  not  the  free-stream  value,  corresponding  to 
zero  pressure  coefficient,  but  is  a  value  determined  by  existing  cross 
flow.  Howarth^O  has  shown  pressure  coefficient  distributions  about 
cylinders  as  measured  by  Ferri.  From  these  can  be  deduced  the 
asymptotic  values  of  pressure  coefficient  on  a  very  long  cylinder. 

The  principal  pertinent  results  shown  by  Howarth  are  reproduced  here 
in  Fig.  1.  The  pressure  coefficient  on  a  cylinder  in  cross  flow  is 
shown  as  a  function  of  azimuthal  angle  for  four  different  Mach  numbers. 

The  pressure  coefficient  is  found  to  be  insensitive  to  Mach  number  and 

Reynolds  number  in  supersonic  crossflow,  so  if  the  crossflow  Mach 

number  M^  =  M^  sin  a  >  1.0,  the  assumption  is  made  that  the  curve  for  M^  =  1.85 

is  a  good  approximation.  If  the  cross  flow  is  subsonic  (M^  <  1.0), 

then  the  pressure  coefficient  may  be  quite  sensitive  to  Mach  number,  but 
sensitive  to  Reynolds  number  only  at  very  low  Mach  number.  For  simplicity, 
it  is  assumed  that  the  dependence  on  Mach  number  is  quadratic: 
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C  =  A  +  BM  +  CM-. 2  (4) 

P  1  1  V  ' 

The  coefficients  A,  B,  C  are  different  for  each  azimuthal  angle  cj)  . 

After  normalizing  the  curves  to  C  =  1  at  d>  =  0 ,  the  value  of  (C  ) 

p  p  supersonic 

and  the  coefficients  A,  B,  and  C  are  as  follows,  for  thirty-degree 

intervals : 


<p 

egrees) 

(c  ) 

p  supersonic 

A 

B 

C 

0 

1.00 

1.00 

0.0 

0.0 

30 

0.75 

0.24 

0.458 

-0.042 

60 

0.27 

-0.19 

-0.617 

1.583 

90 

-0.08 

-0.85 

0.029 

0.229 

120 

-0.14 

-1.00 

0.467 

-0.083 

150 

-0.16 

-1.07 

0.671 

-0.229 

180 

-0.17 

-1.07 

0.671 

-0.229 

These  values,  when  multiplied  by  the  Newtonian  pressure  coefficient  from 
Equation  1,  determine  the  asymptotic  pressure  coefficients  along  the 
body. 


OGIVE-CYLINDER  AT  LOW  SUPERSONIC  MACH  NUMBER 

No  single  theory  will  provide  a  simple  method  of  estimating  the 
pressure  distributions  on  a  body  for  all  Mach  numbers,  all  angles  of 
incidence,  and  all  body  shapes.  We  therefore  break  down  the  problem 
into  different  categories  and  try  to  make  these  categories  as  general 
as  possible.  One  simple  body  shape  is  the  tangent-ogive-cylinder.  We 
shall  consider  this  in  some  detail  at  both  low  and  high  supersonic  Mach 
numbers . 

A  good  example  of  the  ogive-cylinder  at  low  Mach  number  is  found  in 
the  measurements  at  Mach  2,  reported  by  Perkins  and  Jorgensen-^.  As 
was  mentioned  before,  this  example  was  chosen  by  Vendemia  to  test  his 
method,  which  was  found  to  be  good  at  a  =  5°  and  not  as  good  at  a  =  10°. 

Vendemia  had  used  Newtonian  and  Shock  Expansion  methods  for  computing 
pressure  along  the  nose.  It  was  decided  to  use  the  Van  Dyke  method 
instead,  since  this  method  was  designed  for  low  supersonic  Mach  numbers. 
The  computations  for  the  ogival  nose  were  carried  out  according  to  the 
method  used  by  Dunn21?  an  interpretation  of  the  Van  Dyke  method  for 
practical  computation.  For  the  cylindrical  body,  the  pressure  coefficient 
was  allowed  to  approach  asymptotically  the  expected  pressure  coefficients 
for  a  cylinder.  This  expected  value  is  found  from  the  Newtonian  Law 
for  the  windward  meridian  (see  Eq .  1),  and  the  values  at  other  meridians 
related  to  this  by  ratios  determined  from  Howarth^,  using  the  existing 
crossflow  Mach  number. 
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The  above  method  depends  heavily  on  the  Van  Dyke  second  order  theory, 
which  is  valid  only  if  the  hypersonic  similarity  parameter  K  (Eq .  3) 
is  less  than  unity.  If  this  parameter  is  greater  than  unity,  the  flow 
may  be  called  truly  hypersonic,  and  one  is  tempted  to  simply  replace 
the  Van  Dyke  theory  by  Newtonian  flow  theory.  This  works  moderately 
well  in  some  cases ,  but  fails  on  the  leeward  side  of  the  nose  ( (J)  =  180 

degrees),  since  the  angle  between  the  body  surface  and  the  free  stream 
direction  passes  through  zero  and  becomes  negative. 

The  method  used  by  Vendemia  works  well  at  small  angles  of  incidence, 
but  not  too  well  at  larger  angles.  It  also  becomes  quite  involved  at 
the  larger  angles,  not  only  from  the  several  numerical  tables  that 
must  be  used  for  interpolation,  but  also  the  complication  of  choosing 
a  matching  point  between  the  Newtonian  theory  and  the  Prandtl-Meyer 
expansion  method. 

The  method  used  here  is  to  find  the  Newtonian  pressure  coefficients 
at  the  tip  of  the  nose  and  use  Prandtl-Meyer  expansion  all  the  way  to 
the  nose-body  junction.  This  requires  only  a  minimum  number  of  tabular 
values,  and,  surprisingly,  gives  as  good  results  as  Vendemia* s  method 
at  moderate  angles  of  incidence  (up  to  10  degrees) ,  without  the  complica¬ 
tions  . 

The  pressure  coefficients  along  the  body  are  found  just  as  before. 

The  value  on  a  given  meridian  at  the  nose-body  junction  is  allowed  to 
decay  exponentially  to  the  equilibrium  value  on  that  meridian,  as 
determined  in  the  table  derived  from  Howarth. 


CONE-CYLINDERS 

The  cone  has  one  of  the  simplest  geometries  of  nose  shapes  for  a 
rocket  or  guided  missile,  and  as  such  has  been  the  subject  of  more  thorough 
aerodynamic  study  than  any  other  shape.  Extensive  tables  (8,  9,  10)  of 
aerodynamic  parameters  for  cones  have  been  in  use  since  1949,  and 
refinements  to  the  method  of  using  these  tables  have  been  proposed, 
particularly  by  Ferri22. 

There  are  two  main  disadvantages  to  using  the  tables  for  engineering 
approximation  to  pressure  distribution,  such  as  we  are  discussing  here. 

The  principal  objection  is  that  the  use  of  numerous  tables  for  search  and 
interpolation  in  a  computer  program  may  lead  to  storage  problems.  Also, 
the  tables  are  limited  with  respect  to  Mach  number  and  cone  angle,  and 
may  not  include  some  of  the  cases  of  interest.  Extrapolation  from  these 
tables  is  not  very  accurate. 

The  best  high  Mach  number  pressure  data  for  cones  at  incidence  were 
used  by  Ferri22 ,  and  taken  from  Cooper  and  Robinson22.  The  test  bodies 
consisted  of  a  20-degree  cone  on  a  cylindrical  body  of  four-diameter 
length,  and  were  tested  at  a  Mach  number  of  6.86.  The  method  of  Ferri 
predicts  the  pressure  on  the  conical  nose  quite  well,  but  the  Newtonian 
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theory  of  Grimminger  is  even  better,  except  on  the  leeward  side  at  high 
angles,  as  will  be  noted  later.  Since  the  hypersonic  similarity  para¬ 
meter  is  greater  than  2.4,  the  Newtonian  theory  was  selected  as  suitable 
for  the  conical  nose.  It  is  also  much  simpler  to  program  for  computation 
than  the  Ferri  method,  or  the  direct  use  of  Kopal's  tables.  The  compu¬ 
tation  is  especially  simple  for  the  conical  nose,  since  pressures  are 
constant  along  each  meridian,  by  conical  flow  theory. 

To  compute  the  pressure  on  the  cylindrical  body,  one  must  first  find 
the  pressure  jump  for  each  meridian  by  using  a  Prandtl-Meyer  expansion 
at  the  cone-cylinder  junction.  This  starting  value  on  the  cylindrical 
body  is  then  allowed  to  decay  exponentially  along  the  body  to  an 
equilibrium  value  for  that  meridian  in  exactly  the  same  way  as  was  done 
for  the  ogive-cylinder. 

At  low  supersonic  Mach  number  the  procedure  is  similar.  The  principal 
difference  is  that  the  theory  of  Van  Dyke  is  used  instead  of  Newtonian 
theory  for  the  cone. 


COMPARISON  WITH  EXPERIMENT 


A  few  examples  have  been  chosen  to  illustrate  the  results  of  the 
present  method  at  moderate  angles  of  incidence.  In  some  cases,  other 
methods  are  included  for  comparison. 

OGIVE-CYLINDER  AT  LOW  SUPERSONIC  MACH  NUMBER 

The  results  of  this  calculation  for  M  =  2  for  the  windward  meridian 
( cf)  =  0°)  are  shown  in  Fig.  2  for  angles  of  incidence  of  5°,  10°  and  15°. 
In  addition  to  the  Perkins  and  Jorgensen  data  and  curves  from  the  present 
method,  curves  are  shown  for  the  method  of  Vendemia,  and  for  the  Van 
Dyke  method  along  both  head  and  body,  as  interpreted  by  Perkins  and 
Jorgensen.  Figure  3  shows  corresponding  results  for  a  meridian  at 
azimuth  of  <f>  =  120  degrees. 

In  comparison  with  the  method  of  Vendemia,  the  present  method  gives 
comparable  results  for  azimuth  less  than  90  degrees.  For  greater 
than  90  degrees,  Vendemia’s  results  are  not  as  good,  particularly  at 
a  =  10  degrees.  It  may  be  expected  that  his  method  would  show  even  less 
favorable  results  at  ot  =  15  degrees.  For  the  cylindrical  body,  the 
fallacy  of  allowing  the  pressure  coefficient  to  approach  zero  is  evident, 
since  the  experimental  values  tend  toward  non-zero  asymptotes  in  many 
cases . 

The  interpretation  of  Van  Dyke  by  Perkins  and  Jorgensen  tends  to 
disagree  with  the  present  method  at  the  nose  tip,  for  some  unknown 
reason.  It  seems  evident  from  the  data  shown,  as  well  as  from  others 
not  shown  that  the  present  method  predicts  the  pressure  distribution 
on  a  tangent-ogive  cylinder  at  low  supersonic  Mach  numbers  better 
than  any  other  simple  engineering  method. 
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As  mentioned  before,  the  results  of  this  method  were  not  expected  to 
be  precise,  but  are  better  than  expected.  In  order  to  check  this  method, 
a  few  data  points  have  been  excerpted  from  a  test  by  Lord  and  Ulmann^. 
One  of  their  models  was  an  ogive-cylinder  with  nose  fineness  ratio  of 
about  3  and  Mach  number  of  about  4,  which  gives  a  hypersonic  similarity 
parameter  of  4/3.  Comparison  with  theory  is  made  for  azimuth  meridians 
of  0,  45,  90  and  180  degrees. 

Figures  4-7  show  that  the  present  method  gives  very  good  results 
at  5  and  10  degrees  incidence,  except  for  the  case  a  =  10  degrees, 

(j)  =  180  degrees.  This  method  fails  at  higher  angles  of  incidence,  but 

could  be  used  for  approximate  results  at  a  =  15  degrees. 


CONE-CYLINDERS 

A  comparison  with  experiment  is  made  in  Fig.  8-10.  The  Mach  number, 
as  mentioned  above,  is  6.86  and  angles  of  incidence  are  6.7,  14  and  20 
degrees.  Azimuth  meridians  are  <p  =  0,  60  and  180  degrees.  The  agree¬ 
ment  with  experiment  is  excellent  except  for  cj)  =  180  degrees,  where 
Newtonian  theory  is  not  expected  to  be  valid  for  a  =  14  and  20  degrees, 
since  the  angle  of  incidence  is  greater  than  the  cone  semi-angle,  and 
these  meridians  lie  in  the  aerodynamic  shadow  of  the  cone.  It  is 
obvious  from  Fig.  10  the  C^  =  0  is  a  good  approximation  for  this  case. 

It  is  recommended,  on  this  basis,  that  pressure  coefficients  be 

computed  by  the  above  method,  and  then  C  =0  substituted  on  the  cone 

P 

for  any  position  in  the  aerodynamic  shadow.  The  criterion  for  this  is 
that  sin  6  <  0,  where 

sin  6  =  sin  0^  cos  a  +  cos  0^  sin  a  cos  cp  . 

Good  data  for  cones  at  low  supersonic  Mach  number  (hypersonic 
similarity  parameter  less  than  one)  are  scarce.  The  best  are  shown 
by  Ferri22  for  a  Mach  number  of  1.6,  but  for  cones  only,  with  no  cylindrical 
afterbody.  In  keeping  with  the  need  for  simplicity,  we  prefer  to  avoid 
the  method  used  by  Ferri,  or  the  direct  use  of  Kopal's  tables.  The  obvious 
method  is  to  use  Dunn!s  interpretation  of  the  Van  Dyke  theory,  as  was 
done  for  ogive-cylinders  at  low  supersonic  Mach  number. 

A  comparison  is  made  between  the  Ferri  and  Dunn  methods  in  Fig.  11. 

Since  data  are  available  only  on  the  cone.  Fig.  11  shows  the  theoretical 
results  and  experimental  data  for  a  =  6  degrees  at  azimuths  from  0  degrees 
(windward)  to  180  degrees  (leeward) .  This  comparison  shows  that  the 
Ferri  and  Dunn  methods  both  represent  the  data  equally  well. 

In  Fig.  11,  additional  curves  show  the  results  of  using  Kopal’s  tables 
directly.  Data  from  Ref.  9,  for  yawing  cones,  and  Ref.  10,  for  cones 
at  large  yaw,  are  shown.  Both  show  very  poor  agreement  with  the 
experimental  data  and  with  the  Ferri  and  Dunn  methods.  This  disagreement 
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becomes  progressively  worse  at  larger  angles  a.  Since  the  Dunn  method 
is  simpler  to  use  than  either  Ferri  or  the  Kopal  tables,  it  appears  to 
be  the  best  method  for  our  engineering  approximation  for  the  case  of 
hypersonic  similarity  parameter  less  than  unity. 


CONCLUSION 


Computer  programs  have  been  written  for  the  prediction  of  pressure 
distribution  on  ogive-cylinders  and  cone-cylinders  at  incidence.  These 
may  readily  be  modified  for  use  as  subroutines  in  other  programs.  These 
programs  are  applicable  for  all  supersonic  and  hypersonic  Mach  numbers 
and  for  angles  of  incidence  up  to  about  15  degrees .  Two  sources  of 
error  tend  to  inhibit  accuracy  above  this  angle  of  incidence.  First, 
most  theoretical  methods  give  poor  results  on  the  nose  for  larger  angles. 
Second,  the  presence  of  shed  vortices  along  the  body  tends  to  distort 
the  pressure  distribution  at  larger  angles.  An  attempt  to  correct  for 
the  presence  of  these  vortices  has  failed. 

The  approximation  used  here  for  the  cylindrical  body  may  also  prove 
useful  for  cylinders  with  blunt  nose  shapes  provided  the  pressure 
coefficient  is  known  at  the  nose-body  junction. 
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FIG.  2  -  Comparison  of  Theoretical  and  Experimental  Pressure 
Distribution  on  a  Tangent-Ogive  Cylinder  at  Various 
Angles  of  Incidence:  M  =  2,  (p  =  0  degree 
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Angles  of  Incidence:  M  =  2,  <j>  =  120  degrees 
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FIG.  4  -  Comparison  of  Theoretical  and  Experimental  Pressure 

Distribution  on  a  Tangent-Ogive  Cylinder  at  Two  Angles 
of  Incidence:  M  =  4.04,  (f>  =  0  degree 
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FIG  5  Comparison  of  Theoretical  and  Experimental  Pressure 

Distribution  on  a  Tangent-Ogive  Cylinder  at  Two  Angles 
of  Incidence;  M  =  4.04,  <j>  =  45  degrees 


145 


8th  Navy  Symposium  on  Aeroballistics 


Vol.  1 


FIG.  6  -  Comparison  of  Theoretical  and  Experimental  Pressure 

Distribution  on  a  Tangent-Ogive  Cylinder  at  Two  Angles 
of  Incidence;  M  =  4.04,  <j>  =  90  degrees 
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FIG.  7  -  Comparison  of  Theoretical  and  Experimental  Pressure 

Distribution  on  a  Tangent-Ogive  Cylinder  at  Two  Angles 
of  Incidence;  M  =  4.04,  <j)  =  180  degrees 
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FIG.  8  -  Comparison  of  Theoretical  and  Experimental  Pressure 
Distribution  on  a  Cone-Cylinder  at  Various  Angles  of 
Incidence,  M  =  6.86,  (p  =  0  degrees 
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FIG.  9  -  Comparison  of  Theoretical  and  Experimental  Pressure 
Distribution  on  a  Cone-Cylinder  at  Various  Angles  of 
Incidence;  M  =  6.86,  <p  =  60  degrees 
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FIG.  10 


Comparison  of  Theoretical  and  Experimental 
Distribution  on  a  Cone-Cylinder  at  Various 
of  Incidence;  M  =  6.86,  (J>  =  180  degrees 


Pressure 

Angles 


150 


8th  Navy  Symposium  on  Aeroballistics 


FIG.  11  -  Comparison  of  Experimental  Pressure  Distribution  with 
Theoretical  Curves  for  Various  Azimuth  Angles  cf)  at 
Angle  of  Incidence  a  =  6  Degrees 
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Paper  No.  5 

COMPACT  GAS-TRANSPIRATION  COOLING  SYSTEM  ANALYSIS 

(U) 

(Paper  UNCLASSIFIED) 
by 

R.  W.  Allen  and  R.  W.  Newman 
Applied  Physics  Laboratory 
The  Johns  Hopkins  University 
Silver  Spring,  Md.  20910 


ABSTRACT  (U)  Gas- transpiration  cooling  is  capable  of  providing 
a  controllable  and  re-usable  means  of  protecting  critically  heated 
surfaces  of  high-speed  flight  vehicles.  In  future  system  installations 
made  under  tight  space  restrictions,  lateral  coolant  motion  will  be 
superimposed  on  the  coolant  motion  toward  the  permeable  wall  of  the 
plenum  chamber.  We  call  this  space-restricted  system  ’'compact1'  in 
order  to  distinguish  it  from  systems  equipped  with  large  plenum 
chambers.  A  meaningful  thermal  analysis  must  couple  the  variables  of 
the  compact  gas- transpiration  system  to  external  radiation  and  aero¬ 
dynamic  heating  of  the  prescribed  flight.  This  is  accomplished 
through  a  flexible,  multi-noded  computer  program  which  allows  external 
heating  influences  to  act  on  lateral  coolant  pressure  and  temperature 
in  the  compact  plenum  chambers  of  candidate  compact  transpiration 
systems.  Internal  coolant  pressure  and  temperature  distributions  in 
turn  govern  the  blowing  effects  on  aerodynamic  heating  and  permeable 
wall  heat  conduction.  The  necessary  terms  are  included  in  the  governing 
equations  processed  by  the  computer  program.  Special  compact-plenum 
experiments  were  made  to  determine  plenum-coolant  pressure  drop  and 
momentum-change  parameters  associated  with  lateral  flow.  By  Reynolds 
analogy,  the  corresponding  compact  coolant  channel  heat  transfer  from 
the  wall  to  the  gas  was  inferred. 
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INTRODUCTION 

Although  gas -transpiration  cooling  has  been  researched  in  con¬ 
siderable  detail  (Reference  1),  its  use  in  thermal  protection  service 
is  limited.  It  can  be  said  however,  that  gas -transpiration  cooling 
systems,  like  ducted  cooling  systems,  have  promise  where  cooling  rate 
control,  system  re-usability,  or  geometry  preservation  are  significant 
considerations.  When  compared  to  ducted  cooling  systems,  gas- 
transpiration  cooling  systems  are  inherently  capable  of  making  more 
effective  use  of  on-board  coolant  reserve  because  of  the  external 
blocking  effect. 

Up  to  the  present  time  it  appears  that  existing  and  proposed  gas- 
transpiration  systems  have  employed  amply-sized  plenum  chambers  coupled 
to  a  sufficient  number  of  feeder  lines  so  that  spatial  variations  in 
coolant  pressure  and  temperature  in  the  plenum  were  virtually  precluded. 
However  in  future  installations  (e.g.  leading  edges)  where  critically 
heated  structural  elements  will  provide  less  and  less  space  for  plenum 
passages,  the  necessity  of  designing  for  and  operating  with  reduced 
plenum  size  must  be  considered.  At  the  same  time,  tighter  space 
restrictions  will  undoubtedly  limit  the  number  of  coolant  feeder  lines. 
When  plenum  chambers  are  sized  and  operated  under  these  space 
restrictions,  characteristic  pressure  and  temperature  distributions 
will  develop  in  the  plenum  coolant  and  the  thermal  task  will  revolve 
around  the  determination  of  these  distributions  and  their  interaction 
with  the  external  heating  load. 

The  foregoing  considerations  have  led  to  a  design  concept  called 
the  compact  gas-transpiration  cooling  system  (Reference  2)  wherein 
passages  simultaneously  serve  as  plenum  chambers  and  as  ducts  or 
channels  to  convey  coolant  to  more  distance  downstream  points.  For 
analytical  purposes,  a  computer  program  has  been  devised  which  handles 
compact  plenum  heat  and  mass  flows  in  conjunction  with  external 
thermal  radiation,  aerodynamic  heating  and  structural  heat  conduction. 
The  program  is  multi-noded  and  provides  for  a  variety  of  thermal 
interconnections  between  the  flight  environment,  the  flight  vehicle, 
and  coolant  conditions  within  the  compact  plenum.  The  following 
account  describes  the  finite-difference  engineering  heat  and  mass 
transfer  relations  used  in  formulating  the  computer  program.  Also 
described  are  special  compact- plenum  experiments  which  were  performed 
to  determine  the  effect  of  negative  blowing  on  channel  flow  in  a  bench- 
test  model.  These  experiments  were  necessitated  by  the  lack  of 
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published  information  on  the  channel-flow  case.  Future  plans  are 
indicated  at  the  conclusion  of  this  paper. 


NOMENCLATURE 


blowing  parameter 

Biot  number  hS /k 

specific  heat 

the  rma 1  c  ap  ac i t anc  e 

hydraulic  diameter 

skin-friction  coefficient 

configuration  factor 

Fourier  number  CY0/&2 

transpiration  parameter  mc^6/kA 

heat  transfer  coefficient 

enthalpy 

the  rma 1  c  ond  u  c  t i v i ty 

permeance 

length 

ratio  of  coolant  mass  flux  to  free  stream  mass  flux 

transpiration  mass  flow 

channel  mass  flow 

molecular  weight 

a  constant 

pressure 

heat  flow 

recovery  factor 

Reynolds  number 

Stanton  number 

tempera ture 


stream  velocity 
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v 

w 

X 

Ax 


velocity  normal  to  walls 
width  of  transpiring  section 
external  distance  downstream 
axial  length  of  control  volume 


GREEK  LETTERS 

a  a  constant  or  thermal  diffusivity 

B  flow-geometry  factor 

6  slab  thickness 

e  emittance 

0  elapsed  time 

T|  film  cooling  effectiveness 

p L  absolute  viscosity 

$  9if  functions 

p  density 

a  Stefan-Boltzmann  constant 

T  shearing  stress 


SUBSCRIPTS 
a  air 

aw  adiabatic  wall 

b  bulk 

c  average  coolant  gas  at  external  porous  surface 

f  friction 

g  coolant  gas 

H  enthalpy 

in  inside  surface  or  entrance 

i5j  entrance  and  exit  of  compact-plenum  gas  control  volume 

m  matrix 

o  no  blowing  or  no  film  cooling 

out  external  surface  of  porous  wall 

r  radiation 

s  subsurface 

St  Stanton  number 
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st 


static 


w 


wall 


w  ,a 


wall,  apparent 


SUPERSCRIPTS 


o 


no  blowing 


THE  COMPACT  GAS -TRANSPIRATION  SYSTEM 


In  order  to  establish  the  physical  nature  of  a  typical  problem,  we 
consider  a  hypothetical  leading  edge  in  flight.  The  cross  section 
view  of  Figure  1  shows  an  oncoming  external  flow  passing  through  a 
shock  wave  and  forming  an  aerodynamic  boundary  layer  on  a  flight  body. 
The  severely  heated  part  of  the  leading  edge  in  Figure  1  is  fabricated 
from  a  permeable  porous  metal,  with  the  permeable  segment  backed  up 
by  a  compact  plenum.  The  latter  is  shown  as  a  simple  cylindrical 
flow  passage  running  in  the  spanwise  direction  along  the  leading  edge. 
Portions  of  the  stream  of  gaseous  coolant  flowing  axially  in  the 
compact  plenum  are  forced  into  and  through  the  permeable  wall  by  the 
pressure  being  maintained  in  the  plenum.  The  passing  of  coolant  mass 
from  the  main  stream  into  the  wall  is  usually  called  suction  but  the 
term  negative  blowing  generally  leads  to  a  clearer  understanding  of 
the  process  with  respect  to  channel  flow  in  the  plenum.  After  coolant 
penetrates  the  permeable  wall,  it  transpires  through  it  and  emerges 
on  the  external  side.  Here  it  produces  the  well-known  blowing  effect 
which  tends  to  block  the  incoming  aerodynamic  heat  load.  This  action 
will  be  called  positive  blowing.  Blown  coolant  merges  with,  and 
develops  along  the  external  permeable  surface  up  to  its  juncture  with 
the  solid  external  wall.  There,  a  film-cooling  region  begins  and 
thereafter  the  cooling  effect  dies  out  as  the  boundary  layer  is 
followed  further  aft. 

Overall,  a  compact  gas  transpiration  system  is  made  up  of  com¬ 
ponents  shown  in  Figure  2.  Coolant  flow  can  be  traced  from  the 
reservoir  into  supply  piping,  through  the  coolant  flow-control  device, 
and  into  the  compact  plenum  where  it  flows  spanwise  inside  the  leading 
edge.  Of  foremost  concern  is  the  fact  that  spanwise  variations  in  the 
positive  blowing  rate  and  spanwise  variations  in  temperature  will 
occur  as  the  blown  coolant  emerges  from  the  leading  edge.  They  will 
be  functions  of  internal  coolant  pressure  and  temperature  distribution 
along  the  compact-plenum  flow  axis  as  well  as  the  external  spanwise 
pressure  and  heat  load  variation.  The  primary  goal  of  the  present 
analysis  is  to  develop  a  method  for  determining  the  internal  pressure 
and  temperature  distributions  and  their  effect  on  external  thermal 
protection  performance. 
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PROBLEM  DEFINITION  AND  FINITE  DIFFERENCE  SCHEME 


A  proposed  compact  gas -transpiration  system  layout  together  with 
the  local  structural  configuration  of  the  flight  body  constitutes  the 
region  of  analytical  interest.  About  this  region  of  interest,  specified 
unsteady  flow  and  thermal  conditions  are  given  in  terms  of  flight  time. 

The  initial  thermal  condition  of  the  flight  body  is  prescribed.  A 
specification  is  made  of  exactly  how  the  coolant  mass-flow  control 
device  is  to  operate  and  the  gas  flow  resistance  characteristics  of 
the  porous  wall  are  assumed  to  be  known.  Structural  material  and 
coolant  thermal  properties,  specified  as  a  function  of  temperature, 
complete  the  definition  of  the  problem. 

In  analogic  terms  we  can  say  that  the  mathematical  model  consists 
of  an  array  of  network  nodal  plug-in  points  and  electric  conductance 
elements  which  can  be  plugged  in  between  any  pair  of  nodes.  Some,  but 
not  all,  nodes  are  permanently  connected  to  grounded  capacitors.  These 
capacitors  are  set  to  the  thermal  capacitance  values  of  corresponding 
regions  of  flight-body  material  while  conductances  are  chosen  to  agree 
with  the  thermal  conductance  value  for  the  internodal  path  in  question. 
Provisions  are  made  for  temperature  dependency  of  both  thermal  capacitance 
and  thermal  conductance.  Many  other  nodes  are  not  connected  to  grounded 
capacitors.  These  are  used  for  surface,  external  air,  radiation,  and 
internal  coolant  gas  temperatures.  No  capacitance  is  required  for  air 
duct  coolant  because  gases  are  invariably  found  to  have  negligible 
thermal  capacitance  when  compared  to  the  flight-body  material.  External 
aerodynamic  heat  transfer  is  introduced  in  an  engineering  form  by 
calculating  the  local  thermal  conductance  corresponding  to  the  local 
external  convective  heat  transfer  coefficient  of  the  boundary  layer. 

At  those  points  where  external  heat  transfer  takes  place  between  the 
boundary  layer  and  a  permeable  surface,  the  heat  transfer  coefficient 
is  reduced  in  proportion  to  the  blocking  effect  of  positive  blowing. 

Nodes  for  coolant  gas  temperatures  in  flow  passages  receive  a  net 
coolant  enthalpy  influx.  They  are  also  connected  to  thermal  con¬ 
ductances  corresponding  to  local  pipe-flow  heat  transfer  coefficients. 

At  points  on  the  permeable  surface  of  the  compact  plenum,  the  con¬ 
vective  heat  transfer  coefficient  is  increased  in  proportion  to  the 
effect  of  negative  blowing. 

In  addition  to  heat  transfer  calculations,  coolant  mass-flow 
distribution  over  the  positively  blown  external  surface  is  determined 
by  satisfying  the  continuity  equation  and  momentum  equation  of  fluid 
mechanics  throughout  the  coolant  piping  system  and  permeable  medium. 
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A  grid  of  pressure  nodes  is  established  in  the  coolant  piping  system. 
Coolant  pressure  values,  consistent  with  the  governing  equations,  are 
computed  in  an  iterative  process. 


EXTERNAL  FLOW  AND  SURFACE  HEAT  TRANSFER 

External  flow  conditions  are  specified  in  terms  of  the  Mach 
number,  pressure,  and  temperature  of  the  oncoming  flow  (Figure  1)  as 
a  function  of  flight  time.  Local  flow  conditions  between  the  shock 
wave  and  boundary  layer  are  calculated  at  points  opposite  each  surface 
node.  This  provides  the  adiabatic  wall  enthalpy  needed  to  make  the 
energy  balance  on  a  non-capacitive  surface  area.  Figure  3  depicts 
the  energy-balance  procedure  schematically  with  respect  to  surface 
(1)  located  on  the  permeable  blowing  surface.  Letting  subscript  1 
denote  a  condition  evaluated  at  the  surface,  we  have  the  following 
energy  equations  with  respect  to  surface  (1),  Figure  3 


'  ,  .  =  OA  SF  (T4  -  T4) 

rad ,  in  i  n  r  i  7 

(1) 

convection  in  ^aw,i 

V 

(2) 

=  k  —  (T  - 

conduction  out  is  L  l 

V 

(3) 

i  s 

External  boundary- layer  heat  transfer  coefficient  hR  is  the  net  value 
after  the  positive  blowing  effect  has  been  included  in  accordance 
with  Eckert Ts  (Reference  3)  recommendation 
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Likewise,  the  adiabatic  wall  enthalpy  Haw  is  computed  from  a  recovery 
factor  which  is  the  net  value  after  the  positive  blowing  effect  has 
been  included  (Reference  3). 
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Both  functions  §  and  \|t  usually  have  a  linear  dependency  on  the  positive 

blowing  rate  m  .  Later  on  we  shall  see  how  the  positive  coolant  blowing 
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rate  m  is  computed  when  continuity  and  momentum-law  principles  are 
applied  to  coolant  flow  in  the  compact  plenum  and  in  the  porous  matrix. 

Following  the  external  boundary  layer  along  the  porous  surface  to 
the  solid  surface,  we  encounter  a  region  of  film  cooling  where 


H  =H°  -  T)  (H  °  -  H) 

aw  aw  H  aw  c 


(6) 


gives  the  necessary  change  from  the  non-film-cooled  adiabatic  wall 
enthalpy  Haw  to  the  film  cooled  value  Haw.  According  to  Goldstein, 
Shavit,  ana  Chen  (Reference  4)  the  film-cooling  effectiveness  down¬ 
stream  from  a  porous  section  is  most  accurately  predicted  by  the 
Kutateladze  and  Leont'ev  (Reference  5)  relation 
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In  the  film  cooling  method,  no  adjustment  is  made  in  the  local  heat 
transfer  coefficient  computed  from  local  flow  conditions. 


INTERNAL  HEAT  CONDUCTION  AND  TRANSPIRATION 


An  imaginary  element  buried  in  the  transpired  porous  leading  edge 
receives  net  heat  flow  by  net  heat  conduction  and  net  coolant  enthalpy 
flow  into  the  control  region  defining  the  element.  Figure  3  depicts 
a  porous  element  (2)  and  the  schematic  energy  flows  with  respect  to 
its  control  surfaces.  Thermal  energy  is  stored  in  the  element  by 
virtue  of  the  thermal  capacitance  of  the  metal  matrix.  Turning  to 
Figure  3  we  can  show  the  energy  balance  on  permeable  element  (node) 

(2)  in  more  detail.  One-dimensional  heat  flow  sufficiently  portrays 
the  basic  terms  as  follows 
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where  subscripts  1,  2,  and  3  correspond  to  three  adjacent  permeable 
elements  numbered  in  ascending  order  in  the  direction  of  gas  flow  and 
where  gas  enters  element  2  at  Tag  and  leaves  at  T^g .  We  use  an  explicit 
form  with  nodal  temperature  Tp  occurring  at  the  beginning  of  time 
interval  At  and  T1  occurring  after  At  has  elapsed.  It  will  be  noted 
that  we  have  assumed  that  the  coolant  gas  (subscript  g)  does  not 
contribute  to  heat  conduction  or  thermal  capacitance  but  does  give 
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rise  to  the  net  enthalpy  inflow  m  cpg  (T  -Tbg)#  Only  in  exceptional 
cases  will  porous  passages  be  large  enough  to  require  that  we  differ¬ 
entiate  between  the  local  temperature  of  the  transpiring  gas  on  the  one 
hand  and  the  local  temperature  of  the  passage  walls  (i.e.  matrix  tem¬ 
perature)  on  the  other.  That  is,  in  the  majority  of  cases 


T 

ag 


=  T 

a 


and  T, 
bg 


(9) 


where,  for  equally  sized  elements,  entering  and  leaving  coolant  tem¬ 
peratures  are  the  average  internodal  temperatures  of  the  matrix 
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V  T3 
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(10) 


For  unequally  sized  elements  we  introduce  geometric  weighting  coefficients  . 
The  energy  balance  on  impermeable  wall  elements  has  the  same  form  as 
equation  (8)  with  the  coolant  term  removed. 


COOLANT  PRESSURE  AND  TEMPERATURE  IN  PLENUM  AND  PIPING 


The  continuity  condition  on  gaseous  coolant  flow  in  the  compact 
plenum  is  illustrated  with  respect  to  a  finite  control  volume  in 
Figure  4.  The  negative-blowing  mass  flux  (3)  leaving  the  side  wall 
of  cylindrical  control  volume  (4)  depletes  the  main  stream  flowing  in 
the  compact  plenum.  Letting  subscripts  i  and  j  denote  channel  flow 
conditions  at  entrance  and  exit  circular  faces  of  control  volume  (4), 
respectively,  we  have  the  finite-difference  equation 


M. 

l 


-  M. 


(11) 


Where  m  is  negative  blowing  with  respect  to  control  volume  (4).  As 
in  boundary  layer  flow,  mass  storage  in  the  control  volume  is  negligible 
for  typical  density  change  rates.  The  negative  blowing  mass  flux  enters 
porous  wall  segment  (2),  Figure  4,  in  the  pores  of  which  transpiration 
flow  is  also  essentially  steady.  In  most  systems,  transpiration  through 
the  porous  wall  is  nearly  one -dimensional ,  hence  it  suffices  to  divide 
the  segment  into  prescribed  channels  leading  from  the  compact  plenum 
to  the  external  surface.  Thus*  by  continuity  the  negative  blowing  mass 
flux  m  is  equal  to  the  sum  of  flows  to  those  channels  fed  by  the  side 
wall  of  control  volume  (4). 

The  momentum  equation  of  channel  flow  in  a  compact  plenum  with 
negative  blowing  must  include  the  effect  on  wall  friction  and  momentum, 
produced  by  motion  of  gas  toward  the  porous  wall.  Figure  4  depicts  the 
principal  terms •  The  corresponding  finite-difference  equation  is 
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The  bar  over  the  wall  friction  factor  indicates  the  averaging  of 
friction  effects  over  both  solid-wall  and  porous-wall  zones  about  the 
cylindrical  periphery  of  the  control  volume.  The  skin  friction  effect 
in  porous  wall  zones  is  physically  increased  by  the  action  of  negative 
blowing  which  causes  the  stream-wise  velocity  profile  of  turbulent 
flow  to  fill  out  near  the  wall.  That  is,  the  coefficient  f  in  equation 
(12)  is  numerically  greater  than  the  coefficient  f°  without  blowing 
and  at  the  same  channel-flow  Reynolds  number.  Experiments  to  be 
discussed  later  show  that  a  modified  form  of  an  equation  due  to 
Mickley,  et  al  (Reference  6),  namely 


f_ 

f° 


At  a  given 
channel  Re 


(13) 


is  applicable.  We  call  B£  the  blowing  parameter  for  friction  and 
define  it  as  a  measure  of  the  ratio  of  the  blowing  velocity  v^  (which 
has  a  negative  value  for  negative  blowing)  to  the  main  stream  bulk 
velocity  u^  by  the  relation 
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As  the  blowing  parameter  becomes  small,  a  linear  approximation  to 
equation  (13) 
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can  be  found.  We  would  then  say  that  the  dimensionless  friction  blowing 
factor  was  1/2  in  turbulent  channel-flow. 

Negatively  blown  coolant  gas  passes  into  the  porous  matrix  (2), 
Figure  4,  where  the  pressure-drop  associated  with  coolant  flow 
through  interconnected  pores  depends  primarily  on  the  permeability 
coefficient  of  the  matrix.  If  flow  through  interconnected  pores  is 
laminar,  incompressible,  and  isothermal,  the  coefficient  of  permeability 
is  a  constant  and  the  pressure  field  of  transpiration  flow  obeys  Darcy’s 
law.  In  practice  however,  flow  through  interconnected  pores  is 
turbulent  and  compressible.  Flow  is  also  non- iso thermal  due  to  spatial 
and  time  variations  in  temperature  associated  with  the  transient  heat 
transfer  processes.  Under  these  complex  conditions  it  is  appropriate 
to  make  use  of  the  fact  that  transpiration  flow  is  predominantly  one 
dimensional  and  to  employ  a  modified  form  of  Darcy’s  law 
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relating  local  transpiration  mass  flow  m  to  the  overall  pressure  drop 
in  a  prescribed  local  flow  direction.  The  wall  permeance  and  the 
exponent  a  are  determined  from  available  information  on  the  matrix 
in  question  including  laboratory  permeance  measurements. 


The  energy  equation  of  channel  flow  in  the  compact  plenum  accounts 
for  the  net  total-enthalpy  flow  into  a  finite  control  volume  (4), 

Figure  3,  and  heat  transfer  across  the  cylindrical  portion  of  the 
control  surface.  The  latter  is  taken  to  be  in  contact  with  the  plenum 
wall  and  bar  is  placed  over  the  Stanton  number  to  indicate  averaging 
over  the  cylindrical  portion  of  the  control  surfaces.  Thus  the 
energy  equation  for  control  volume  (4)  is 
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axial  mean 


where  denotes  the  average  enthalpy  of  negatively  blown  coolant 
leaving  the  control  volume  and  entering  the  porous  wall.  As  in  the 
case  of  the  skin  friction  coefficient  f,  the  Stanton  number  St  is 
evaluated  in  the  presence  of  negative  blowing.  Although  channel-flow 
data  are  not  available  on  the  Stanton  number  with  negative  blowing, 
Reynolds*  analogy  suggests  that  a  slightly  modified  form  of  Mickley's 
(Reference  6)  equation 
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is  applicable.  We  call  the  parameter  Bgt 
transfer  and  define  it  by  the  relation 
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the  blowing  parameter  for  heat 
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where  St°  is  the  Stanton  number  of  unblown  flow, 
which  Reynolds  analogy 


In  gas  flows  for 


(20) 
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holds  we  conclude  that  B^  and  B  ^  are  numerically  the  same.  As  noted  in 
the  discussion  of  friction,  a  linear  approximation 


St 


(21) 


is  applicable  at  low  blowing  rates.  We  would  then  say  that  the  dimension¬ 
less  heat  transfer  blowing  factor  was  1/2  in  turbulent  channel-f low . 
Because  of  Reynolds1  analogy,  the  value  is  the  same  as  for  friction. 

The  value  of  1/2  for  turbulent  channel  flow  can  be  compared  to  Eckert Ts 
(Reference  3)  recommended  value  of  0.37  for  turbulent  boundary- layer 
flow. 


COMPACT- PLENUM  FLOW  EXPERIMENTS 


Although  Mickley  et  al  (Reference  6)  experimented  with  negative 
blowing  effects  on  boundary- layer  flow  and  Tennekes  (Reference  7)  has 
recently  obtained  improved  results,  there  is  no.  literature  dealing  with 
negative  blowing  effects  on  the  skin  friction  and  heat  transfer  of 
channel  flows.  For  this  reason  a  somewhat  idealized  but  basic  compact- 
plenum  bench  test  apparatus  was  devised.  As  shown  in  Figure  5  in  part 
(a)  the  basic  apparatus  was  a  porous  tube  test  section.  The  tube  was 
0.4  inch  o.d.,  0.160  i.d.  porous  sintered  stainless  steel,  of  1  to 
1  1/2  micron  pore  size  and  15%  porosity.  The  tube  shown  in  part  (a) 
of  Figure  5  was  about  8  inches  long.  It  was  supplied  from  the  left 
with  metered,  clean,  dry,  room- temperature  nitrogen  gas.  A  sliding 
pressure  probe  on  the  axis  was  used  to  determine  the  streamwise 
static  pressure  distribution.  A  sliding  collection  chamber,  part  (b) 
of  Figure  5,  was  used  to  measure  the  streamwise  variation  in  the 
negative  blowing  rate. 

Tests  were  first  conducted  with  the  tube  wrapped  to  prevent 
negative  blowing.  Basic  no-blowing  skin  friction  coefficients  f°  were 
determined  by  computing  the  net  pressure  force  and  momentum  change. 

The  resulting  plot  of  f°  versus  Reynolds  number  indicated  friction 
was  about  10%  above  the  smooth  annulus  f°  value.  This  difference 
was  attributed  to  surface  roughness. 

Tests  were  next  conducted  with  the  wrapping  removed.  The  measured 
streamwise  variations  of  negative  blowing  rate  and  static  pressure  were 
used  in  calculations  of  the  coefficient  of  skin  friction  f  at  ten 
local  points  along  the  flow  axis.  The  local  skin  friction  ratio  f/f° 
was  then  formed  and  plotted  against  each  of  three  different  negative 
blowing  (suction)  parameters.  The  three  parameters  were  taken  from 
the  boundary  layer  literature  and  modified  by  replacing  free-stream 
velocity  with  bulk  velocity.  The  three  parameters  are 
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(Olsen  and  Eckert  (Reference  8)  positive  blowing 
pipe-flow  experiments) 

(Tewfik  (Reference  9),  boundary  layer  experiments, 
and  Rannie  (Reference  10),  sublayer  theory) 

(Mickley  et  al  (Reference  6),  boundary  layer  experi¬ 
ments) 

8  show  a  progressive  improvement  in  correlation  of 
the  results  as  f°  plays  a  stronger  role  in  the  blowing  parameter. 

Figure  8  shows  how  well  the  modified  film  theory  equation  of  Mickley 
fits  the  results  of  these  tests.  For  engineering  purposes,  we  con¬ 
clude  that  negative  blowing  effects  on  axisymme trie  channel  flow  are 
covered  by  the  modified  Mickley  equation 


It  has  already  been  pointed  out  that  equation  (22)  reduces  to  the 
linear  form 


w 


Figures  6,  7,  and 


=  1  -  7  B  (23) 

small  B^ 

at  low  negative  blowing  rates.  This  relation  is  obtained  by  differ¬ 
entiating  equation  (22)  to  find  the  slope  at  B^  =  0.  Large  negative 
blowing  rates  cause  the  exponential  term  to  drop  out  of  the  denominator 
of  equation  (22)  thus  f/f°  is  asymptotic  to  -  B^,  i.e. 


Large  negative  B^ 


(24) 


This  also  follows  if  we  accept  the  premise,  based  on  turbulent  boundary 
layer  experimental  evidence,  that  large  negative  blowing  rates  cause 
the  u-velocity  profile  of  two-dimensional  channel  flow  to  approach 
the  slug  profile  shape.  For  such  an  asymptotic  channel-flow  velocity 
profile  shape,  the  sidewise  momentum  transport  to  the  channel  wall  is 


-pv^u^  per  unit  wall  area, 
defined  by  the  relation 


The  apparent  shearing  stress  T  is 


w.a 
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T 

w,a 


£  K 

2 


(25) 


Assuming  that  apparent  shear  overshadows  viscous  shear  we  find  the  skin 
friction  coefficient  with  negative  blowing, 


% 


(26) 


Thus 


^  Asymptotic 
slug  profile 


2v 

w 


(27) 


where  the  right  hand  side  is  -  B^,  by  definition. 

The  heat  transfer  relation  follows  from  the  above  discussion  by 
establishing  the  asymptote  based  on  sidewise  energy  transport 
-pv  c  (Tu-T  )  per  unit  wall  area  and  the  definition 

W  p  D  W 

-  st  (fVp)  <T„  -  V  <28> 


By  equating  the  energy  transport  to  and  dividing  this  by  St°  we  find 


Asymptotic 
slug  profile 


v 

w 


(29) 


By  similar  arguments,  the  overall  trend  of  channel  flow  experiments  can 
be  expected  to  follow  the  prediction  of  Mickley  et  al  (Reference  6)  as 
modified  for  channel  flow,  i.e. 


BSt 

e  -  1 


(30) 


The  experimental  friction  apparatus  can  be  instrumented  as  shown  in 
parts  c  and  d  of  Figure  5  to  test  this  theory  on  an  electrically  heated 
apparatus.  Such  work  is  planned  for  the  future. 
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At  present  there  are  no  test  data  available  by  means  of  which  the 
entire  computer  program  might  be  checked  out.  It  is  essential  that 
this  be  done  and  future  work  is  aimed  at  that  goal.  Meanwhile,  as 
each  portion  of  the  program  has  been  completed,  its  successful  operation 
has  been  proven  out  by  check  runs  confined  to  the  portion  of  the 
program  in  question. 

The  first  example  of  program  check  out  is  the  computation  of 
transient  temperatures  in  a  transpiration  cooled  porous  matrix.  An 
initially  isothermal  transpiring  plane  porous  slab,  Figure  9,  is 
supplied  at  one  face  with  coolant  at  a  fixed  temperature.  The  sudden 
application  of  convective  heating  at  the  other  face  causes  the  outer 
face  temperature  to  rise  and  reach  an  equilibrium  value,  Figure  10. 
Temperature  rise  of  the  outer  face  is  plotted  against  the  Fourier 
number  of  the  slab  for  one  value  of  the  Biot  number  and  two  values  of 
the  transpiration  parameter  g.  Blotted  points  portray  computer  results 
which,  because  of  the  actual  temperatures  used,  reflect  the  reference 
temperature  effect  on  the  otherwise  constant  convective  coefficient. 

When  Schneider's  (Reference  11)  theoretical  results  are  adjusted  down¬ 
ward  to  reflect  the  changing  Biot  number  the  agreement  is  seen  to  be 
very  good.  On  this  basis  the  finite-difference  energy  portions  of  the 
program  involving  transpiration  effects  were  judged  to  be  working 
satisfactorily  0 

The  interaction  of  the  coolant  continuity  equation,  the  momentum 
equation,  and  the  modified  version  Darcy's  law  were  checked  out  in  a 
computer  run  utilizing  actual  test  data  from  the  compact  plenum 
apparatus  previously  described.  In  this  run,  the  known  local  permeance 
characteristics  of  the  wall  of  the  porous  tube  were  fed  into  the 
program  along  with  the  modified  Mickley  equation  governing  the  negative 
blowing  effect.  The  program  run  then  computed  the  negative  blowing 
(radial  mass  flow  m)  from  Darcy's  law  and  the  static  pressure  dis¬ 
tribution  in  the  compact  plenum  chamber.  In  Figure  11,  the  close 
agreement  between  computed  radial  mass  flows  (negative  blowing  rates) 
and  measured  values  is  due  to  the  fact  that  input  permeance  values 
were  based  on  the  test  measurements.  The  close  prediction  of  static 
pressure  distributions  in  the  compact  plenum  is  an  indication  of  the 
reliability  of  the  modified  Mickley  equation  and  of  the  reliability 
of  the  finite  differencing  and  iterative  schemes  involved  in  the 
channel-flow  portions  of  the  program*. 
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The  bulk  temperature  of  the  coolant  main  stream  in  a  compact  plenum 
cannot  be  measured  locally  without  upsetting  conditions  further  down¬ 
stream.  However  some  unreported  external  porous  surface  temperatures 
have  been  measured  in  the  laboratory  and  calculated  by  means  of  the 
computer  program.  The  results  indicate  that  the  computer  program  is 
performing  consistently  and  as  designed.  Discrepancies  between  computed 
and  measured  values  are  undoubtedly  due  to  the  difficulty  of  developing 
meaningful  surface  temperature  measurement  techniques  for  porous  wall 
models.  More  intensive  work  in  this  area  is  planned  for  the  future. 

Film  cooling  computer- prog ram  procedures  have  been  checked  out 
against  the  downstream  measurements  of  Goldstein,  Shavit  and  Chen 
(Reference  4)  using  the  effectiveness  equation  of  Kutateladze  and 
Leont'ev  (Reference  5)  and  the  finite-difference  model.  Results  for 
the  model  shown  in  Figure  12  are  presented  in  Figure  13  and  indicate 
satisfactory  performance  of  the  computer  program. 


FUTURE  PLANS 


In  order  to  prove  out  the  computer  program,  an  instrumented  wind- 
tunnel  test  model  of  a  compact  system  will  be  designed  and  tested 
under  controlled  conditions.  The  controlled  conditions  will  correspond 
to  those  described  under  ’’Problem  Definition”.  They  will  provide  the 
necessary  inputs  to  the  computer  program  along  with  model  specifications. 
Instrumentation  will  give  readings  of  key  body  temperatures  and 
coolant  pressures  against  which  computed  values  will  be  compared. 
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(U)  Fig.  1  A  LEADING-EDGE  COMPACT-PLENUM  INSTALLATION  IN  FLIGHT 
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(U)  Fig.  2  BASIC  ELEMENTS  OF  A  COMPACT  GAS-TRANSPIRATION  SYSTEM 
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(U)  Fig.  3  BASIC  APPROACH  TO  ENERGY  RELATIONS 
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(U)  Fig.  4  BASIC  APPROACH  TO  COOLANT  CONTINUITYAND  MOMENTUM  RELATIONS 
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(U)  Fig.  5  TEST  ARRANGEMENTS  FOR  TRANSPIRATION  EXPERIMENTS 
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(U)  Fig.  7  LOCAL  SUCTION  RATE  PARAMETER  NO.  2 
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(U)  Fig.  8  LOCAL  SUCTION  RATE  PARAMETER  NO.  3 
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(U)  Fig.  10  TEMPERATURE  RESPONSE  AT  x  =  0  OF  A  POROUS  PLATE,  0  <  x  <  6,  AFTER 
EXPOSURE  TO  A  UNIFORM  TEMPERATURE  CONVECTIVE  ENVIRONMENT  Tc 
AT  x  =  0  AND  COOLED  BY  A  STEADY  FLOW  OF  TRANSPIRATION  FLUID 
THROUGH  THE  PLATE  ENTERING  AT  x  =  8  WITH  T  =  T0  AT  x  »  8,  FOR 
VARYING  BIOT  NUMBERS  (Bi)  AND  COOLANT  FLOW  PARAMETERS  (g) 
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(U)  Fig.  11  COMPARISON  OF  THEORETICAL  PREDICTIONS  AND  EXPERIMENTAL 
RESULTS  OF  PRESSURE  AND  RADIAL  MASS  FLOW  IN  A  POROUS  TUBE 
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(U)  Fig.  12  EXPERIMENTAL  MODEL  USED  TO  CHECK  FILM  COOLING 
COMPUTER  PROGRAM 
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Paper  No.  6 

EXPERIMENTAL  INVESTIGATION  OF  TRANSPIRATION 
COOLING  NEAR  THE  STAGNATION  POINT  OF  A  CYLINDER 

(U) 

(Paper  UNCLASSIFIED) 
by 


Richard  L.  Humphrey 
U.S.  Naval  Ordnance  Laboratory 
White  Oak,  Silver  Spring,  Md.  20910 


ABSTRACT.  (U)  A  series  of  tests  were  made  with  a 
cylindrical  model  equipped  with  transpiration  cooling  over 
a  region  ±15  degrees  from  the  stagnation  line.  The  model 
was  instrumented  with  thermocouples  to  measure  wall  and 
reservoir  temperatures.  Air  was  used  as  a  cooling  medium. 
The  tests  were  made  at  a  nominal  Mach  number  of  6.5  and 
stagnation  point  gas  temperatures  of  2500  to  3500°K  and  a 
Pitot  pressure  of  0.14  atm.  The  dimensionless  temperatures 
agree  within  5  to  7  percent  with  predictions  while  heat- 
transfer  rates  agree  within  15  percent. 

(U)  The  effect  of  the  external  pressure  distribution 
on  the  local  coolant  flow  rate  and  wall  temperature  is 
discussed.  Varying  wall  thickness  and  permeability  may  not 
be  the  whole  answer. 
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NOMENCLATURE 

C  =  velocity  gradient 

0^  =  specific  heat  at  constant  pressure 

H  =  enthalpy 

m  =  mass  flow  of  coolant 

N  =  Nusselt  number 

u 

Pr  =  Prandtl  number 
q  =  heat  transfer  rate 
Re  =  Reynolds  number 
t  =  time 

T  =  temperature 

W  =  weight  of  matrix  along  coolant  path 
0  =  angle  from  stagnation  line 

y  =  viscosity 

p  =  density 

Subscripts 

c  =  cavity  (reservoir) 
e  =  edge  of  boundary  layer 
k  =  side  wall  conduction 
m  =  matrix 
o  =  stagnation  point 
a  =  ambient 
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oc 


stagnation  point  cooled 


ouc 


stagnation  point  uncooled 


sh 


shielding  from  transpiration 


t^& =  time  1  and  time  2 


w 


2 


outside  wall  of  porous  section 
inside  wall  of  porous  section 


INTRODUCTION 


(U)  Transpiration  cooling  is  an  attractive  means  of 
providing  protection  for  surfaces  subjected  to  moderate  to 
high  heat-transfer  rates  for  long  periods  of  time.  This 
would  be  the  case  for  a  ramjet  vehicle  or  a  hypersonic 
transport  operating  at  normal  flight  altitudes.  Highest  heat 
loads  are  encountered  on  stagnation  lines  such  as  engine 
air  inlets,  aerodynamic  surface  leading  edges  and  body 
noses.  Most  of  these  surfaces  are  two  dimensional  and 
close  to  cylindrical. 

(U)  When  the  present  studies  were  initiated,  a 
considerable  amount  of  information  was  available  on  flat- 
plate  and  three-dimensional  (conical  and  spherical) 
transpirational  cooling.  However,  two-dimensional  data 
particularly  near  the  stagnation  line  were  missing.  Most 
leading  edges  proposed  for  aircraft  wings  and  air  inlets 
are  very  small  in  radius,  in  the  case  of  ramjet  air  inlets 
as  small  as  0.050  inch.  This  led  us  to  attempt  as  small 
a  model  as  was  practicable  to  manufacture  and  instrument. 

It  was  recognized  that  the  strong  pressure  gradients  over 
the  cylindrical  surface  would  result  in  a  variation  of  mass 
injection  rate  and  wall  temperatures  near  the  stagnation 
line  and  the  design  was  modified  to  minimize  these  effects. 

DESCRIPTION  OF  MODELS 

(U)  Considerations  of  available  materials,  ease  of 
instrumentation,  tunnel  blockage  requirements,  and  wind 
tunnel  conditions  led  us  to  choose  a  cylinder  0.500  inch 
in  diameter  as  the  basis  of  the  model.  The  model  design  is 
detailed  in  Fig.  1.  The  model  consists  of  a  holder  for 
mounting  in  the  tunnel-injector  mechanism,  a  gas  reservoir, 
a  porous  tube  holder,  the  porous  tube  and  an  external  steel 


sheath . 
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(U)  The  porous  tube  is  nominally  one-half  inch  in 
diameter  with  0.125-inch  walls  and  is  sintered  from  10  or  20 
microns  type  316  stainless-steel  spheres.  This  size  wall 
thickness,  material,  and  particle  size  are  standard  filter 
elements  manufactured  by  the  Mott  Metallurgical  Company. 

The  porous  tube  is  held  in  a  cradle  and  the  two  parts  are 
mounted  in  a  holder  that  covers  the  ends  of  the  tube  and  is 
as  thick  as  the  tube  diameter. 

(U)  The  exterior  of  the  model  is  covered  by  a  sheet  of 
0.003-inch  type  304  stainless  steel.  The  sheet  has  a  cutout 
two-inches  long  by  0.125-inch  wide  along  the  stagnation  line 
of  the  model.  The  sheet  is  electron  beam  welded  airtight 
along  all  exposed  edges  and  the  cutout.  Around  the  cutout 
the  sheet  is  welded  directly  to  the  porous  tube.  Thus  the 
transpiring  fluid  is  ejected  from  the  model  only  along  the 
stagnation  line  and  ±15°  from  it. 

(U)  The  coolant  fluid  is  brought  into  the  model 
through  a  single  tube  opening  into  a  rear  reservoir.  It 
passes  through  a  series  of  connecting  passages  to  the 
cavity  inside  the  porous  tube.  From  here  it  filters  through 
the  porous  tube  and  out  the  cutout  along  the  stagnation 
line.  The  mass  flow  of  coolant  air  (m)  was  measured  with  a 
variable  area  glass -tube  flowrater.  The  air  was  obtained 
from  the  house  air  system  (100  psi;  -50°F  to  -140°F  dew¬ 
point)  ,  filtered  and  reduced  to  the  desired  pressure  with 
a  regulator.  The  flow  rate  of  coolant  was  independent  of 
the  external  pressure  on  the  model  over  the  range  of  400  y 
Hg  to  200  mmHg. 

(U)  Three  thermocouples  are  mounted  in  the  model. 

The  chromel-alumel  thermocouples  are  of  number  36  (0.005  inch) 

wire  and  are  located  on  the  outside  wall  (Tw) ,  inside  wall 
(T2)  and  in  the  cavity  (Tc) .  The  outside  wall  thermocouple 
passes  through  the  wall  in  two  insulator  tubes  0.125  inch 
apart  and  the  bead  is  formed  between  them.  The  thermocouple 
bead  is  in  a  region  where  the  flow  of  the  coolant  is 
undisturbed  and  thus  is  more  likely  to  indicate  the  actual 
cooled  wall  temperature. 

TEST  PROCEDURES 

(U)  The  tests  were  conducted  in  the  Naval  Ordnance 
Laboratory  (NOL)  3  Megawatt  Arc  Tunnel.  The  nominal  Mach  6.4 
nozzle  was  used.  When  tunnel  pressures  had  stabilized,  a 
quick  injection  mechanism  with  a  transit  time  of  150  ms  was 
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used  to  insert  the  model  into  the  gas  stream.  The  model 
remained  in  the  stream  for  five  or  seven  seconds.  Coolant 
flow  in  the  model  was  established  before  starting  the  tunnel. 

(U)  Conditions  for  the  runs  with  the  first  model  are 
given  in  Table  1.  Arc  heater  and  Pitot  pressures  remain 
very  stable  during  tunnel  operation  but  the  heater  tempera¬ 
ture  reaches  a  steady  state  slowly.  Thus  there  is  a  variation 
of  stagnation  point  temperature  and  Mach  number  during  the 
time  the  model  is  in  the  hot  gas.  This  is  reflected  in  the 
steadily  increasing  wall  temperatures  recorded  on  the 
models,  for  example  Fig.  2.  However,  the  model  does  achieve 
a  steady  state  with  respect  to  the  stagnation  point  tempera¬ 
ture  . 


(U)  After  a  run,  the  surface  of  the  model  was  found  to 
be  coated  with  a  fine  oxide  deposit.  This  deposit  did  not 
appear  to  cause  any  nonuniformity  in  the  coolant  distribution. 
This  is  probably  because  the  pressure  drop  across  the  wall 
is  large  compared  to  the  external  pressure.  Wiping  the 
surface  with  solvent  restored  the  appearance. 

(U)  At  the  lowest  mass  flow  of  coolant  the  stainless 
steel  sheath  failed  at  the  stagnation  line  near  the  ends. 

Since  the  porous  window  does  not  extend  to  the  ends  of  the 
model,  this  region  is  uncooled  except  by  internal  conduction. 
The  damaged  area  appears  to  have  ruptured  outward  from  the 
pressure  inside  the  wall.  Examination  of  the  model  revealed 
that  the  sheath  had  remained  welded  to  the  porous  tube  all 
around  the  window.  All  three  thermocouples  and  the  Pitot 
tube  were  still  working  and  the  porous  surface  was  unharmed. 

RESULTS  AND  ANALYSIS 

Temperatures 

(U)  Figure  2  shows  typical  temperature  history  results 
for  the  transpiration  cooled  models.  The  temperatures  show 
a  steady  increase  during  the  test.  After  the  first  three 
seconds  they  follow  the  trend  of  the  stagnation  point 
temperatures.  An  indication  of  steady-state  temperatures 
is  the  difference  between  Tw  and  T2 ,  the  gradient  across 
the  porous  wall.  After  1.5  seconds  Tw  -  T2  does  not  change 
significantly.  The  thermal  conductivity  of  the  porous  wall 
material  is  only  a  weak  function  of  temperature  and  the 
constant  temperature  drop  may  be  taken  as  an  indication  of 
steady-state  cooling. 
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(U)  The  dimensionless  temperature  function  (Te-Tw)/ 
(Te-Tc)  given  in  Fig.  3  shows  the  rapid  approach  to  steady- 
state  temperatures.  Here  the  values  of  wall  temperatures 
predicted  from  the  tunnel  flow  conditions  are  shown  for  the 
three  runs  with  the  first  model.  While  the  predicted  and 
experimental  values  do  not  agree  exactly,  after  about  three 
seconds  the  trend  of  the  data  is  in  good  agreement. 

(U)  The  analysis  of  Roberts-*-  is  followed  for  the 
calculation  of  the  predicted  values  of  the  runs.  A  somewhat 
simplified  treatment  of  Roberts'  material  is  given  by 
Truitt^.  Roberts  looks  at  mass  transfer  cooling  near  the 
stagnation  point  and  considers  both  two- and  three-dimensional 
flow  situations.  It  is  assumed  that  the  velocity  parallel 
to  the  wall  is  a  linear  function  of  the  distance  along  the 
wall  from  the  stagnation  point.  The  temperature,  velocity, 
and  concentration  boundary  layers  are  assumed  to  be  of 
constant,  although  different,  thicknesses .  Gradients  in 
these  layers  are  taken  to  be  a  function  of  the  distance  from 
the  wall  only  and  the  coolant  is  assumed  to  be  injected  at 
equal  velocity  along  the  wall.  It  is  considered  that  the 
boundary  layer  will  increase  in  thickness  because  of  the 
injection;  but  the  analysis  implies  that  the  extra  thickness 
will  be  small  enough  to  avoid  disturbance  of  the  flow  field 
around  the  body. 

(U)  The  dimensionless  temperature  function  for  a  two- 
dimensional  body  is  found  by  Roberts  to  be 


T  -T 
e  w 


where 


0.4 


P 


0.4 


rw 


The  various  quantities  are  evaluated  at  the  cooled 
wall  temperature  and  at  the  temperature  behind  the 
normal  shock  in  front  of  the  body.  In  so  far  as 
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possible  real  gas  thermodynamic  and  transport  properties 
were  used  in  the  evaluation.  The  coolant  flow  rate  was 
assumed  to  be  uniform  over  the  porous  surface. 

(U)  The  dimensionless  temperature  function  as  calculated 
for  the  three  runs  with  the  first  model  is  shown  on  Fig.  3. 
With  the  exception  of  the  first  run,  the  predicted  values  of 
wall  temperature  are  lower  than  the  experimental  data.  The 
calculated  temperature  function  for  runs  2  and  3  agree  with 
the  experimental  values  within  5  to  7  percent.  The  different 
behavior  of  the  first  run  is  believed  due  to  an  error  in  the 
coolant  mass  flow  or  to  some  error  in  tunnel  flow  conditions 
but  has  not  been  found  at  this  time. 

(U)  The  dimensionless  temperature  function  as  a  function 
of  the  mass  flow  parameter  m/(pw  yw  c) is  given  in  Fig.  4. 
Both  experimental  and  calculated  values  are  shown.  Again 
with  the  exception  of  the  first  run,  agreement  is  within  5 
to  7  percent.  The  range  of  mass  flows  covered  is  not  large 
but  represents  relatively  large  absolute  values. 

(U)  The  heat-transfer  parameter  Nuw/Rew^-/2  as  a  function 
of  the  relative  mass  rate  is  given  in  Fig.  5.  The  values 
calculated  from  the  experimental  results  are  compared  with 
curves  for  Pr  =  0.7  and  1.0  taken  from  Roberts'  report. 

The  Pr  of  the  experiments  is  about  0.75.  From  this  it  seems 
that  the  heat-transfer  function  at  these  large  mass  rates  is 
either  constant  or  is  increasing  with  increasing  mass 
injection.  These  large  mass  rates  may  disturb  the  flow 
field  around  the  body  and  may  result  in  the  flow  near  the 
stagnation  point  becoming  turbulent  (ref.  3  for  instance). 

(U)  The  analysis  of  Roberts  does  not  tell  us  what  is 
happening  inside  the  wall  but  predicts  only  the  outer  wall 
temperature.  The  temperature  profiles  within  the  wall  are 
vital  to  the  proper  choice  of  material,  wall  thickness, 
etc.,  to  allow  the  wall  to  sustain  the  aerodynamic  and 
mechanical  loads  while  meeting  the  cooling  requirements. 

For  these  experiments  it  is  important  to  know  if  the  coolant 
and  matrix  or  wall  temperatures  are  equal  at  any  location  in 
the  wall.  If  this  is  true  the  wall  thermocouples  will  be 
measuring  both  the  coolant  and  matrix  temperatures  correctly. 

(U)  Koh  and  del  Casal^  discuss  the  processes  occurring 
within  the  wall  and  develop  expressions  for  determining 
coolant  and  matrix  temperature  and  pressure  profiles. 

These  are  based  on  the  coolant  mass  flow  and  the  coolant 
matrix  properties.  For  the  model,  coolant  and  flow  rates 
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used  in  these  tests  and  using  their  relations  the  coolant 
and  matrix  temperatures  are  found  to  be  equal  at  every 
point  through  the  wall. 

Heat  Transfer 


(U)  The  steady-state  heat  transfer  to  the  model  consists 
basically  of  two  components,  q>  -  the  heat  transfer  blocked 
by  transpiration,  and  qw  -  the  heat  transferred  to  the  wall 
and  raising  the  coolant  temperature.  Since  the  tunnel 
conditions  are  changing  during  the  tests  there,  is  a  term 
qm  “  the  rate  of  heat  storage  in  the  matrix  as  the  wall 
temperatures  increase.  So  the  total  stagnation  point  heat 
transfer  to  the  model  is 


and 


^oc  ^sh  +  %  +  ^m 

q  .  =  (1  -  4  p  _0’6)  (H  -  II  )  m 

sh  3  rw  e  w 

%  “  (Hw  -  V* 

(Tt2-Ttl> 

^m  m  Cpm  At 


(U)  A  breakdown  of  the  average  experimental  heat-transfer 
rates  for  the  runs  is  given  in  Table  2.  The  largest  item  is 
qsh  which  amounts  to  70  to  80  percent  of  the  total. 

Eighteen  to  28  percent  of  the  total  reaches  the  wall  raising 
the  coolant  temperature.  The  small  remainder  is  the  heat 
storage  in  the  matrix.  In  a  steady  state  this  storage  term,  of 
course,  would  vanish. 

(U)  The  stagnation  point  heat  transfer  to  an  uncooled 
wall  was  calculated  using  standard  methods.  With  the  excep¬ 
tion  of  the  maverick  first  run  the  experimental  and  calcu¬ 
lated  stagnation  point  heat  transfer  rates  agree  within 
15  percent.  This  is  consistent  with  previous  calorimeter 
tests  made  in  the  tunnel. 

(U)  The  last  point  in  the  heat- transfer  discussion  is 
the  term  q^  listed  in  the  table.  This  term  is  given  by 


^k 


=  m(H  -  II  ) 
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and  gives  the  heating  of  the  coolant  fluid  as  it  passes 
through  the  model  to  the  matrix.  This  term  while  small  is 
not  negligible.  The  body  of  the  model  is  not  cooled  and 
qj,  gives  an  indication  of  the  heat  transferred  to  the 
cavity  through  the  uncooled  sidewalls. 

(U)  The  agreement  between  the  uncomplicated  theory 
of  Roberts  and  the  experiments  is  good.  The  dimensionless 
wall  temperatures  could  be  calculated  within  5  to  7  percent 
of  the  experimental  results.  Agreement  of  the  experimental 
and  calculated  stagnation  point  heat-transf er  rates  is 
within  normal  limits. 

(U)  The  conditions  of  the  tests  and  the  design  of  the 
model  were  such  as  to  minimize  problems  that  must  be  faced 
for  flight  vehicles.  The  most  important  of  these  is 
"coolant  starvation". 

Coolant  Starvation 

(U)  The  external  pressure  distribution  around  a  body 
near  the  nose  of  an  axisymmetric  shape  or  the  leading  edge 
on  a  two-dimensional  body  may  have  a  profound  effect  on  the 
local  coolant  flow  rate  near  the  nose.  If  the  porous  wall 
is  of  uniform  thickness  and  permeability  and  the  plenum  or 
reservoir  pressure  is  constant  the  coolant  flow  rate  will 
vary  with  changes  in  the  external  pressure.  At  the  stagnation 
point  the  flow  for  these  conditions  would  be  a  minimum  and 
would  increase  away  from  the  stagnation  point  as  the  external 
pressure  drops. 

(U)  This  change  in  coolant  flow  is  dramatically  shown 
in  Fig.  6.  Two  flight  conditions  and  the  model  test 
environment  are  compared.  In  each  case  the  wall  thickness 
and  permeability  were  uniform  and  the  internal  pressure  was 
chosen  to  provide  the  desired  coolant  flow  at  the  stagnation 
point.  Away  from  the  stagnation  point  a  one-dimensional 
flow  analysis  with  a  form  of  the'  Darcy  pressure  drop 
equation4  gave  the  local  coolant  flow.  In  the  model  tests 
the  ratio  of  cavity  pressure  to  stagnation- point  pressure 
is  high  while  in  the  sea  level  flight  case  it  approaches 
unity.  The  model  will  have  only  a  25  percent  increase  in 
coolant  flow  15  degrees  from  the  stagnation  line  while  the 
sea  level  flight  coolant  flow  will  increase  to  six  times  the 
stagnation-point  value  only  5  degrees  from  the  stagnation 
line . 
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(U)  Recently  Schneider  '  and  others  have  described 
this  flow  variation  in  some  detail  and  discuss  the  under¬ 
cooling,  over-cooling,  and  coolant  starvation  that  will 
result  if  average  coolant  flow  rates  are  used.  Usually  the 
stagnation  region  would  be  under-cooled  or  have  "coolant 
starvation"  while  regions  away  from  the  stagnation  point 
will  be  over-cooled  or  have  more  coolant  than  needed.  As 
the  internal  pressure  approaches  the  stagnation  point 
pressure  the  external  flow  may  be  ingested  into  the 
stagnation  region. 

(U)  Schneider  discusses  the  effects  of  varying 
permeability  and  wall  thickness  to  control  the  distribution 
of  coolant  over  the  body.  However,  if  the  body  is  small, 
as  the  leading  edge  of  a  wing,  it  may  not  be  practicable 
to  vary  thickness  or  permeability  sufficiently  to  obtain 
uniform  cooling.  In  this  case  the  system  would  need  to  be 
designed  to  confine  the  injection  to  an  area  close  to  the 
stagnation  line.  The  coolant  flow  rate  would  have  to  be 
adequate  at  the  stagnation  line  and  still  allow  for  the 
extra  flow  away  from  that  point. 

(U)  In  the  model  tests  described  here  the  variation 
in  coolant  flow  is  not  large  but  does  cause  the  wall 
temperatures  to  change.  Figure  7  gives  an  estimate  of  the 
variation  in  temperature  over  the  cooled  portion  of  the 
model.  At  the  edge  of  the  cooled  portion  the  wall  tempera¬ 
ture  is  70  percent  of  the  stagnation  point  value. 

SUMMARY  AND  CONCLUSIONS 

(U)  Agreement  between  the  experiments  and  the  analysis 
of  Roberts  is  good  even  with  the  very  high  coolant  mass 
rates  and  the  low  wall  temperature  to  stagnation-point 
temperature  ratio.  Compensation  for  the  variation  in 
coolant  flow  across  the  porous  section  would  improve  the 
agreement . 

(U)  In  a  flight  vehicle  with  a  small  diameter  leading 
edge  it  may  not  be  feasible  to  adjust  the  wall  thickness  and 
permeability  to  allow  the  desired  coolant  distribution. 

In  such  a  case  it  would  be  essential  that  the  irregular 
distribution  be  considered  in  order  that  the  stagnation 
point  not  be  starved.  An  alternative  would  be  to  go  to  a 
slot  injection  scheme. 
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(U)  The  uncomplicated  approach  of  Roberts  can  be 
expected  to  give,  satisfactory  results  if  proper  consideration 
is  given  to  the  effect  of  external  pressure  distribution 
on  coolant  flow  rates. 
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(U)  TABLE.  1 


Test  Conditions  for  Runs  of  Model  No.  1 


Item 

Run  No.  1 

Run  No.  2 

Run  No.  3 

Heater  Temp 
(°K) 

3340-3880 

2500-3700 

2600-3550 

(°R) 

6012-6984 

4500-6660 

4680-6390 

Heater  Press 
(atm) 

21.41 

20 . 71 

20.84 

Model 

Stagnation  Point 

Temp 
( °K) 

2950-3250 

2400-3180 

2470-3090 

( °  R) 

5310-5850 

4320-5724 

4446-5562 

Model  Stag 

Point  Press 
(atm) 

0.142 

0.140 

0.139 

Mach  Number 

7.1-6.65 

7 .69-6 . 74 

7 . 65-6 . 91 

Model  Immersion 

Time  (sec) 

5.03 

7.05 

7 . 01 

Coolant  Mass 

Flow  (lb/sec) 

2.42xl0_4 

3. 07xl0"4 

2. 12xl0-4 
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(U)  FIG.  1  Details  of  Two-Dimensional  Stagnation  Line  Transpiration 
Cooled  Model 
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(U)  FIG.  2  Typical  Model  Temperatures  and  Temperature 
Differences  during  Test  Run 
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(U)  FIG.  3  Dimensionless  Wall  Temperatures  Compared  with  Predicted  Values 
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(U)  FIG.  4  Dimensionless  Wall  Temperatures  for  the  Mass  Flow  Function, 
Experimental  and  Predicted 
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(U)  FIG.  5  Comparison  of  Nusselt- Reynolds  Number  Function  from 
Test  Results  with  Results  of  Roberts 
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(U)  FIG.  6  Estimated  Effect  of  External  Pressure  Distribution  on 
Local  Mass  Injection  Rate  for  Three  Cases 
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(U)  FIG.  7  Estimated  Effect  of  External  Pressure  Distribution  on 
Wall  Temperatures  for  Model  Tests 


200 


8th  Navy  Symposium  on  Aeroballistics 


Vol.  1 


Paper  No.  7 

BALLISTICS  RANGE  EXPERIMENTS  ON  THE  EFFECT  OF 
UNIT  REYNOLDS  NUMBER  ON  BOUNDARY-LAYER  TRANSITION 

(U) 

(Paper  UNCLASSIFIED) 
by 

Norman  W.  Sheetz,  Jr. 

U.S.  Naval  Ordnance  Laboratory 
White  Oak,  Silver  Spring,  Md.  20910 


ABSTRACT.  (U)  It  has  been  observed  in  numerous  experi¬ 
mental  investigations  that  the  unit  Reynolds  number  appears 
to  have  a  significant  effect  on  the  boundary- layer  transition 
Reynolds  number.  Many  wind-tunnel  investigations  have  indi¬ 
cated  that,  for  an  increase  in  the  unit  Reynolds  number  of 
a  factor  of  10,  the  transition  Reynolds  number  can  increase 
by  a  factor  of  1.5  to  4,  depending  upon  the  particular 
experiment.  Pate  and  Schueler  have  published  a  correlation 
that  attempts  to  explain  this  effect  in  terms  of  the  turbu¬ 
lent  boundary  layer  that  existed  on  the  wind-tunnel  walls 
and  the  aerodynamic  noise  that  is  radiated  into  the  test 
section  from  these  boundary  layers.  However,  tests  have 
been  performed  by  Potter  in  a  ballistics  range  which  also 
produced  a  unit  Reynolds  number  effect  that  cannot  be 
explained  by  an  analysis  such  as  that  of  Pate  and  Schueler 
which  would  predict  no  effect  in  a  ballistics  range. 

(U)  In  order  to  help  define  the  effect  of  varying  the 
unit  Reynolds  number  on  transition,  a  test  is  currently  in 
progress  in  the  NOL  ballistics  range.  To  date,  these  results 
indicate  the  existence  of  a  unit  Reynolds  number  effect,  but 
suggest  that  it  may  not  be  as  strong  as  observed  in  many 
previous  investigations. 
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INTRODUCTION 

(U)  The  importance  of  boundary-layer  transitions  in 
many  aerodynamic  phenomena  is  well  known.  These  include 
as  some  of  the  more  obvious  the  viscous  drag,  aerodynamic 
heating  rates,  location  of  boundary- layer  separation,  and 
the  base  flow  characteristics.  It  is  also  generally 
accepted  that  a  large  number  of  parameters  affect  boundary- 
layer  transition.  Some  of  these  parameters  include  Mach 
number,  surface  roughness,  heat-transfer  rate,  mass  injec¬ 
tion  of  ablation  rate,  pressure  gradient,  nose  bluntness, 
cone  angle,  and  unit  Reynolds  number.  In  order  to  predict 
conditions  under  which  transition  will  occur,  it  is  first 
necessary  to  understand  the  role  that  each  of  the  above- 
mentioned  parameters  play  in  controlling  transition. 

(U)  Judging  by  the  volume  of  work  that  has  gone  into 
examining  boundary- layer  transition,  Morkovin  lists  265 
references  in  a  recent  comprehensive  survey  article  (Ref.  1), 
one  might  suspect  that  the  problem  has  been  solved.  Quite 
the  contrary!  The  contradicting  conclusions  that  can  be 
obtained  from  many  of  the  investigations  have  often  tended 
to  complicate  the  picture,  rather  than  to  enlighten.  In 
many  cases,  this  has  been  caused  by  experiments  in  which  a 
number  of  parameters  were  varied  in  such  a  manner  that  it 
was  not  possible  to  separate  and  define  the  effect  of  the 
individual  parameters.  It  is  the  purpose  of  the  present 
paper  to  help  gain  insight  on  the  effect  of  one  parameter, 
unit  Reynolds  number,  that  has  been  observed  by  many  to 
influence  transition. 

(U)  Potter  and  Whitfield  (Refs.  2,  3),  in  a  series  of 
wind-tunnel  investigations,  found  that  the  transition 
Reynolds  numbers  that  they  measured  on  6-degree  and  10- 
degree  cones  changed  by  more  than  a  factor  of  two,  with  a 
change  of  approximately  an  order  of  magnitude  in  the  unit 
Reynolds  number.  Further,  it  was  noted  that  this  variation 
was  dependent  on  nose  bluntness.  Unit  Reynolds  number 
effects  were  also  reported  by  Pate  and  Brillhart  (Ref.  4). 
Their  results  were  obtained  on  sweptwing  plan  forms. 

Van  Driest  and  Blumer  (Refs.  5*  6)  and  Jack,  et  al  (Ref.  7), 
show  increases  in  the  transition  Reynolds  number  for  a 
10-degree  cone  with  increases  in  the  unit  Reynolds  number 
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while,  in  an  earlier  investigation.  Van  Driest  and  Boison 
(Ref.  8)  found  little  or  no  effect  under  certain  test 
conditions. 

(U)  While  the  majority  of  the  above  wind-tunnel  work, 
along  with  a  number  of  others,  strongly  indicated  the 
existence  of  a  unit  Reynolds  number  effect,  many  remained 
skeptical.  They  pointed  to  the  strong  dependence  of  the 
turbulence  level  in  the  wind  tunnel  on  the  unit  Reynolds 
number  and  suggested  that  this  was  the  real  mechanism  by 
which  transition  was  affected.  Pate  and  Schueler  (Ref.  9)j 
pursuing  this  line,  have  published  a  correlation  that  col¬ 
lapses  data  obtained  from  a  variety  of  wind  tunnels  over  a 
wide  range  in  Mach  number  and  unit  Reynolds  number  into  a 
single  curve.  The  data  are  correlated  in  terms  of  the 
tunnel  diameter  and  the  mean  turbulent  0™  and  5*  of  the 

tunnel  wall  boundary  layer  at  the  test  section.  However, 
tests  have  been  performed  by  Potter  (Ref.  10)  in  a  ballistics 
range  that  cannot  be  explained  by  an  analysis  such  as  that 
of  Pate  and  Schueler,  since  the  ballistics  range  does  not 
have  a  wall  boundary  layer  that  can  radiate  noise  into  the 
test  section.  Potter's  tests  were  made  on  10-degree  sharp 
cones  at  a  local  Mach  number  of  approximately  4.3.  He 
observed  an  increase  in  the  transition  Reynolds  number  by 
a  factor  of  approximately  four,  for  an  order  of  magnitude 
increase  in  the  unit  Reynolds  number.  The  data  obtained  in 
the  present  paper  were  obtained  on  5-degree  sharp  cones  at 
a  higher  local  Mach  number,  approximately  6.9,  and  suggest 
a  somewhat  smaller  effect  at  these  conditions. 


DESCRIPTION  OF  TESTS 

(U)  The  tests  were  conducted  in  the  NOL  Pressurized 
Ballistics  Range  (a  3-foot-diameter  tube,  270  feet  long) 
and  in  the  NOL  1000-foot  Hyperballistics  Range  (Ref.  11) 

(a  10-foot-diameter  tube,  1000  feet  long).  The  ranges  are 
eauipped  with  27  and  37  pairs  of  shadowgraph  stations, 
respectively.  In  addition,  they  have  X-ray  stations  at 
various  positions  along  the  range  to  monitor  changes  in 
model  contour  that  may  occur  during  the  flight  due  to 
aerodynamic  melting  or  damage  during  launch.  The  test 
models  were  5-degree  half-angle  cones  with  a  tip  radius  of 
approximately  0.001  inch.  The  base  diameter  varied  from 
0.5  to  1.75  inches.  The  model  tips  were  machined  from  a 
tantalum  alloy  (TalOW)  and  screwed  onto  a  titanium  after¬ 
body.  The  outer  surface  was  finished  by  grinding,  with  a 
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resultant  surface  finish  of  less  than  10  microinches  rms. 

In  order  to  make  the  models  more  stable  and  reduce  the 
effects  of  angle  of  attack  on  transition,  the  aft  portion 
of  the  models  was  hollowed  out  and  an  internal  ballast  was 
placed  in  the  forward  section.  In  addition,  a  finned 
cylindrical  afterbody  was  placed  on  the  smaller-diameter 
models  to  heln  increase  their  stability.  A  photograph  of 
the  models  is  shown  in  Pig.  1. 

(U)  The  location  of  transition  was  determined  opti¬ 
cally  from  the  range  shadowgraphs.  A  shadowgraph  of  a 
model  in  free  flight  is  shown  in  Pig.  2.  It  can  be  seen 
that  the  range  photographs  have  sufficient  clarity  and 
detail  at  these  test  conditions  to  allow  accurate  estimates 
of  the  location  of  boundary- layer  transition  to  be  made. 

The  transition  Reynolds  number  was  obtained  by  measuring 
the  length  of  laminar  flow  and  determining  the  flow  proper¬ 
ties,  such  as  velocity,  temperature,  and  pressure,  from  test 
conditions  and  trajectory  measurements  in  the  range. 

(U)  The  unit  Reynolds  number  was  adjusted  by  controlling 
the  ambient  pressure  in  the  range.  The  higher  unit  Reynolds 
numbers  were  obtained  in  the  Pressurized  Ballistics  Range. 
This  facility  can  be  pressurized  to  approximately  five  atmo¬ 
spheres.  Tests  were  conducted  in  this  facility  at  range 
pressures  as  high  as  two  atmospheres.  At  these  pressures, 
transition  was  occurring  approximately  one  inch  from  the 
model  tip.  It  was  decided  not  to  investigate  higher  unit 
Reynolds  numbers  due  to  the  difficulty  in  determining 
shorter  lengths  of  laminar  flow  from  the  range  shadowgraphs. 
The  smaller-diameter,  thus  shorter-length,  models  were  used 
for  the  higher  unit  Reynolds  number  tests.  The  lower 
Reynolds  number  tests  were  conducted  in  the  1000-foot  Hyper¬ 
ballistics  Range.  The  launchers  that  are  currently  available 
for  use  with  this  facility  made  it  possible  to  launch  con¬ 
siderably  larger  models,  thus  providing  an  opportunity  to 
measure  longer  lengths  of  laminar  runs.  To  date,  at  lower 
Mach  numbers,  cones  with  mechanically  sharp  tips  have  been 
tested  in  the  ballistics  ranges  at  NOL  with  no  evidence  of 
melting  due  to  aerodynamic  heating  (Ref.  12).  However,  at 
the  present  test  Mach  number  of  nearly  8,  the  tips  were 
purposely  blunted  to  a  radius  of  approximately  0.001.  This 
was  sufficiently  blunt  to  eliminate  the  heating  problem, 
but  still  allowed  the  local  flow  properties  to  recover 
rapidly  to  sharp  cone  properties.  Fig.  3  shows  the  cal¬ 
culated  Mach  number  and  Reynolds  number  distribution  over 
one  of  the  test  models.  It  can  be  seen  that  the  bluntness 
effects  due  to  the  finite  tip  radius  are  quickly  "swallowed," 
reaching  99  percent  of  cone  values  within  the  first  0.4  inch 
along  the  surface.  The  local  flow  properties  were  calculated 
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by  a  momentum-integral  method  described  by  Wilson  (Ref.  13). 
This  method  takes  into  account  the  curved  bow  shock  wave 
that  exists  for  slightly  blunted  slender  bodies,  and  allows 
a  variation  in  total  pressure  along  the  outer  edge  of  the 
boundary  layer  in  the  conical  portion  of  the  body.  It 
assumes,  however,  that  the  static  pressure  along  the  sur¬ 
face  is  constant  and  equal  to  the  inviscid,  sharo-nosed 
cone  value. 


DISCUSSION  OF  RESULTS 


(U)  Since  the  purpose  of  the  present  test  is  to  study 
the  effect  of  unit  Reynolds  number  on  boundary-layer  tran¬ 
sition,  considerable  attention  was  given  to  the  other 
parameters  that  are  known  to  also  affect  transition  to 
minimize  or  eliminate  their  contribution.  The  tests  were 
planned  for  a  Mach  number  of  8.  However,  the  average  test 
Mach  number  for  five  of  the  tests  was  7.o  with  a  variation 
from  7.3  to  7.8.  One  additional  data  point  was  included 
at  a  Mach  number  of  8.3.  Since  the  ambient  temperature  in 
the  range  was  5^-0°R  for  all  of  the  tests,  and  the  models 
were  not  preheated,  the  heat-transfer  rate  varied  with  the 
Mach  number.  To  eliminate  model  geometry  effects,  all 
models  were  made  with  the  same  cone  angle  and  tip  radius. 
The  variation  in  length  was  not  significant,  since  transi¬ 
tion  was  never  measured  near  the  base.  The  surface  of  the 
models  was  ground  to  finishes  of  better  than  10  micro¬ 
inches  rms.  The  maximum  size  of  a  discrete  roughness  pro¬ 
tuberance  found  on  the  cones  was  approximately  120  micro¬ 
inches.  An  analysis  (Ref.  12)  has  shown  that  this  is  not 
sufficient  to  affect  transition  at  the  present  test  condi¬ 
tions  . 

(U)  One  of  the  most  undesirable  elements  in  a  ballis¬ 
tics  range  transition  urogram  is  control  of  the  angle  of 
attack  of  the  models.  As  mentioned  earlier,  the  models 
were  mass  stabilized  by  hollowing  out  the  base  and  using 
internal  ballast  in  the  nose  section,  and  in  some  cones  by 
adding  finned  cylindrical  afterbodies.  While  in  general 
this  was  sufficient  to  keep  the  models  at  small  angles  of 
attack,  occasionally  the  angles  did  get  sufficiently  large 
to  affect  drastically  the  location  of  boundary-layer  tran¬ 
sition.  To  limit  the  effect  of  angle  of  attack  on  the 
present  tests,  only  data  obtained  at  angles  of  less  than 
1.5  degrees  are  presented.  In  general,  the  angle  is  less 
than  1.0  degree. 

(U)  At  combinations  of  high  Mach  number  and  high  range 
pressures,  the  aerodynamic  heating  rates  in  a  ballistics 
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range  become  quite  high.  If  the  heating  rate  becomes  suf¬ 
ficiently  high,  it  can  cause  portions  of  the  model  to  melt, 
burn,  or  ablate,  depending  upon  the  model  material.  The 
problem  becomes  more  acute  for  configurations  with  sharp 
tips,  such  as  those  used  for  the  present  tests.  The  result¬ 
ing  effect  of  a  melting  tip  on  boundary-layer  transition 
would  be  difficult,  if  not  presently  impossible,  to  predict. 
Fortunately,  both  the  Pressurized  Ballistics  Range  and  the 
1000-foot  Hyperballistics  Range  have  in  excess  of  25  pairs 
of  shadowgraph  and  schlieren  stations  located  along  the 
length  of  the  ranges.  By  observing  the  results  of  all  of 
the  data  stations,  any  change  in  the  transition  location  as 
the  temperature  of  the  tip  region  increases  can  be  detected. 
If  a  considerable  change  in  transition  location  occurred 
during  the  flight,  only  data  obtained  prior  to  the  change 
are  used. 

(U)  Transition  Reynolds  numbers  were  obtained  for  each 
flight  by  using  data  from  as  much  of  the  flight  as  nossible. 
Approximately  27  orthogonal  pairs  of  shadowgraphs  are 
obtained  for  each  flight  in  the  Pressurized  Ballistics 
Range,  and  37  pairs  from  flights  in  the  1000-foot  Hyper¬ 
ballistics  Range.  Since  the  range  pressure  and  temperature 
are  uniform  along  the  entire  length  of  the  ranges,  and  the 
change  in  velocity  of  the  model  is  relatively  small  due  to 
the  high  ballistics  coefficient,  the  variation  in  Mach  num¬ 
ber  and  unit  Reynolds  number  during  a  flight  is  small. 
Therefore,  to  obtain  a  transition  Reynolds  number  for  a 
flight,  all  low  angl e-of-attack  readings  of  transition 
location  are  averaged  until  any  evidence  of  tip  melting  is 
observed.  Subseauent  data  are  not  used. 

(U)  To  date,  data  have  been  obtained  from  six  launch¬ 
ings  at  a  nominal  Mach  number  of  7.5  over  a  range  in  local 
unit  Reynolds  number  of  approximately  1  x  10b  to  1.5  x  in' 
per  inch.  A  summary  of  the  test  conditions  and  transition 
results  are  shown  in  Table  1.  The  data  are  also  shown  in 
Fig.  4.  It  can  be  seen  that  the  data  strongly  suggest  that 
there  is  a  stabilizing  effect  of  increasing  the  unit 
Reynolds  number.  In  particular,  it  appears  that,  as  the 
unit  Reynolds  number  is  increased  by  a  factor  of  ten,  the 
transition  Reynolds  number  increases  by  a  factor  of  approxi¬ 
mately  1.7.  The  oresent  data  are  also  compared  with  data 
collected  by  a  number  of  other  investigators  in  Fig.  5- 
While  there  is  a  considerable  amount  of  variation  in  the 
trends  shown  in  Fig.  5,  in  general, the  present  data  show 
less  of  an  effect  than  the  majority  of  the  other  data  shown. 
It  again  should  be  pointed  out  that  the  present  data  were 
obtained  in  a  ballistics  range  and  are,  therefore,  immune 
from  unit  Reynolds  number  sensitive  "tunnel  flow  noise"  that 
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is  found  in  the  wind-tunnel  investigations.  However,  the 
other  range  data,  that  of  Potter  (Ref.  10)  shown  in  Fig.  5, 
also  show  a  much  stronger  effect  of  unit  Reynolds  number 
than  the  present  data. 

(U)  Caution  should  be  taken  in  comparing  transition 
data  from  various  experiments  at  different  conditions.  An 
earlier,  comprehensive  parametric  study  in  the  NOL  ballis¬ 
tics  ranges  of  the  effects  of  Mach  number,  temperature  ratio, 
and  body  geometry  on  boundary-layer  transition  shows  that 
they  can  have  a  very  strong  influence.  Also,  the  effect  of 
these  parameters  can  vary,  depending  upon  the  magnitude  of 
the  various  parameters. 

(U)  Additional  tests  have  been  planned  for  the  NOL 
1000-foot  Hyperballistics  Range  to  extend  the  present  tests 
to  considerably  lower  unit  Reynolds  numbers.  The  tests  will 
be  made  with  models  4  inches  in  diameter,  20  inches  long. 

If  the  trend  observed  in  the  present  tests  continues  at  the 
lower  unit  Reynolds  numbers,  it  will  be  possible  to  measure 
transition  at  unit  Reynolds  number  as  low  as  0.25  x  10°  per 
inch.  This  is  well  within  the  range  of  test  conditions  that 
can  be  obtained  in  wind  tunnels  and  would  provide  an  oppor¬ 
tunity  for  a  direct  comparison  of  values  measured  in  wind 
tunnels  and  ballistics  ranges. 
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also  show  a  much  stronger  effect  of  unit  Reynolds  number 
than  the  present  data. 

(U)  Caution  should  be  taken  in  comparing  transition 
data  from  various  experiments  at  different  conditions.  An 
earlier,  comprehensive  parametric  study  in  the  NOL  ballis¬ 
tics  ranges  of  the  effects  of  Mach  number,  temperature  ratio, 
and  body  geometry  on  boundary-layer  transition  shows  that 
they  can  have  a  very  strong  influence.  Also,  the  effect  of 
these  parameters  can  vary,  depending  upon  the  magnitude  of 
the  various  parameters. 

(U)  Additional  tests  have  been  planned  for  the  NOL 
1000-foot  Hyperballistics  Range  to  extend  the  present  tests 
to  considerably  lower  unit  Reynolds  numbers.  The  tests  will 
be  made  with  models  4  inches  in  diameter,  20  inches  long. 

If  the  trend  observed  in  the  present  tests  continues  at  the 
lower  unit  Reynolds  numbers,  it  will  be  possible  to  measure 
transition  at  unit  Reynolds  number  as  low  as  0.25  x  10°  per 
inch.  This  is  well  within  the  range  of  test  conditions  that 
can  be  obtained  in  wind  tunnels  and  would  provide  an  oppor¬ 
tunity  for  a  direct  comparison  of  values  measured  in  wind 
tunnels  and  ballistics  ranges. 
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FIG.  1.  Range  Models 
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FIG.  Z.  Shadowgraph  of  Model  in  Flight 
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FIG.  3  DISTRIBUTION  OF  LOCAL  PROPERTIES  ALONG  MODEL 
SURFACE. 
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FIG.  4  EFFECT  OF  UNIT  REYNOLDS  NUMBER  ON 
CURRENT  TESTS. 
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Paper  No.  8 

CALCULATION  OF  BLUNT  BODY  FLOWS  USING  PADE  FRACTIONS 
AND  THE  METHOD  OF  CHARACTERISTICS 
(U) 

(Paper  UNCLASSIFIED) 
by 

Andrew  H.  Van  Tuyl 
Naval  Ordnance  Laboratory 
White  Oak,  Silver  Spring,  Md.  20910 


ABSTRACT.  (U)  A  procedure  is  given  for  calculating  the  axially 
symmetric  flow  of  a  perfect  gas  past  a  blunt  body  of  revolution,  using 
a  method  involving  Pade  fractions  in  the  subsonic  region,  and  the  method 
of  characteristics  in  the  suDersonic  region.  The  method  of  calculation 
in  the  subsonic  region  uses  the  Taylor  expansion  of  the  stream  function 
in  the  neighborhood  of  the  nose  of  the  shock.  It  is  an  extension  of  one 
published  earlier,  generalized  to  include  calculation  of  an  arbitrary 
number  of  terms  of  the  Taylor  exoansion.  Calculations  are  made  in  sev¬ 
eral  special  cases,  using  a  method  of  characteristics  program  written 
for  the  IBM  7090  by  R.  H.  Thompson  at  the  Naval  Ship  Research  and  Devel¬ 
opment  Center.  The  examples  comDuted  include  sphere-cones,  sphere- 
cylinders,  spheroid-cones,  and  spheroid-cylinders.  The  extent  to  which 
the  accuracy  in  the  supersonic  region  is  affected  by  the  number  of  terms 
of  the  Taylor  expansion  used  is  investigated. 
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INTRODUCTION 


(U)  A  number  of  methods  for  calculating  the  axially  symmetric 
supersonic  flow  past  a  blunt  body  of  revolution  have  appeared  in  the 
literature  (Refs.  1  and  2,  for  example).  These  include  inverse  methods, 
in  which  the  bow  shock  is  given  and  the  body  which  would  produce  it  is 
calculated,  and  direct  methods,  in  which  the  body  is  given.  When  the 
nose  of  the  body  is  given  by  a  single  analytic  equation  throughout  the 
subsonic  region  and  beyond,  it  is  often  possible  to  solve  the  direct 
problem  in  a  region  containing  the  sonic  line  by  iteration  of  the 
solution  of  the  inverse  problem.  If  the  entire  flow  over  a  given  body 
is  desired,  these  methods  must  usually  be  supplemented  by  another  method, 
such  as  the  method  of  characteristics.  Computer  programs  for  calculating 
axially  symmetric  flow  by  the  method  of  characteristics  are  readily 
available,  as  described,  for  example,  in  Refs.  3  and  4. 

(U)  Procedures  for  finding  the  Taylor  expansion  of  the  stream 
function  in  the  neighborhood  of  the  nose  of  the  shock  in  the  case  of  a 
perfect  gas  have  been  given  by  Lin  and  Shen  (Ref.  5)  and  Cabannes  (Refs. 

6  and  7).  Calculations  by  Van  Dyke  (Refs.  2  and  8)  indicate  that  this 
Taylor  expansion  diverges  at  the  body  for  all  free  stream  Mach  numbers, 
and  therefore  cannot  be  used  directly  to  calculate  the  flow.  As  shown 
in  Ref.  9,  the  use  of  Pade  fractions  obtained  from  the  Taylor  expansion 
is  a  means  of  obtaining  convergence.  Other  methods  for  removing  the 
divergence  of  the  Taylor  expansion  at  the  body  have  been  given  by  Lewis 
(Ref.  10^  Leavitt  (Ref.  11),  and  Sanematsu  and  Chapkis  (Ref.  12). 

Moran  (Ref.  13)  also  uses  Pade  fractions,  but  his  method  differs  con- 
siderably  from  the  present  one. 

(U)  In  Ref.  14,  a  method  for  calculating  the  axially  symmetric  flow 
of  a  perfect  gas  past  a  blunt  body  is  given  in  which  terms  of  the  Taylor 
expansion  up  to  and  including  degree  8  are  used.  In  the  present  paper, 
the  method  of  Ref.  14  is  generalized  to  include  calculation  of  an 
arbitrary  number  of  terms  of  the  Taylor  expansion  of  the  stream  function, 
and  a  procedure  is  given  for  calculation  of  initial  values  for  use  with 
the  method  of  characteristics.  Calculations  are  carried  out  in  several 
special  cases,  using  a  method  of  characteristics  program  written  for  the 
IBM  7090  by  R.  H.  Thompson  at  the  Naval  Ship  Research  and  Development 
Center  (Ref.  4).  The  examples  chosen  include  sphere-cones,  sphere- 
cylinders,  spheroid-cones,  and  spheroid-cylinders.  The  number  of  terms 
of  the  Taylor  expansion  used  in  these  calculations  is  varied  in  order 
to  investigate  how  the  downstream  accuracy  is  affected.  It  is  of 
interest  to  determine  how  many  terms  of  the  Taylor  expansion  are  required 
for  reasonable  accuracy,  since  the  computing  time  increases  rapidly  as  the 
number  of  terms  increases. 
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(U)  As  in  Ref.  14,  let  x  and  r  be  cylindrical  coordinates  with 
origin  at  the  nose  of  the  body.  We  will  consider  a  perfect  gas  with 
ratio  of  specific  heats  y  The  flow  ahead  of  the  shock  is  assumed 
uniform  and  parallel  to  the  x-axis,  with  Mach  number  M  ,  density  p  , 
velocity  magnitude  q^,  and  pressure  p^  =  M^y.  ^et  xq  be  the  shock 

detachment  distance,  and  let  R  and  R  be  the  radii  of  curvature  of  the 
shock  and  body  noses,  respectively.  Then  for  small  r,  we  will  assume 
that  the  equations  of  the  shock  and  body  have  expansions  of  the  forms 


respectively.  In  agreement  with  existing  notation,  the  coefficients 

and  are  called  the  bluntnesses  of  the  shock  and  body,  respectively. 
We  will  assume  that  Rg  and  R^  are  finite  and  non-zero. 

(U)  The  method  of  calculation  to  be  described  is  a  direct  general-- 
ization  of  that  of  Ref.  14  to  an  arbitrary  number  of  terms  of  the  Taylor 
expansion  of  the  stream  function.  The  Taylor  expansions  of  the  stream 
function  and  density  are  first  obtained  in  cylindrical  coordinates, 
starting  from  Bernoulli!s  equation  and  the  vorticity  equation,  as  in 
Ref.  5.  The  Taylor  expansion  of  the  stream  function  is  then  transform¬ 
ed  to  a  new  system  of  orthogonal  coordinates  x  and  £,  as  shown  in  Fig.  1, 
in  which  the  shock  is  the  coordinate  surface  x  =  0.  Rational  expressions 
are  obtained  for  the  flow  on  the  body  and  in  the  shock  layer,  as  in 
Ref.  14.  Calculations  are  carried  out  using  subroutines  for  manipulation 
of  power  series.  The  subroutines  use^.  for  raisin?  a  single  or  double 
power  series  to  an  arbitrary  power  are  based  on  algorithms  due  to  Leavitt 
(Ref.  15). 

(U)  In  the  direct  problem,  a  shock  is  found  by  iteration  from  a  two- 
parameter  family  of  rational  shock  equations.  The  flow  alon?  the  body 
is  obtained,  as  a  function  of  the  arc  length  in  the  inverse  problem,  and 
the  same  variation  with  respect  to  arc  length  is  assumed  to  hold  along 
the  given  body.  Thus,  the  calculated  body  in  the  solution  of  the  in¬ 
verse  problem  is  never  used  directly  in  the  calculation  of  the  flow, 
but  only  as  a  means  of  deciding  when  a  satisfactory  shock  has  been 
found . 

(U)  Finally,  initial  values  for  use  with  the  method  of  character¬ 
istics  are  calculated  along  a  suitable  curve  of  the  form  £  =  constant. 
The  curve  £  =  constant,  and  the  flow  quantities  along  it,  are  approxi¬ 
mated  by  means  of  cubics  in  x  interpolated  through  four  values. 
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SOLUTION  OF  THE  INVERSE  PROBLEM 


THE  TAYLOR  EXPANSIONS  IN  CYLINDRICAL  COORDINATES 


(U)  If  we  take  p  =  q  =  1  and  make  the  substitution  r  =  z 

oo  oo 

Bernoulli’s  equation  and  the  vorticity  equation  become 

2 


1/2 
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+  4T  2  +  F(T)p 

z  z  y-I 
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\  z  zz 


p  Y 

X  X 


+  p 


2Cp2  =  o 
Y+2  f'(?) 
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=  0 


respectively,  where 


1 


C  =  ir  + 


2  (y-l)  m2 

CD 

and  F  (T)  depends  on  the  shock  shape.  As  in  Ref.  5,  we  have 

F(t)  '  sia2°  -  £l)  Al  +  (Y+l)ki„2a 


then 


(3) 


(4) 


(5) 


(6) 


along  the  shock,  where  a  is  the  angle  between  the  shock  and  the  free 
stream.  Denoting  the  equation  of  the  shock  by  z  =  z(x) ,  we  have 

z"2  (x) 


sin2o  = 


4z(x)  +  z  (x) 

The  values  of  T  and  p  on  the  shock  are  given  by 


(7) 


and 


V  =  “  Z 


}  +  2 _ 

P  l  y+l  (y+l)M2sin2c 


(8) 


(9) 


In  order  to  find  the  Taylor  expansions  of  T  and  p  when  the  expansion  of 
equation  (1)  is  given,  it  is  convenient  to  set  R  =  !•  We  can  then 
find  the  Taylor  expansions  in  the  forms 
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respectively.  The  first  step  in  the  calculation  of  the  coefficients  is 

to  expand  the  right-hand  sides  of  equations  (6)  and  (9)  in  powers  of  V 

and  x  +  respectively.  The  former  expansion  is  obtained  from  (6), 

(7),  and  (8),  and  the  latter,  from  (7)  and  (9).  On  substituting  (10) 

and  (11)  in  (3)  and  (4)  and  collecting  coefficients  of  tX (x  +  x^)3,  we 

obtain  algebraic  relations  A—  =  0  and  B^_.  =  0,  respectively,  between 

the  coefficients  of  (10)  and  (11).  The  A^-j  and  are  functions  of  ¥ 

and  p  for  values  of  in  £  i  and  n  i  .  Similarly,  substituting  (1)  ana 

(10)  in  (8), and  (1)  and  (11)  in  (9)  and  collecting  coefficients  of 

(x  +  Xq)^,  we  obtain  algebraic  relations  =  0  and  =  0.  As  before, 

C-r  andUDj  are  functions  of  ¥  and  for  values  of  in  <  i  and  n  <  j  . 

These  four  sets  of  relations  are  sufficient  to  determine  all  the 

and  p . . . 
ij 

(U)  The  coefficients  ¥  ^  p^  Q,  and  pQ  ^  are  found  by  solv¬ 
ing  the  four  simultaneous  eqdation£  A^  5=  0,  A^  0,  -  0,  and 

=  0.  The  determination  of  these  co£rf icients ’ is  special,  since  it 
involves  the  choice  of  the  root  of  a  quadratic  equation.  For  k  >  1,  the 
equations  A^  q  =  0,  ]_  =  0,  =  0,  and  =  0  are  four  simultaneous 

linear  equations  for  the  coefficients  q,  ±  o>  an^  pk-l  1* 

All  other  coefficients  and  p^_.  such  that  2i  +  j  =  Ik  -  1  or  2k  are 

then  found  by  recursion. 

(U)  All  the  substitutions  iust  described  are  carried  out  using 
subroutines  for  power  series  manipulations.  When  K-l  terms  of  the  series 
for  the  shock  equation  are  given,  we  can  calculate  the  and  p  for 
0  <  i  <  K  -  1  and  0  <  j  _<  2K  -  2i  -  2 .  In  the  following,  the  number  of 

terms  used  will  be  specified  by  the  preceding  integer  K.  Thus,  for  a 
given  value  of  K,  K-l  terms  of  the  series  in  equation  (1)  are  used. 

(U)  Finally,  we  can  rewrite  (10)  and  (11)  in  terms  of  arbitrary 
values  of  the  reference  quantities  in  the  forms 


and 


(12) 


(13) 


respectively . 


A  COORDINATE  TRANSFORMATION 

(U)  As  in  Ref.  14,  let  the  equation  of  the  shock  curve  be  given 
by 


219 


8th  Navy  Symposium  on  Aeroballistics 


Vol.  1 


*  +  *c 

R 


=  f 


(14) 


Then  a  transformation  to  coordinates  x  and  E,  will  be  defined  by  the 
relation 


x  +  x 


R 


-  +  i  ^  =  u(c). 


where 


U(C)  =  5  +  f(-?2) 


(15) 


(16) 


and  ?  “  t  +  i£.  This  transformation  is  orthogonal  and  isometric,  and 
is  such  that  the  shock  coincides  with  the  coordinate  surface  x  =  0. 

(U)  It  follows  from  equation  (15)  that 
00 
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(IS) 


i  =  0 

where  the  superscripts  denote  differentiation. 

(U)  In  order  to  express  the  Taylor  expansion  of  V  in  terms  of  x 
and  we  substitute  (17)  and  (13)  in  (12)  by  means  of  subroutines  for 
power  series  manipulations.  We  will  write  the  resulting  series  as  a 
single  power  series  of  the  form 


p  q  R  * 

00  1oo  g 


=  ¥ 


- 1  h 

i  =  1 


(T)5 


2i 


where 


WT>  ■  2  tht1 

3  -  f> 


(19) 


(20) 


When  K  1  terms  of  the  series  for  the  shock  equation  are  used,  the  series 
for  ¥  (x)  is  known  through  the  term  of  degree  2K-2i. 

(U)  We  note  that  it  is  also  possible  to  write  equation  (20)  as  a 
single  power  series  in  x  with  coefficients  which  are  functions  of  . 
This  is  one  of  two  possibilities  which  are  considered  by  Moran  in 
Ref.  13.  The  preceding  form  was  chosen  in  the  present  paper,  however, 
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because  of  the  fact  that  coefficients  of  the  expansions  of  the  body 
equation  and  flow  quantities  in  power  of  £  are  given  exactly  by  finite 
combinations  of  the  ¥.(x)  and  their  derivatives.  In  the  calculation  of 
the  flow  behind  a  given  shock,  approximations  for  the  T.(x)  and  their 
derivatives  need  to  be  calculated  only  once.  1 

PADE  FRACTIONS 

(U)  Given  a  power  series 


oo 


i 


(21) 


c .  z 
1 


with  cq  ^  0,  the  Pade  fraction  f^  n(z)  ,  k  >_  0,  n  >_  0,  is  defined  as 
follows:  it  is  a  rational  fraction  with  numerator  and  denominator  of 
degrees  less  than  or  equal  to  n  and  k,  respectively ,  such  that  the 
Taylor  expansion  of  f^  n(z)  at  the  origin  agrees  with  equation  (21)  to 
more  terms  than  that  of  any  other  rational  fraction  with  numerator  of 
degree  <  n  and  denominator  of  degree  <_  k  (Ref.  16).  It  is  determined 
by  the  first  n  +  k  +  1  coefficients  of  equation  (21),  and  is  unique. 

When  k  is  not  too  large,  f ^  n (z)  can  be  calculated  conveniently  as  the 
quotient  of  two  determinant^  of  order  k  +  1  or  less  (Ref.  17,  equation 
(75)).  For  larger  values  of  k,  the  OD  algorithm  of  Rutishauser  (Ref.  18) 
is  a  convenient  method  of  calculation. 

(U)  Pade  fractions  with  k  and  n  equal  or  nearly  equal  often  give 
good  results  for  relatively  small  values  of  n,  even  when  the  given  power 
series  does  not  converge.  A  proof  of  convergence  is  not  available  in 
the  present  problem,  but  numerical  results  indicate  that  the  sequences 
of  Padd  fractions  obtained  converge  at  the  body,  while  the  original 
power  series  diverge. 

(U)  Finally,  the  following  property  of  Pade  fractions  is  of  impor¬ 
tance  in  the  calculation  of  the  flow  in  the  shock  layer:  Let  f (z)  be 
a  rational  function  reduced  to  lowest  terms  with  numerator  of  degree  q 
and  denominator  of  degree  p.  Then  f ^  ^(z)  is  identical  with  f(z)  when 
k  >_  p  and  n  _>  q.  This  result  is  discussed  in  Ref.  17,  page  22. 

RATIONAL  APPROXIMATIONS  FOR  THE  7(t)  AND  THEIR  DERIVATIVES 

(U)  We  will  obtain  rational  approximations  for  the  T.(x)  and  their 
derivatives _by  forming  Padd  fractions  from  the  corresponding  partial 
sums.  The  ¥.(x)  and  their  even  derivatives  are  replaced  by  rational 
functions  with  numerator  and  denominator  of  equal  degrees,  and  all  odd 
derivatives  are  replaced  by  fractions  in  which  the  degree  of  the 
denominator  exceeds  that  of  the  numerator  by  one. 

(U)  The  accuracy  of  these  rational  approximations  is  greatest  at 
the  shock,  and  decreases  as  the  body  is  approached.  When  K  -  1  terms 
of  the  series  for  the  shock  are  given,  the  rational  approximations  for 
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the  Ti(x)  are  sufficiently  accurate  at  the  body  only  for_  i  less  than 
about  K/2.  Similarly,  the  order  of  the  derivatives  of  ^(t)  must  be 
less  than  about  K/2  -  i  at  the  body. 


CALCULATION  OF  THE  BODY 

(TJ )  In  coordinates  t  and  £,  the  equation  of  the  body  is  given  by  a 
series  of  the  form 


t  -  tq  =  a 1?2  +  a2 ^  +•••,  (22) 

where  Tq  is  the  smallest  positive  root  of  the  eauation  =  0.  When 

K  -  1  terms  of  the  serie_s  for  the  shock  equation  are  given,  the  numera¬ 
tor  and  denominator  of  ^(r)  are  of  degree  K  -  1.  We  obtain  the  co- 
efficients  in  equation  (22)  successively  by  substituting  (22)  in  (19) , 
collecting  like  nowers  of  £2,  and  equating  the  coefficients  to  zero. 
Substituting  (22)  in  (17)  and  (IS)  and  collecting  like  nowers  of  ^  we 
obtain  x/Rs  and  r/Rs  along  the  body  as  power  series  in  together  with 
values  of  Xq/R-^  and  R^/Rs.  Finally,  eliminating  E,  between  these  equa¬ 
tions,  we  obtain  the  equation  of  the  body  in  the  form  shown  in  equation 
(2). 

_(U)  The  coefficients  of  equation  (2)  calculated  in  this  way  involve 
the  ^(t)  and  their  derivatives  evaluated  for  t  =  Tq.  When  K  -  1  terms 
of  equation  (1)  are  given,  we  can  calculate  only  about  K/2  terms  of 
equation  (2)  with  sufficient  accuracy. 


(U)  We  obtain  increased  accuracy  in  the  calculation  of  the  body  by 
replacing  the  right  hand  side  of  equation  (2)  by  a  Padd  fraction.  "When 
K  =  9,  10,  and  11,  the  rational  approximation  found  in  this  way  is  of 
the  form 


1  +  b1(x/P^)  +  \)2  (x/R^)  2 
1  +  b^(x/Rb)  +  b^(x/Rb)2 


(23) 


(U)  We  also  need  an  expansion  for  £  in  powers  of  the  arc  length  s 

along  the  body.  From  equation  (18)  and  the  expansion  of  the  arc  length 

in  Dox^ers  of  r,  and  using  the  calculated  value  of  R,  /R  .  we  obtain 

b  s 


When  K  =  11,  we  find  a  rational  approximation  of  the  form 
r2  1  h  2  /S_ \2  1  +  HL(s/Rh)2  +  H2(s/Rbp 
1  (Rs  J  1  H-  H3(s/FLb)2  +  H4(s/Rb)4 


(24) 


(25) 


CALCULATION  OF  THE  FLOW  ON  THE  BODY 

(U)  We  obtain  the  density  on  the  body  from  the  equation 
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(r/Rs)2|u^(?)|2 


=  2C 


'PSt\2 


St 


) 


Y+l 


Vol.  1 


(26) 


where  p  is  the  density  at  the  stagnation  point.  We  have 


st 


(v  -  Pc 

yf(0) 


1 

Y-l 


(27) 


Substituting  the  expansion 


-9-  •  =  1  +  d.£2  +  d_£4  +• 
Q  1  Z 

pst 


(28) 


in  equation  (26)  and  collecting  like  powers  of  ^ ,  we  can  determine  the 
coefficients  of  equation  (26)  successively.  Finally,  from  (24)  and  (28), 
we  have 


=  1  + 


St  “  \  \ 


+  "2  K 


+• 


(29) 


We  obtain  expansion  for  the  pressure  and  velocity  magnitude  on  the  body 
by  substituting  (28)  in  the  equations 


Y 


and 


st 


Si.  - 


st 


(30) 


2C 


(31) 


respectively,  where 


P  .  =  p  q2  F(0)  (  -Si.  ]  Y 

■  St  ^oo-oo  W\  / 


is  the  pressure  at  the  stagnation  point. 


(32) 


(U)  As  in  the  calculation  of  the  body,  it  is  found  that  convergence 
is  improved  when  the  preceding  partial  sums  are  replaced  by  Padd  frac¬ 
tions.  When  K  =  11,  we  have 


a_  _ 


i  (S/V  ■  2 


+  n  r~/v>  \ 3 


(s/I lb)- 


1  +  03  (s/R^)' 


+  04  (s/R^1 


(33) 
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and 


i  +  p-^s/y2  +  ys/y4 

Pst  1  +  P3(s/Rt)2  +  P4(s/Rb)4 
with  an  approximation  of  the  same  form  as  (34)  for  p/pst- 


(34) 


CALCULATION  OF  THE  FLOW  IN  THE  SHOCK  LAYER 

(u)  Immediately  behind  the  shock,  the  quantities  u,  p,  p,  c  ,  and 
q2 5  where  u  is  the  x-component  of  velocity  and  c  is  the  speed  of  sound, 
are  rational  functions  of  cos2a  such  that  the  degrees  of  the  numerators 
and  denominators  are  less  than  or  equal  to  2  and  1,  respectively  (Ref.  19). 
As  in  Fig.  1,  let  a  at  a  given  point  be  the  angle  between  the  tangent  to 
the  curve  x  =  constant  at  that  point  and  the  positive  x-axis.  Since  a 
is  equal  to  a  at  the  shock,  it  follows  from  the  special  property  of 
Pade  fractions  stated  earlier  that  if  we  expand  the  preceding  flow 
quantities  in  powers  of  cos9a  in  the  shock  layer,  then  all  Pade  frac¬ 
tions  formed  from  these  series  with  k  >_  1  and  n  >_  2  are  exact  at  the 

shock. 

(U)  In  order  to  find  these  expansions,  we  first  obtain  expansions 
in  powers  of  £2.  In  coordinates  x  and  5,  Bernoulli !s  equation  becomes 

?  2  +  ?  2 

x _ § _ 

(r/Rs)2|u'fe ) |2 

Substituting  the  expansi 


+ 


2yFffl 

Y-l 


=  2C 


(35) 


=  R  (x)  +  R0(t)52  +  R.(tH4  +  •••  (36) 

p  1  2  J 

oo 

in  equation  (35),  collecting  like  powers  of  £2,  and  equating  the  coeffi¬ 
cients  to  zero,  we  first  obtain 


R/+1  -  2CE  2  +  ^  ,  o  (37) 

Y-l  1  1  „'4 

For  a  given  value  of  x,  equation  (37)  is  easily  solved  for  R^  by  Newton's 
method.  After  finding  Rj_,  we  calculate  the  remaining  coefficients  in 
(36)  by  recursion. 

(U)  We  obtain  expansions  of  the  same  form  as  equation  (36)  for 
p/p^q^s  (c/q ot):  ,  and  (q/q^)2  by  using  the  relations 


F 

. 

p  q 

OO  00 


(38) 
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and 


t 


3- 

i  . 


I  =  .E/.p.oo,q,o 

'  Y  p/p„ 


+■ 


Y-X 


=  2C 


Vol.  1 

(39) 

(40) 


in  the  order  given.  Finally,  we  obtain  an  expansion  for  u/q  by  expand¬ 
ing 


PU 

P  u 

OO  OO 


-(x/R  )  4*  +  (x/R  ) 

_ S  4  T _ ST  g 

(r/R  ) IlMO I2 


(41) 


in  powers  of  £2  and  dividing  by  (36) 
(U)  We  have 


cos2a  = 


[(x/R  )  J2 


s  g 


u'  (O  I2 


(42) 


from  which  we  find 


=  E^(x)  cos2a  +  E^x)  cost+a  +*  *  * 


(43) 


On  substituting  (43)  in  the  preceding  expansions,  we  obtain  the  desired 
expansions  in  powers  of  cos?a. 

(U)  When  K  =  11,  we  obtain 


P  q 

OO^OO 


P0(T)  +  p^(x)cos2a  +•  •  •  4-  p^(r)cos1(^a 
1  4-  p^(x)cos2a  +• ■ •  4-  p^Q(x)cos10a 


(44) 


with  rational  approximations  of  the  same  form  for  p/p  ,  u/q  ,  (q/q  )2 

-  f  I  \  p  oo  oo  ■L  ^-oo  7 

and  (c/q^;  .  These  approximations  are  exact  at  the  shock,  and  remain 
accurate  throughout  most  of  the  shock  layer.  Near  the  body,  however, 
higher  order  coefficients  in  the  expansions  become  inaccurate,  and 
better  accuracy  may  be  obtained  by  forming  Pade  fractions  from  fewer 
terms  of  the  series.  When  x  =  0  in  these  approximations,  the  fractions 
obtained  are  not  necessarily  in  lowest  terms. 


(U)  We  obtain  expansions  for  (x  4*  Xq) /R  and  r/Rs  in  powers  of 
cos2a  by  substituting  (43)  in  (17)  and  (18),  respectively.  When  K  =  11, 
we  obtain  the  rational  approximations 


x  +  x0 
R 

s 


z  (t)  4-  z^(x)cos2a  +• • • 

1  4*  z, (x) cos2a  +* • * 

0 


+z^ (t) cos1 °a 
+  z^Q(x)cos10a 


(45) 
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and 


r 


(46) 


R 


1  +  y^(T)cos2a  +• • •  +  yg(T)cos10a 


s 


The  coordinate  transformation  depends  only  on  the  shock  equation,  and 
it  is  not  necessary  to  use  the  same  number  of  terms  of  equation  (1)  in 
the  calculation  of  x  and  r  as  in  the  calculation  of  the  flow. 

(U)  Finally,  we  will  need  a  rational  approximation  for  £  in  terms 
of  cosa.  Taking  the  square  root  of  (43)  and  forming  a  Pade  fraction, 
we  obtain 


(47) 


1  +  e^(x)cos  a  +•  •  •  4-  eg(x)cos10a 


Since  £  becomes  equal  to  r/R  when  t  =  0,  it  follows  that  (46)  and  (47) 
are  identical  at  the  shock. 


SOLUTION  OF  THE  DIRECT  PROBLEM 


(U)  In  the  direct  problem,  we  use  the  two-parameter  family  of 
shock  equations  given  by 


(48) 


The  Taylor  expansion  of  the  right-hand  side  of  (48)  at  r  =  0  begins 
with  the  terms 


Comparing  with  equation  (1),  we  see  that  A0  =  B  /2. 

(U)  A  value  of  A3  is  first  chosen,  and  A2  is  found  by  the  method 
of  false  position  so  that  the  bluntness  of  the  calculated  body  matches 
that  of  the  given  one  to  a  prescribed  number  of  figures.  New  values  of 
A^  are  then  chosen  until  the  calculated  body  intersects  the  given  one 
near  some  prescribed  point.  Usually  two  or  three  choices  of  A3  are 
sufficient  to  match  the  given  body  closely  in  the  subsonic  region. 


226 


8th  Navy  Symposium  on  Aeroballistics 


CALCULATION  OF  THE  SONIC  LINE 


Vol.  1 


(U)  Sonic  speed  behind  the  shock  is  given  by 


Replacing  the  left-hand  side  of  (50)  by  its  rational  approximation  in 
terms  of  cos^a,  we  solve  for  cos2a  by  Newton’s  method  for  a  given  value 
of  t.  We  then  find  the  corresponding  values  of  the  cylindrical  coordi¬ 
nates  from  equations  (45)  and  (46) .  Flow  quantities  at  that  point  can 
be  found  from  their  rational  approximations  in  terms  of  cosa. 

(U)  As  in  Ref.  14,  we  find  the  sonic  point  on  the  body  from  the 
relation 


JR 


P 


st 


(51) 


Replacing  the  left-hand  side  by  equation  (34),  we  solve  for  (s/R-^)^  at 
the  sonic  point  by  Newton’s  method.  Finally,  we  interpolate  cubics  in 
x  through  the  values  of  r/R^  and  any  desired  flow  quantities  evaluated 
at  four  values  of  x.  These  values  usually  include  the  points  on  the 
shock  and  body. 


CALCULATION  OF  INITIAL  VALUES  FOR  THE  METHOD  OF  CHARACTERISTICS 

(U)  A  convenient  initial  curve  to  use  for  starting  the  method  of 
characteristics  is  a  coordinate  curve  E,  ~  constant.  For  a  given  value 
of  x ,  we  can  find  the  value  of  cosa  on  a  curve  £  =  constant  by  solving 
equation  (47)  by  Newton’s  method.  The  corresponding  values  of  the 
cylindrical  coordinates  and  flow  quantities  are  then  obtained  from  their 
rational  approximations  in  terms  of  cosa.  We  calculate  the  flow  at  the 
intersection  of  the  initial  curve  with  the  body  by  solving  equation  (25) 
for  (s/R^)2  and  substituting  the  latter  in  the  rational  approximations 
for  the  flow  on  the  body. 

(U)  As  in  the  case  of  the  sonic  line,  the  curve  £  =  constant  and 
flow  quantities  along  it  can  be  approximated  accurately  by  cubics  in  x. 
These  cubics  are  interpolated  through  values  of  r/R-^  and  the  flow 
quantities  at  four  values  of  x  exactly  as  before. 

(U)  In  the  program  described  in  Ref.  4,  the  variables  required 
along  the  initial  curve  are  tan6 ,  i/m2"  -  1,  and  S/cv,  where  0  is  the 
angle  between  the  streamline  and  the  positive  x-axis,  M  is  the  local 
Mach  number,  S  is  the  entropy,  and  cv  is  the  specific  heat  at  constant 
volume.  A  typical  initial  curve  and  characteristic  net  is  shown 
schematically  in  Fig.  2.  In  the  program  of  Ref.  4,  the  calculations 
proceed  along  the  right-running  characteristics  from  top  to  bottom. 
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NUMERICAL  RESULTS  AND  DISCUSSION 


(U)  Calculations  have  been  carried  out  on  the  IBM  7090  for  a 
sphere-cone,  sphere-cylinder,  spheroid-cone,  and  spheroid-cylinder, 
where  the  half  angle  of  the  cone  is  20°  and  the  bluntness  of  the 
spheroid  is  4.  All  calculations  were  made  with  y  =  1.4. 

(U)  In  Fig.  3,  the  shock,  sonic  line,  and  a  typical  initial  curve 
are  shown  for  the  case  of  a  sphere  at  M  =4,  with  K  =  7  and  11. 
Corresponding  results  are  shown  for  a  sphere  at  M  =10  in  Fig.  4  and  for 
a  spheroid  of  bluntness  4  at  M^  =  5.98  in  Fig.  5,°°with  K  =  11.  The 
calculated  bodies  after  the  iterations  have  been  completed  are  also 
shown  in  Figs.  3,  4,  and  5  for  K  =  11. 

(U)  Figs.  6  and  7  show  the  pressures  on  the  sphere-cone  and 
sphere-cylinder  for  K  =  11  calculated  with  the  initial  curves  in  Figs. 

3  and  4,  respectively,  and  Fig.  8  shows  the  pressures  on  the  spheroid- 
cone  and  spheroid-cylinder  calculated  with  the  initial  curve  of  Fig. 

5  and  K  =  11.  The  pressures  in  Figs.  6  and  7  are  compared  with  cal¬ 
culations  by  Chushkin  and  Shulishnina  (Ref.  20).  In  addition,  the 
pressure  on  the  sphere  in  Fig.  6  is  compared  with  measurements  by  Xerikos 
and  Anderson  (Ref.  21)  at  =  3.975.  The  pressure  on  the  spheroid  in 
Fig.  8  is  compared  with  measurements  by  Pasiuk  (Ref.  22). 

(U)  As  seen  in  Fig.  3,  the  calculated  flows  for  K  =  7  and  K  =  11 
differ  from  each  other  mainly  near  the  sonic  line  in  the  vicinity  of 
the  body.  Even  though  the  pressure  found  on  the  sphere  for  K  =  7  is 
about  5%  less  than  that  for  K  =  11  near  the  initial  curve,  the  calculated 
pressures  a  short  distance  downstream  of  the  initial  curve  are  nearly 
identical  with  each  other.  The  calculations  of  Figs.  3  and  6  have  also 
been  made  for  K  =  9,  and  the  results  are  indistinguishable  from  those 
for  K  =  11  on  the  scale  of  the  figures. 

(U)  These  results  indicate  that  reasonable  accuracy  is  obtained  in 
the  calculations  for  the  sphere-cone  and  sphere-cylinder  when  K  =  7,  and 
that  the  results  for  K  =  9  are  nearly  as  accurate  as  those  for  K  -  11. 

The  time  of  computation  for  one  case  of  the  inverse  problem  is  about 
15  seconds  for  K  =  7,  22  seconds  for  K  =  9,  and  50  seconds  for  K  =  11. 
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(U)  FIG.  3  Sonic  Line  and  Initial  Curve  in  Flow  Past  a  Sphere 
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(U)  FIG.  4  Sonic  Line  and  Initial  Curve  in  Flow  Past 

a  Sphere  at  M  =10 
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(U)  FIG.  5  Sonic  Line  and  Initial  Curve  in  Flow  Past  a  Spheroid  of 

Bluntness  4  at  M  =5.98 
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(U)  FIG.  6  PRESSURE  ON  A  SPHERE-CONE  AND  A  SPHERE- 
CYLINDER  AT  /V^  =  4 
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(U)FIG.  7  Pressure  on  a  Sphere-Cone  and  a  Sphere-Cylinder 

at  M  =10 
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(U)  FIG.  8  Pressure  on  a  Spheroid-Cone  and  a  Spheroid-Cylinder 

of  Bluntness  4  at  M  =5.98 
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